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Abstract. 


Iy the first part of this paper the duration of contact during 
collision, of two equal spheres of the same material, is 
investigated with a ballistic pendulum, employing the partial 
discharge of a condenser to measure the duration. For 
aluminium, tin, brass, babbitt, and a lead-tin alloy, a law of 
the form ¢ = ¢)+a/v" is obtained, where v is the relative 
velocity at impact, f, a, and n are constants. Zinc follows 
a similar law, but the constants cannot be accurately 
determined. 

In the second part a description of the flattened areas 
remaining on the spheres after collision is given, and an 
explanation of the rim which surrounds them is offered. 
If d is the diameter of these rims, d = b(v—w)”, where 
b, vp, and m are constants different for different spheres. 
From vp and the elastic moduli of the materials the least 
pressure required to start flow is calculated in the various 
cases, and is found to be of the order 10° dynes/em. 


* Communicated by Prof. C. H. Lees, F.R.S, 
Phil. Mag. 8.7. Vol. 8. No. 53. Suppl. Dec. 1929. 3G 


782 Mr. J. P. Andrews. on the 


Introduction. 


Ree into the process of collision between two 

solid bodies has followed two main routes, originally 
laid down by St. Venant and Hertz respectively. Of these, 
the wave theory of St. Venant has been found to apply to 
long bars, and bodies whose dimensions are comparatively 
large. For bodies whose dimensions in the direction of the 
motion at impact are small, Hertz’s theory of perfectly 
elastic collision is known to be more in harmony with 
experiment, No precise criterion can be given, but it is 
known that where the duration of contact is large compared 
with the time required for a wave of compression to travel 
across the body and back again Hertz’s theory applies *. 
Thus, in the case of aluminium, where the velocity of sound 
is 5‘1 x 10° cm./sec., 2//(5°1 x 10°) has to be small compared 
with the duration of impact, / being the appropriate linear 
dimension. This paper is concerned with experiments upon 
spheres small enough for Hertz’s theory to apply if the 
materials are perfectly elastic. 

Investigations of the impact of perfectly elastic spheres 
have, however, already been carried out, and the formule 
deduced by Hertz for the area of compression and the 
duration of contact have been verified, first by Hertz 
himself and then by other investigators+. If v is the 
relative velocity of two equal spheres of the .same material 
at the moment of collision, and ¢ is the duration of their 
contact, ty/> = constant, according to Hertz. 

When, however, we employ materials other than the 
hardest steel the assumption of perfect elasticity is no longer 
justified, and modifications must be introduced. Assuming 
a resistance to advance proportional to the square of the 
velocity, from the moment the spheres touch to the moment 
they begin to recede, T. Péschl has shown how to modify 
the theory of Hertz, and has produced an approximate 
solution of the problem, which shows that tv = constant in 
this case {. It was the original purpose of this research 
to discover whether or not this law were obeyed. 


Method of Investigation. 


A good method for investigating this problem was already 
to hand, and had been used with various modifications by 


* Vide B. W. Tschudi, Phys. Rev. xviii. p. 423 (1921); J, E. P. Wag- 
staff, Proc. R.S. ev. p. 544 (1924), 

+ A. Dinnik, Jour. d. Russ. Phys. Ges. xli. p. 57 (1909); W. Miiller, 
Weiner Ber. exxiii. p. 2157 (1914), and others. 

T. Pésehl, Zs. fiir Phys, xlvi. p. 142 (1927). 
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a number of investigators. Only minor improvements have 
been incorporated here. A condenser is charged, the impact 
of the metal spheres is allowed to partially discharge it, and 
the remaining charge is immediately afterwards measured 
by passing it through a ballistic galvanometer. From 
a knowledge of the original and final charges and the 
resistance of the circuit, the time during which the spheres 
were in contact may be deduced. This method appears to 


Fig. 1, 


Diagram o Apparatus. 


have originated with R. Sabine* in 1876, and has been used, 
with small differences, by others since, notably by Wagstaff. 
The chief pitfall to be avoided is the inclusion of appreciable 
inductance in the circuit, since the formula for calculating 
the duration of contact then becomes far too complicated. 
The innovations in the present experiments are the employ- 
ment of the pendulum suspensions of the spheres as conducting 
leads and an arrangement enabling the observer to turn the 


* KR. Sabine, Phil Mag. i. p: 337 (1876) 
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balls after each collision so as to present a new pair of 
surfaces at each impact. The spheres were each suspended 
by four steel wires so as to form pendulums about 180 em. 
long, which were so adjusted that the bails were in light 
contact at points on their equators when at rest. The balls 
were screwed, as shown, to small bars of brass, to which the 
suspension wires were soldered. The electrical circuit is 
represented above. The condenser C is charged by depressing 
K on the left, while the spheres are separated. This separation 
is secured by drawing one of the spheres aside by means of 
a thread passing under a glass rod, and fastened to a hook 
depending from a platform P. Depressing K to the right 
then discharges the condenser through the ballistic galvano- 
meter G, producing a throw @). The condenser is charged 
anew, platform P, which is controlled electro-magnetically, 
is removed by sudden rotation in the direction of the arrow 
so as to let the hook and displaced sphere fall perfectly 
freely, and the spheres collide. Immediately after the 
collision, the condenser is discharged through G, a deflexion 0 
being obtained. In order to prevent a further collision, the 
second sphere, in moving back, actuates a very easily moved 
switch, and thus causes a flag of mica covered with soft felt 
to move into the space between the two spheres, insulating 
them. The masses of the two spheres are adjusted to equality 
by adding a small weight to one of them if necessary. 


Theory of the Method and Preliminary Experiments. 


(1) The Duration of Contact. 


If self-inductance is absent, and the total resistance of 
the circuit through which the condenser discharges during 
collision is R, then during this period 


dq 
where g is the charge on the condenser. ‘The solution of this 
equation is 


t 
G=qQMe Re, (1) 


C being the capacity of the condenser. If, however, self- 
inductance is present, but only to a small extent, such that 


ae is small, then the equation of motion, which is-now 


strictly 


tka 0, 
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may be written approximately 


Rat (Gt 


since from (1) 
dt? 


and the new solution is 


L 


In actual practice g was measured by discharging through 
a ballistic galvanometer whose deflexions @ were » proportional 
to the quantity discharged, so that (2) may be written 


log @ = log alc 
Hence, provided self-inductance is zero, and there is no 
other perturbing quantity, z will be a linear function of 
log 0, the time ¢ being kept constant. 

This matter may be tested by separating the balls to the 
same distance on each occasion, so that the velocity of 
approach is always the same, and trying the collision experi- 
ment with different values of Rin the cireuit. But asecond 
consideration enters here. When the spheres first touch, 
the resistance at their junction must be extremely large, then 
fall rapidly, reaching a minimum, and finally increase to an 
infinite value. It may be considered as having an unknown 
average value r, which may be a not altogether negligible 
quantity. Equation (3) ought, therefore, to be ‘written 
approximately 


log 0 = log (4) 
OR(1+ 


From the form of this equation, provided a wide range of 
values of R is considered, the two terms in the brackets cannot 
always compensate one another unless each is separately 
unity. An experiment was arranged on these lines, the 
spheres being rotated through a small angle after each 
collision, using values of R ranging from 113 to 10,000 ohms, 


and a perfectly rectilinear graph was obtained when a was 


log g = log yo— R\CG CR? 
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plotted against log @. It follows that r = 0 and L = 0 very 
closely. The usual value of R was about 500 ohms. (See 


fig. 2.) 
(2) The Velocity of Approach, 


When the bob of a pendulum of length / is withdrawn 
a horizontal distance x, it may be shown that the velocity it 


Fig. 2. 


acquires in falling freely to the lowest point of its path is 


given by 
v= in} 
so long as # is small compared with J. If appreciable 


damping occurs, and 2 is the logarithmic decrement, then 


the above expression must be multiplied by (1— *), giving 
finally 
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The value of «?/4/ never exceeded 1 per cent., and \ varied 
between 0°012 for aluminium and 0:002 for tin. The 
corrections were therefore considered negligible. 


The Initial Separation of the Spheres. 


As it was specially desirable not to contaminate the 
surfaces of the spheres, the following methods were employed 
to measure their initial separation. For small separations 
up to about 20 em.a large lens was arranged to produce 
a diminished real image of the spheres, and the horizontal dis- 
tance between these images was measured in a good travelling 
microscope. It is estimated that by this means the separation, 
and therefore the velocity of approach, could be found to 
about 2 per cent. accuracy, with the exception of the 
smallest velocities, which were sometimes in error owing to 
slight draughts. After each set of observations the readings 
of the microscope were calibrated by observing the apparent 
length of 10 cm. of a metre scale placed in the position 
previously oceupied by the spheres. For larger distances, 
where aberration in the lens may have affected the results, 
the distances were measured directly by a metre scale with 
stiff paper pointers. For the same degree of accuracy it 
was only necessary to be certain of the reading to within 
4 mm., and it was easy to read itto 1 mm. For simplicity, 
all the results in this paper are reduced to the microscope 
scale of the first method. 


The Condenser employed. 


This was a Siemens Standard 1 mfd., selected because of 
its small absorption and leakage. In the course of the 
experiment, as already outlined, the condenser was charged 
about half a second before the collision, and then dis- 
charged about half a second after, so that the chances of any 
spurious effects due to absorption and leakage were reduced 
to entirely negligible quantities. The short-circuiting plug 
of the condenser was inserted for some minutes between 
collisions. The e.m.f. of the charging battery was carefully 
watched, and kept at 4 volts. 


The Galvanometer. 

In the majority of experiments a fluxmeter was employed, 
this instrument having the advantage of a large periodic 
time (46 sec.), which enables the ballistic throw to be read 
off quite accurately. A separate experiment proved that the 
throw was accurately proportional to the charge suddenly 
passed through the coil over the whole range employed. 
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T he Spheres. 


These were first cast and then turned in the Instrument 
Shop of East London College. The equatorial zone was 
turned with special care to a truly spherical figure, since 
here the collision was to take place; and the whole then 
polished. In some instances the spheres were subsequently 
annealed, but this was not possible in the case of zine and 
aluminium, as non-conducting films were thereby formed. 
In all cases the spheres were 4:0 cm, in diameter. 


Part I. 
Results. 


Six substances have so far been examined, viz., aluminium, 
tin, zinc, an alloy of lead and tin containing 60 per cent. 
tin, and a babbitt alloy. Typical results follow. 


TABLE I. 
Tin, Hapt. 2. 


Temperature of the air 10°2 C. 

Resistance R, including that of the steel leads =513°2 ohms- 
Length of the pendulum=184°5 cm, 

Capacity of condenser=1 microfarad. 

10 cm. of metre scale=0°785 cm. on microscope. 

65 = 32°77 em. 


Separation of Throw of t (sec.) 9 v t cale 
spheres, in em. of galv. 0, 1:18 log, =29'408 sec. 
microscope, =8. cm. em,/sec. 

0:050 10 35 5°91 1:47 5°91 
0:120 12:18 5:09 3°53. 5:01 
O177 13:02 £73 5:20 4°71 
0285 12:95 4:76 8:38 4°38 
0°485 14°60 4:15 12°78 4:15 
0-900 15°49 3°85 26°5 
1-420 16°15 3°62 41:8 3°65 
1°875 16°58 3°50 55:0 3°57 
2°63 16-77 3°45 173 3°48 
0:032 8°85 6°71 0°94 6°49 


These results and those in Table I]. are represented graphi- 
cally in fig. 3. Curves representing the results for the other 
materials are also given in figs. 4 and 5. 
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II. 
3. 


Aluminium. 
Temperature of the air 18°C. 
Resistance R, including that of the steel leads, =525 ohms. 
Length of pendulum=182:0 cm. 
Capacity of condenser=1 microfarad, 
10 cm. of metre scale=0°800 on microscope. 
cm. 


Observed ¢ (sec.) 


Separation of Throw of 6 v t eale 
spheres, in cm. of galy. =1-212log., =29058 sec, 
microscope, =S. em. x10-* em./sec. 
11:99 291 2°46 3°28 
0150 11:53 2°65 4:50 2°47 
0°230 13:14 1:96 6°68 2°13 
0°335 13°09 1:99 973 1:90 
0°515 13°66 176 14:96 171 
0-710 14:10 1:60 20°6 162 
0°935 13°89 168 27°2 1:56 
14-73 1:37 34°3 151 
1*425 14°50 1°44 41:3 1:49 
1:700 1461 1°41 49°2 1°47 
10°20 3°30 2°90 3°04 
0:065 8:35 5°34 1°89 3°86 
2-000 14°45 1°46 58°0 1:45 
1:020 14:00 1:60 29°6 1:54 
Fig, 3 
7-0 
Observations 
6-0 X Calculations 
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Fig. 4. 


Babbitt I. 
Zinc 


tx 10*sec. 
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Vv cm/sec. 


Fig. 5. 
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The answer to the original question may now be sought, 
namely, do these results agree with the assumption of 
a resistance to distortion proportional to the square of the 
relative velocity? Do they, that is to say, follow a law of 
the form tv=constant ? 

In no case is this law obeyed. In some instances a law 
of the form (¢—t))v=constant would be a first approximation; 
but after an extensive investigation, the best agreement is 
found to be with a formula of the type 


t=t)+ a/v”, 


where t, a, and 7 are all constants which differ for the 
different materials. The values of these constants are 


TaBLe III. 
Substance. Expt aX 104, n. Remarks. 
2 1°38 20°7 1°58 Points scattered. 
LS 1:34 4°52 0-92 Re-turned. 
Tin 2 2°82 3°58 0°39 
3 3°40 6°88 0°92 Re-turned. 
1:89 3:18 0°32 
Babbitt 2 2°43 661] 0:98 Re-turned. 
1:93 5°88 0°59 Re-turned and 
annealed. 
1 0:16 3:79 0:19 
Lead-Tin ‘ 
} 1 150 500 0:30 
Zi No formula found to cover the whole range, but the 
curve runs in the same general way. 


tabulated in Table III., and the crosses on the curves of 
figs. 3, 4, and 5 are points calculated from the formula to 
show how nearly they agree with observations. In one 
instance, that of zinc, no formula was found to cover the 
whole range, and repetition gave erratic results, which may 
have been due to marked differences at different parts of the 
surface, and to the formation on collision of flattened portions 
of an unusual character. 

In calculating these constants, the usual method (as given, 
for instance, in Running, ‘ Hmpirical Formulas’) was found 
to exaggerate small errors at one end of the range. The 
following method was tlierefore employed. The slopes of 
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tangents drawn at chosen points on a smooth curve through 
the observations were recorded. Then, since 


dt\? 

(a) 
dt 


2 lo 


log (n?a?) + 2,.1—n) log », 


: 
so that (a) if log—is a linear function of log v, the formula 


dv 
= t,+a/v" follows, and (b) n and a are obtained from the 
sia and intercept of the straight-line graph indicated. 
These slopes and intercepts, it shoul Id be noted, were dclnite 
calculated by the method recently given by J. H. Awbery * 


Discussion. 

According to the formula, the duration of contact tends 
towards a steady value ty, which is greater than zero, as the 
velocity at impact increases. That to would be other than 
zero is to be anticipated, because at high speeds permanent 
deformation occurs, and during the process ef deformation 
the spheres must be in contact. It appears that at the 
highest speeds the diminution of the duration of contact 
which ought to follow from the larger elastic forces called 
into play on collision is just counterbalanced by the extension 
of this period, due to the flow of the metal. A conclusion 
reached by J. H. Vincent t, after experiments in which 
a steel ball was dropped upon a lead plate, namely, that the 
time taken by any one sphere toe produce a dent was inde- 
pendent of the velocity of impact, is another example of this 
compensating action. 

At very small speeds t=a/v” very nearly. This is remin- 
iscent of Hertz’s formula ¢v'/’=constant. In one instance, 
that of brass, n is almost exactly 1/5; while in the cases of 
lead-tin alloy, babbitt, and tin, the value of n in the first 
experiments, although greater, is not widely different from 
1/5, the difference being possibly due to working on the 
lathe. Of zine nothing definite can be said in this respect, 
while aluminium appears to be so susceptible to mechanical 
working, and so readily forms a non-conducting layer when 
annealing is attempted, that the results cannot be regarded 
as evidence on this point. It therefore seems likely that at 
very small speeds of approach Hertz’s law, or a law very 
similar to it, may apply. Further experiment is required to 


* J. H. Awbery, Proc. Phys. Soc. xli. (1929). 
+ J. H. Vincent, Proc. Camb. Phil. Soe. x. p. 382 (1898-1900). 
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establisli this point, but it may be noted that some experi- 
ments by C. V. Raman* upon slow collisions have shown 
that the coefficients of restitution for lead, brass, and bronze, 
all approach unity at small speeds of appreach. Similar 
results have been obtained by Okubo f fora number of other 
substances. 


The Effect of Mechanical Working and Annealing. 


From Table III. it appears, as a rule, that all three 
constants, to, a, and n, increase after the spheres have been 
re-turned and re-polished ; and decrease after annealing. 
n seems to be the most sensitive in this respect, and f the 
least. These facts may be understood if t) depends upon 
the properties of the metal in bulk, while n depends on the 
properties of the outside layer of the material. This is 
probably the case, for the theory of Hertz assumes that large 
stresses occur only in the immediate neighbourhood of the 
area of contact, and the elastic properties of the surface 
layers are therefore of prime importance in this theory. The 
present experiments provide a value of » similar to Hertz’s 
when the temporary deformation is confined to the surface 
and the surface has not been over worked. As soon, 
however, as comparatively large permanent deformations 
appear, metal below the surface, and possibly out of the 
influence of the lathe tool, begins to have a preponderating 
effect. This is the region in which ¢o is the chief term of 
the formula. 

Finally, these results do not support the view that, after 
contact is made, the spheres approach one another under 
a resistance proportional to the square of the velocity. 


Parr II. 


The Permanent Deformation at the Surfaces 
of the Spheres. 


When the speed of approach is fairly high, small flattened 
portions, circular in shape, are visible on the surfaces of the 
spheres. To a cursory view they appear planes, as though 
a small spherical cap had been cut off; but under the 
microscope they each appear surrounded by a pronounced 
rim, higher than the flat portion. A magnified section of 
the sphere through one of these would probably have the 
appearance of fig. 6 (a) in the majority of cases. When light 


. V. Raman, Phys. Rev. xii. p. 442 (1918). 
. Okubo, Tohéku Univ. Sci. Rep. p. 454 (1922), 


* 
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is reflected obliquely from the “flat,” it is possible to 
photograph the rim. PI. XIII. (a) shows this. It is not here 
a question of one harder sphere imposing a concave impression 
on the other, for as a rule both spheres exhibit the same 
deformation. To show this clearly, the four photographs, 
(b), (c), (d), (e) of Pl. XIII. are views of the left and right 
sides respectively of two “flats,” one upon each sphere, 
formed during a single collision. On some occasions the 
portion usually plane was itself undulating, probably having 
a section, as in fig. 6 (b). This is illustrated in (/) and (q) 
of Pls. XIII. & XIV. ° These remarks apply to all materials, 
but the case of zine requires amplification. Whereas in the 
other cases the rim was perfectly circular, or only slightly 
distorted, in the case of zinc it was nearly always irregular, 
and only roughly circular. ‘The flattened portion within 
the rim, moreover, was traversed by a network of subsidiary 


Fig. 6. 


A (c) (d) 


ridges which, although readily seen by the eye, were found 
diffienlt to record photographically. (h), (i), and (7) of 
Pl. XIV. are, however, attempts which do show how regions 
of the “flat” are cast into shadow when viewed by light re- 
flected obliquely. The irregularity of the outline is, in faet, 
exaggerated in these pictures, the outside rim, although com- 
plete, being partly hidden in shadow. The first thought was 
that the ridges were boundaries of crystals, but in this case 
the section would have been as in fig. 6(c), stepped on the 
flattened part. Illumination from the other side ought then 
to have made the bright lines dark, and shown up other 
bright lines. In point of fact the same lines were brought 
out. The conclusion is that the section must have been as 
in fig. 6 (d). 

An explanation on the following lines of these rims and 
ridges suggests itself. During the beginning of the collision, 
the metal in the neighbourhood of the point of contact is 
compressed, the pressure being greatest at the centre *. As 


* Vide Hertz, ‘ Miscellaneous Papers,’ English ed., p. 173. 
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the centres of the spheres continue to approach, the pressure 
becomes so great that, if the initial momentum be sufficient, 
breakdown occurs. The area over which the metal has 
“fused” wil] continue to expand until the balls begin to 
recede. ‘he continued approach of the centres involves 
the “fusion” of more material at the centre of the compressed 
area, and this material has to be forced aside. When the 
spheres cease to approach, and begin to recede, the pressure 
on this extruded material is relaxed, and it re-sets. ‘The 
sharply-defined rim suggests that it has been extruded beyond 
the area of contact, for here it would solidify on emergence, 
and immediately be covered by more material, thus building 
up a projecting ring which would be constrained to lie in the 
plane of contact until the pressure was released. On relaxation 
the layers first deposited would recover from the strain of 
subsequent pressure, pulling the projection into a rim. 

It may be asked why, if “fusion” occurs over a common 
area, the spheres do not adhere. This would be easier to 
explain than the apparent adhesion of the zinc spheres, for 
the surfaces are covered with exceedingly fine scratches from 
the lathe tool, each of which provides an air-film to prevent 
the coalescing of the two parts. But in the case of zinc 
transitory adhesion does appear to take place, and the ridges 
within the rim are probably the lines along which the 
adhesion was most protracted. This offers a ready explana- 
tion of the erratic nature of the results for ‘zinc, but its 
exceptional behaviour remains unexplained. 

A very similar phenomenon has been observed by 
C. V. Raman * when steel balls are dropped upon a glass 
surface. He hesitates, however, to ascribe the rim to flow 
of fused material, owing to the time available being so small. 
When, however, it is recalled that in my experiments the 
gradient parallel to the surface, of the normal pressure, is of 
the order 10! dynes per sq. cm. per cm. or more, and the 
mass of material to be moved of the order 10~* om., it does 
not seem at all impossible that the required velocity of flow 
(less than 1000 em./sec.) could be acquired in a few 
ten-thousandths of a second, even allowing for powerful 
forces due to viscous resistance. 


Size of the Permanent Deformation in Relation to the 
Velocity at Impact. 
The diameters of the circular rims were measured micro- 
scopically, and the results recorded with the velocities of 


* ©. V. Raman, J. Opt. Soc. Amer. xii. p. 887 (1926). 


796 Mr. J. P. Andrews on the 


approach in Table IV.,a typical example. These results are 
plotted in fig. 7, where d? is taken as the ordinate to bring 
LV. 
Tin Spheres. Eup. (3). 


Diameter of spheres 4:00 em. 
Mass of each sphere 240 gm. 


Diameters of rims of deformations. 


Diameters 
em. cm. cm, 
0:083 0°160 0°095 0-081 
11°93 07110 07112 0111 0°106 
16°5 0-117 0:129 0°123 0-131 
21°0 0°153 0'145 0-149 0150 
23°9 07162 0-159 0°160 
26'8 0-174 0°176 0'175 0172 
0°192 0-199 
318 0°195 0:187 071938 
38°6 0:205 0°208 0:206 0 205 
0°230 0:233 0:231 0:225 
60°2 0:235 0°235 0-248 
Fig. 7. 
d?x|0%cm2 

60} 

40} 

30 o Observations 

x Calculations 
2Qr 
tor 
(e) 10 20 30 40 50 60 70 


v cm/sec. 


out more clearly the intersection of the curve with the v axis. 
In some cases the diameters for one sphere of a pair were 
uniformly a little larger than those for the other, but in all 
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cases the mean was taken as the most convenient quantity 
to deal with. The curves for the remaining substances are 
shown in fig. 8. If the curves may be extrapolated to 
intersect the axis of v, as has been done in the figures, we 
learn that no permanent deformation is produced until the 
velocity of approach reaches a definite value vo, which is 
a characteristic of each material. There is nothing new in 
this idea, but the observation is a justification for it. The 
values of v are given in Table V. In some instances the 
diameters could not be measured with sufficient accuracy to 


Fig. 8. 


O 10 20 30 40 50 60 70 80 
Vv em Sec, 


Observations are indicated. 


warrant their inclusion here; but, from the reliable sets ot 
observations a formula was sought to represent the results. 
This was found to be 

d = b(v—w)”, 


and Table V. gives the values of the constants b, m, and v. 


Discussion. 


In the cases of aluminium, lead-tin alloy, brass, and zine 
after annealing, m is approximately 0°5. Thus d! is propor- 
tional to a quantity of kinetic energy k(v—v)?, which, when 
vis sufficiently large, differs only a little from the expression 
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for the kinetic energy of the moving sphere at collision. The 

4 
volume of the small spherical cap whose diameter is d es 
R being the radius of the sphere. Hence at high speeds we 
have the rule, originally given by Vincent for the case of 
spherical indentations in a lead plate, ‘‘ the volume of the 
material apparently removed is proportional to the kinetic 
energy of the moving sphere at impact.” 

This is, however, no longer true in general for small 
values of the kinetic energy, and according to these results 
it is not a principle applicable to all substances. A similar 
rule might be made for the other substances, but it would be 
neither so simple nor so suggestive. 


TABLE V. 
Substance. Expt. cm./seC, Me Remarks, 
Aluminium... 3 0-028 120 05 
Tin 2 0-048 25 0-41 
0:047 50 0:42 After re-turning, 
Babbitt ......43 100 0875 
1 0:045 9-0 0°37 
3 0-027 Very small. 0°50 After annealing and 
then re-turning, 
Alloy. 
Extruded 1 0:020 3-0 0:47 
Brass. 


The Least Pressure Required to initiate Flow. 


==A knowledge of v enables us to calculate this pressure 
for the various substances. The calculation is based on the 
theory of Hertz, and rests on the assumption that the 
materials may be regarded as perfectly elastic when v<v. 
From the first part of this paper, it will be clear that 
considerable justification exists for this in the cases of brass, 
lead-tin alloy, babbitt, and perhaps tin; for zinc and aluminium 
the case is not so good. 

Following Love’s version of Hertz’s work*, we find that 
the total thrast over the whole area of contact is given by 


at 


* Love, “‘ Elasticity,” 4th edn. p. 200, The English translation of 
Hertz’s ‘ Miscellaneous Papers’ contains several serious misprints, 
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where a is the radius of the compressed area, R the radius 
of the spheres, X and pw are elastic moduli commonly 
employed. 

The maximum value attained by a is related to the 
velocity of impact v thus 


max. =R(2) —2e 


where c is the velocity of compressional waves in the material 
of the spheres. 


; p= density of material ; 


a = Poisson’s Ratio ; 


while the pressure at the centre of the compressed area 


VI. 
Tensile 
Substance. Expt. P',dynes/em,?. strength, 
dynes/cm.?, 
Aluminium... {3 26 271 71x 10° 10° 
2 2:3 1°45 x 10° 1:7 
{3 2-3 1:90 10° 
Babbitt ...... 3 17 2°5) 4°3 x 10° 
3°6 2-41 8-7 10° 12x 10° 
13 3-6 aniall 
to measure. 
ead—Tin “ 
Alloy. 1 1:32 2°6x 10° 
Extruded 
— } 1 41 2:30 9°5 x 109 8 x109 


(where it is a maximum), is given when aattains its maximum 


also by 


m(A+2u)Lc?16 1—2o 
= 


where the coefficient of v?/> is denoted by ¢ 

If, now, we insert v as the velocity, P’ will bethe pressure 
just ‘capable of starting flow. Table anne gives the calculated 
values of these pressures. The necessary elastic constants 
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for aluminium, zine, tin, and brass were derived from tables, 
while those for babbitt and the lead-tin alloy were found 
from separate experiments on rods of these materials specially 
cast for the purpose. ‘These pressures may be compared with 
such quantities as crushing and tensile strengths; but the 
former are not recorded to any extent except for ferrous 
metals. The available data on tensile strengths have been 
added in the last column for comparison. It is interesting 
to note that these are of the same order as the pressures 
culculated. 


Summary of Results. 
Part I. 


(a) The duration of contact is related to the velocity of 


approach by the law ¢ = ty + a 

(6) For the alloys brass, lead-tin, and babbitt n is: 

practically 1/5, as in Hertz’s formula tv!” = constant, provided 

the surface is not mechanically worked to too great an extent. 
8 


(c) For the metals tin and aluminium n is definitely 
greater than 1/5, while for zine no value can be given. 


(d) Mechanical working of the surface raises the values 
of to, a, and n, while annealing reduces them again. 


Part 


(a) The circular permanent deformations are found to be 
bounded by a rim. 


(>) The diameters of these rims are related to the velocity 
at impact by the law d = b(v—w)”. 


(c) In the cases of aluminium, lead-tin, and brass m is 
practically 0°5. 

(d) In the case of zinc the deformations are irregular in 
outline and intersected by ridges. 


(e) The pressures necessary to cause the metal to flow are 
of the order 10° or 10! dynes/em.?, that is, of the order of 
the tensile strengths. 


In conclusion, this is the first paper upon a research, still 
in progress, carried out at East London College under the 
encouragement of Prof. C. H. Lees, to whom I am indebted 
for many helpful suggestions. 
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LXXXVII. On the 42-Minute Period in the Frequency of 
the Ajter-Shocks of Earthquakes. By Cuartus Davison, 
SeD., F.G.S.* 


fone time taken by an earthquake-wave to travel from a 

focus near the surface to its antipodes is almost 
exactly 21 minutes. As the crust within and near the focus 
is for some days in a highly sensitive condition, it is possible 
that the return-pulsation may affect the frequency of the 
after-shocks, and my object in this paper is to show that a 
42-minute periodicity does govern their occurrence. In two 
earthquakes, the Mino-Owari earthquake of 1891 and the 
Kwanto earthquake of 1923, I have applied the same 
methods in the hope of tracing periods of 21 and 84 minutes, 
but I have not been able to find any decisive evidence of 
their existence 

Let us suppose that the crust is on the point of motion in 
one direction over a certain area, and that a periodic foree 
of strength sufficient to precipitate motion acts upon it in the 
same direction. The resulting earthquakes would then be 
subject to a period of the same duration, and the maximum 
epoch of the period would coincide approximately with that 
of the variable force. If the latter force were to act in the 
direction opposite to that in which the crust is about to move, 
there might still be the same periodicity in the frequency of 
earthquakes, but the maximum epoch would coincide with 
the minimum epoch of the variableforce. ‘Thus, if the after- 
shocks of an earthquake are due mainly to a continued exer- 
tion of the forces that give rise to the earthquake, then the 
maximum epoch of the 42-minute period, if it exists, should 
coincide with the arrival of the return-movement. On the 
other hand, if most of the after-shocks are due to the settling 
of the displaced crust, the maximum epoch of the period 
should oceur about 21 minutes after that return. 

The only catalogues of after-shocks at my disposal that 
are full enough for the purpose are those of five Japanese 
earthquakes from 1891 to 1927, three European earthquakes 
from 1894 to 1915, and three others depending on personal 
observations only. 

In the following sections the time of occurrence given 
is that of the initial epoch of the movement. As, however, 
the duration of a great earthquake may amount to 4 or 5 


* Communicated by the Author. 
+ On this subject, see H. H. Turner, R. Astro. Soc. Monthly Notices, 
Geoph. Suppl. vol. i. pp. 31-50, 88-99, and vol. ii. pp. 73-76, 
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minutes, the epoch of the strongest movement may be a 
minute or two later. In each earthquake, successive inter- 
vals of 42 minutes were divided into 14 minor intervals 
of 3 minutes each, and 7-interval means were taken of the 
total numbers of after-shocks occurring in corresponding 
intervals. In the tables the date of the maximum epoch is 
the number of minutes after the last complete interval of 
42 minutes measured from the initial epoch of the earth- 
quakes. In estimating the amplitude the average number 
of after-shocks in each 3-minute interval is taken as unity. 


Mino-Owari Earthquake of Oct. 28, 1891. 


The Mino-Owari earthquake occurred at 6.37 a.m. The 
seismographs at Gifu and Nagoya were at once put out of 
action, but the first after-shock was recorded at 1.10 P.m. at 
Nagoya and 1.55 p.m. at Gifu*. In both records the 
intervals mentioned range from Oct. 28, 2.0 p.m., to Nov. 2, 
1.0 p.m., from that date to Noy. 6, 0.12 p.m., and from the 
latter to Nov. 10, 11.24 P.M. 


Girv. Nacoya. 
Me M 
No. of ee No. of Hex 
Interval. epoch, Ampl. epoch, Ampl, 
min, min. 
Oct. 28-Nov.2... 863 183 09 507 33 
Nov. 2-6 201 92 "22 91 93 ‘47 
Noy. 6-10 ......... 177 33 65 214 “41 


Thus, during the first five days, the minima of the 
42-minute period at Gifu, and the maxima at Nagoya, 
coincided approximately with the return-movements. In the 
next four days, the series of means were somewhat irregular, 
and the epochs at both places occupied an intermediate 
position. In the last four days the epochs were almost 
reversed, so that the maxima at Gifu and the minima at 
Nagoya coincided nearly with their return-movements. 

The explanation of the opposition in epoch at Gifu and 
Nagoya is probably to be found in the positions of these 
places with reference to the great fault-line. Along this 
fault the crust on the east side was relatively lowered by 
10 feet or less, except at one place, Midori, where it was 


* F. Omori, Tokyo Imp. Univ. Coll. Sci. Journ. vii, pp. 178-189 
(1894). 
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raised by 20 feet. Gifu lies about 5 miles from the fault- 
line and 14 miles from Midori. Nagoya lies about 18 miles 
from the continuation of the fault-line and 34 miles from 


Midori. 


Hokkaido Earthquake of Mar. 22, 1894. 

The Hokkaido earthquake occurred at 7.56 p.m., and the 
after-shocks were recorded at Nemuro (75 miles from the 
epicentre) without interruption at the start*. The first five 
days range from Mar. 23, 0.0 4.m., to Mar. 27, 11 p.m, and 
the next ten days from the latter time to Apr. 6, 9 P.M. 


Interval. No. of shocks. Max. epi, Amplitude. 
min. 
Mar. 23-27 ...... 264 25 ‘17 
Mar. 27—Apr. 6 ... 55 233 “25 


Thus, in the first interval, the maxima of the 42-minute 
period coincided approximately with the return-movements ; 
in the second interval the minimum coincided nearly with 
those returns. 


San-riku Earthquake of June 15, 1896. 


The epicentre of this earthquake lay near the foot of the 
western slope of the Tuscaroora Deep, about 160 miles from 
the east coast of Japan. It was registered at Tokyo (about 
340 miles from the origin) at 7.34 p.m.t During the five 
days, June 15-19 (36 atter-shocks), the maxima of the 
42-minute period occurred at 224 min. after the return- 
movements (amplitude -34). After June 19 the recorded 
shocks are too few in number to determine the existence of 
the period. 


Kwanto Earthquake of Sept. 1, 1923. 


The great Kwanto earthquake was recorded at Tokyo at 
11h. 58m. 44s.,4.m. Most of the seismographs in the 
central district were overthrown by the shock, but one 
instrument at Tokyo remained in action throughout. ‘l'wo 
lists of after-shocks have been published: one, by A. lmamura 
and K. Hasegawa, of the after-shocks strong enough to be felt. 
in Tokyo from Sept. 1 to Oct. 1; the other, by T. Yasuda, of all 


* I, Omori, Imp, Harthq. Inv. Com, Publ. no. 7, pp. 85-51 (1902). 
+ Brit. Ass. Rep. 1897, pp. 7-8. 
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the after-shocks recorded in that city from Sept. 1 to Jan. 31, 
1924*. The first after-shocks entered in both lists occurred 
at Oh. 1m. 49s., p.m. on Sept. 1. For the remainder of 
Sept. 1, Imamura gives the times of 189 shocks and Yasuda 
of 225. For the whole of September, the corresponding 
numbers are 720 and 1303. 

Making use first of Inamura’s list, the maximum epoch for 
Sept. 1-6 (463 shocks) occurs at 20} min. with an ampli- 
tude of 13. During the second and third five days (114 
and 45 shocks) the existence of the period becomes uncertain, 
the series of means being irregular. For Sept. 16-30 (84 
shocks) the series of means vary regularly, with the maximum 
epoch at 324 min. and the amplitude °30. 

The central district is divided by Imamura and Hasegawa 
into four regions, to each of which belonged its own series of 
after-shocks. The region A includes most of Sagami Bay 
and the land to the north; the region B contains the Boso 
Peninsula and the sea-bed to the east; the region C lies 
entirely on land to the north of Tokyo; while the region D 
includes Tokyo and the surrounding country. Of these 
regions, A and B are by far the most important ; the region 
C was disturbed by comparatively few after-shocks, and the 
region D by a larger number, usually of slight intensity. 

Nearly 24 hours after the great earthquake, at 11h. 46m. 
55s., A.M., on Sept. 2, another violent shock, almost as 
strong as the first, occurred in region B with its epicentre 
lying about 41 miles to the east. It is probable that most 
of the after-shocks in this region were connected with this 
second earthquake, for only one shock on Sept. 1 is assigned 
to region B. 


Max. epoch, 


Region. Interval. No. of shocks. aia Amplitude. , 
Sept. 6. 239 16 
Sept. l- 6. 89 203 ‘47 
Bie Sept. 7. 96 36 "33 


For regions A and D the first five days are reckoned 
from noon on Sept. 1, the second five from 11 a.m. on 
* A, Imamura and K, Hasegawa, Imp. Harthq. Inv. Com. Bull. ii. 


pp. 72-93 (1928); T. Yasuda, Imp. Earthg. Inv. Com, Rep. no. 1004, 
pp. 261-310 (1925). 
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Sept. 6, and the maximum epoch from 11.59 a.m. on Sept. 1. 
For region B the first five days are reckoned from noon on 
Sept. 2, the second five from 11 a.m. on Sept. 7, and the 
maximum epoch from 11.47 a.m. on Sept. 2. 

Thus, in region D, the minima of the period during the 
first five days coincided closely, and in region A roughly, 
with the return-movement ; during the second five days the 
maxima in region A coincided closely with those returns. 
In region B, during the first and second 5-day intervals, 
the maxima and minima, respectively, coincided roughly 
with those returns. 

Yasuda’s list differs from the other in the inclusion of 
after-shocks imperceptible to human beings in Tokyo. It is 
also eontinued for four additional months. The results are 
as follows :— 


Max. epoch 


Interval. No. of shocks. Ath. Amplitude. 
755 234 14 
133 40 “16 
45 41} 38 
» 15-31 . 51 24 “22 


During the intervals Sept. 6-11 and Sept. 11-16 there 
is no irregularity in the series of means, and the date of the 
epoch is probably that of the resultant of two periods with 
maxima near the middle and end of the 42-minute intervals. 
During November there is no trace of a 42-minute period 
in either of the intervals Nov. 1-14 and Nov. 15-30. Thus, 
after two months the throbbing of the earth due to the 
original displacement seems to have died away. 

Lastly, deducting the number of shocks for each interval 
in Imamura’s list from the corresponding number in 
Yasuda’s list, we obtain the numbers of shocks imperceptible 
to human beings in Tokyo, with the following results :— 


Max. epoch, 


Interval. No. of shocks. : Amplitude, 
min. 
Sept, 292 235 28 
174 7 25 
79 4 30 
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Tango Earthquake of Mar. 7, 1927. 


The earthquake occurred at 6.27 p.m. It originated along 
two faults nearly at right angles to one another. Surveys 
made in May-June, August-September, and October— 
November of the same year showed that the movements of 
the crust continued for at least six months of the earthquake, 
and not always in the same directions. The after-shocks 
were recorded at three stations, Maiduru, Kinosaki, and 
Inemura, about 8, 12, and 13 miles, respectively, from the 
line of the principal. fault*. The first after-shock was 
reeorded at Maiduru at 10.12 p.m. on Mar. 11, and at the 
other stations on the following day. The first interval in 
the table (see p. 807) ranges from Mar. 12, 3.12 P.m., to 
Mar. 16, 11 P.M. 

A strong after-shock occurred at 6.7 A.M. on Apr. 1, and 
this was followed by 191 after-shocks on the same day, so 
slight that only 12 were recorded at all three stations, 20 at 
one other station besides Maiduru, and 162 at Maiduru only. 
The list for the latter station gives the following results :— 


Max. epoch, 


Interval. No. of shocks. ses: Amplitude. 
191 13 23 
58 
53 224 “44 
» 8-14 104 


The results given in the above tables seem to depend on 
the complex origin of the principal earthquake. From 
Mar. 12 to 31 the maximum epoch isat about 30 min., in April 
at about 284 min. Then after about seven weeks the epoch 
is reversed (to 74 min.) in May. In June the number of 
shocks recorded is only 17. The list for July seems to 
contain misprints, and ] am unable to make use of it. In 
August there is no trace of the period. The epochs for Apr. 1 
and Apr. 4-7 are reckoned, as before, from the principal 
earthquake at 6°27 p.m. on Mar.7. As the interval between 
the principal earthquakes and the strong after-shock at 6.74.M. 
on Apr. J exceeds a multiple of 42 minutes by 20 minutes, it is 
at first sight doubtful whether the maxima coincide with the 
return-movements due to the principal earthquake, or the 
minima with the retutn-movements from the atter-shock.. 


* N. Nasu, Tokyo Earthq. Res. Inst. Bull. vi. pp. 300-381 (1929). 
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That the former alternative is the correct one seems clear 
from the frequency of the after-shocks during the first three 
42-minute periods succeeding the strong after-shock. Taking 
the number of after-shocks in 3-minute intervals and 
smoothing the curve by taking for each interval the mean of 
three successive intervals, the maxima of the curve occur 
close to the beginning of each of the 42-minute periods. 
Thus, as the first maximum occurs within a few minutes 
after the strong after-shock, the period must be connected 
with the return-movements from the principal earthquake 
rather than with those (if they existed in sufficient strength) 
from the after-shock. It follows, therefore, that the move- 
ment that gave rise to the after-shock was in the opposite 
direction to the movement that caused the great earthquake. 


Locris Greece) Earthquakes of 
Apr. 20 and 27, 1894. 


The first earthquake occurred on Apr. 20 at 6 52 P.m., the 
second and stronger on Apr. 27 at 9.30 p.m. The after- 
shocks were recorded at Athens (about 60 miles from the 
epicentre)*. The first interval ranges from April 20, 7.0 P.m., 
to Apr. 27, 5.36 p.m., the second from Apr. 27, 9.48 P.M., to 
May 12, 6.48 p.m. 


Max. epoch, 


Interval. No. of shocks, : Amplitude, 
min. 
Apr. 21 18 52 
2i—May 12 ... 31 23 


Thus, for the first interval, the maximum epoch agrees 
nearly with the return-movements. 


Messina Earthquake of Dec. 28, 1908. 


The earthquake occurred at 5.20 a.m. The first after-shock 
was recorded at Mileto at 5.36 a.m. and at Catania at 
6.4 a.m.{. In the following table the first interval ranges 
from Dec. 28, 6.0 am., to Jan. 2, 5.0 A.m., and the second 
from the latter time to Jan. 12, 3.0 a.m. 


* D. Eginitis, Athénes Obs. Nat. Ann. ii. pp. 217-229 (1900). 
+ Ital. Sism. Soc. Boll. xv. pt. 2, pp. 566-689 (1911); xvi. pt. 2, 
pp. 1-50 (1912). 
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Max epoch, 


Record. Interval. No. of shocks. ae! Amplitude. 
General ... Dec. 28-Jan. 2. 107 144 28 
Jan 2-12. 38 352 25 
Mileto...... Dec. 28—Jan. 2. 58 114 33 
Jan, 2-12. 34 824 45 
Catania... Dec, 28-Jan. 2. 94 143 36 


Taking the after-shocks as a whole, the minima of the 
42-minute period preceded the return-movements by 64 min. 
during the first five days. In the second interval the epoch 
is reversed. At Mileto the epoch during the first five days 
preceded the return-movements by 94 min.; ee the next 
ten days it was reversed. During the months of Febr uary 
and March there is no distinct trace of the period. 


Marsica Earthquake of Jan. 13, 1915. 


The earthquake occurred at 7.53 a.m. The record of 
after-shocks at Rocea di Papa (47 miles from the epicentre) 
is a very full one, the first after-shock being registered at 
8.27 am.*. The earthquake was unique among those here: 
considered. Though very destructive the area of damage 
contained only 55 sq. miles and the disturbed area not more 
than 5700 sq. miles. ‘Uhe first interval of five days (554 
after-shocks) gives the maximum epoch at 10 min. and 
the amplitude as ‘06. In succeeding intervals of five, five, 
and ten days the series of means are somewhat irregular ; 
there is no sign of a reversal of epoch and no certain trace 
of a 42-minute period after Jan. 18. Making use of shorter 
intervals, it is clear that the reversal occurred very soon 
after the earthquake, as follows :— 


Max. epoch, 


Interval, No. of shocks. Ampl. 
First six hours (6 hr. 18 m.). 177 19 19 
Third sala 42 114 67 
Fourth 32 24 +30 
Rest of first five days......... 251 * 40 09 


* A. Cavasino, Ital. Sism. Soc. Bol. xix. pp. 286-291 (1915), 
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Thus, during the first six hours, the minimum epoch 
coincided nearly with the return-movement; during the 
next twelve hours the period seems to be the resultant of 
two periods with different epochs ; in the fourth six hours 
and during the rest of the first five days the maximum epoch 
coincided nearly with the return-movement. The early 
reversal of epoch and the brief duration of the 42-minute 
pulsations are no doubt connected with the small area affected 
by the earthquake. 


Hawaiian Earthquakes of Mar. 28 and Apr. 2, 1868, 


In these and the remaining earthquakes, the records of 
after-shocks are non-instrumental. The period considered is, 
however, too short to be affected materially by the varying 
conditions of the observers throughout the day. 

The first strong earthquake occurred on Mar. 28 at 
1.28 p.m., the second and more destructive on Apr. 2 at 
3.40 p.m.* The first interval ranges from Mar. 28, 2.6 P.M., 
to Apr. 1, 6.18 a.m., the second from Apr. 2, 4.12 P.m., to 
Apr. 7, 4.36 P.M. 


Max. epoch, 


Tuterval. No. of shocks. Amplitude. 
Mar. 28-Apr. 1... 118 244 23 


Thus, the minima in the first interval and the maxima in 
the second interval occurred shortly after the return- 
movements. 


Riviera Earthquake of Feb. 23, 1887. 


The principal earthquake occurred at 6.20 a.m. The first 
interval ranges from Feb, 23, 9.6 A.M., to Feb. 28, 8.6 a M., 
and the second from the latter time to Mar. 10, 6.6 a.m. f. 


Interval. No. of shocks. Amplitude. 
Feb. 23-28 ......+4 79 203 
Feb, 28-Mar: 10 ... 34 "45 


* H.0. Wood, Amer. Seis. Soc. Bull. iv. pp. 178-183 (1914), 
+ T. Taramelli and G. Mercalli, Ug. Centr. di Meteor. e Geodin., 
Annali, viii. pp. 260-289 (1888). 
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Thus, in the earlier interval, the maximum epoch occurred 
at nearly 21 minutes after the return-movements. In the 
later interval the number of shocks is small, but the epoch 
tends to agree with those movements. After Mar. 10 there 
is no trace whatever of the 42-minute period. 


Zante Earthquake of Apr. 17, 1893. 


The earthquake occurred at 7.4a.m. The first interval 
ranges from Apr. 17, 7.0 a.m., to Apr. 22, 6.0 a.m., the 
second from the latter time to May 2, 4.0 a.m. * 


Max. epoch, 


Interval. No. of shocks. Amplitude. 
Apr. 17-22 ...... 68 214 “45 
Apr. 22-May 2... 46 303 “34 


Thus, for the first five days, the minima of the 42-minute 
period coincided closely with the return-movements, 


Conclusions. 


(i.) In each of the series of after-shocks examined, there 
is evidence of the existence of a 42-minute period, connected 
more or less closely in point of time with the epoch of the 
earthquake. As a rule, the maxima at first occur about 
21 minutes after the returns of the earthquake-movement. 
This periodicity lasts ror about five days—in one earthquake 
for a few hours and in two others for three or more weeks— 
but is continued for several weeks longer with a reversal of 
epoch. 

(ii.) The after-shocks of a great earthquake seem to have 
two distinct origins. Some, including perhaps the stronger, 
after-shocks result from repetitions of the original movement, 
transferred in part to other regions and especially to those 
bordering the area of displacement. Others are due to the 
seitlement of the displaced mass. At first the latter greatly 
predominate in number, but after about five days the former 
begin to assert themselves. The reversal of epoch in the 
42-minute period is due, not to any change in the variable 
force, but to the change in the predominating cause of the 
after-shocks. The change, however, does not take place 
abruptly, and there are intervals, lasting it may be for some 


* Athenes Obs. Nat. Ann. ii. pp. 195-196 (1900). 
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days, in which neither kind of after-shock distinctly pre- 
dominates. At such times either the maximum epoch 
occupies an intermediate position, or the series of means is 
irregular and leads to no definite period *. 

(iii.) Three of the earthquakes here considered were 
followed, after an interval of from one to seven days, by an 
earthquake of nearly the same, or even greater, intensity. 
In the Kwanto earthquakes of 1923 and the Hawaiian 
earthquakes of 1848 the return-movements coincided nearly 
with the minima of the 42-minute period in the earlier 
earthquake and with the maxima of the period in the later, In 
the Locris earthquake of 1894 the change of epoch was less 
pronounced, but the number of shocks recorded were small 
(21and 31). The evidence, so far as it goes, seems to suggest 
that, when two principal earthquakes occur in one series, 
the movement causing the earlier is in the opposite direction 
to that which gives rise to the first after-shocks, while the 
movement causing the later is in the same direction. 

(iv.) So far as I am aware, the 42-minute period is con- 
fined to the after-shocks of great earthquakes. I have 
examined many records of ordinary earthquakes, both for 
whole years and for single days, but have failed to find any 
certain trace of the period, except perhaps for a few of the 
years of the Philippine earthquakes. 

(v.) Next to the existence of the 42-minute period, the 
most interesting point is the length of time for which it 
lasts. After the great displacement that gives rise to an 
earthquake, the earth continues to throb for at least ten days, 
and, in the Kwanto and Tango earthquakes, for not less than 
two and three months, respectively. Now, during the latter 
half of the nineteenth century, Milne chronicles 219 destruc- 
tive earthquakes of the greatest intensity and 368 of the 
next degree, or, on an average, 12 a yearyt. As there are 
many other earthquakes in uninhabited regions of equal or 
greater strength, it would thus seem possible that the earth 
is seldom, if ever, entirely at rest from the quivering that 
follows great earthquake-producing displacements. 


* A similar reversal of epoch occurs in the diurnal periodicity of 
after-shocks (Phil. Mag. xli. pp. 914-915 (1921). Mr, K. Suda has 
also noticed that, after three or four days, the frequency of the after- 
shocks of the Kwanto earthquake of 1923 at various stations. was some- 
what abruptly reduced (Kobe [Japan] Kup. Mac. Obs. Mem. i. pp. 181— 
186, 1924), 

+ Brit. Ass. Rep. 1911, pp, 709-740 
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LXXXVIII. The Explosion of Hydrogen-Air Mixtures in a 
Closed Vessel.—Part I. By B. H. Toorp, 


Introduction. 


URING the last twenty years several research workers 
have been measuring the maximum pressures developed 
in the explosion of gaseous mixtures in closed vessels. The 
purpose of this paper is to survey the results obtained for 
mixtures of hydrogen, oxygen, and nitrogen, and to attempt 
to account for the extraordinary differences between the 
figures given by various workers. A new diaphragm 
indicator will also be described, the development of which 
has emphasized the importance of certain points that have 
generally been neglected. 

As an example of the disagreement which at present exists, 
the following results may be quoted for a 29-6 per cent. 
H,+14°8 per cent. O, +55°6 per cent. N, mixture exploded 
at an initial pressure of one atmosphere :— 


I. 
Ratio of Temperature. Apparent mean 
According maximum to ~ specific heat 
to initial Initial Maximum of steam 
pressure. tm O. 
R. W. Fenning ......... 778 51:0 2397 
Maxwell and Wheeler § .. 15:0 2420, 


The pressure ratios are all corrected to an initial tempera- 
ture of 15° C. to make them comparable (see Appendix C). 
‘Lhese figures have been obtained by interpolation ; the dif- 
ferent workers have seldom used mixtures of identical 
composition, and direct comparison of pressure ratios is 
difficult. The most convenient method seems to be that. 
of comparing the “apparent mean specific heat” of steam 


* Communicated by Prof. W. T. David, Sc.D. 

+ Zeit. f. Elektrochem, xv. p. 536 (1909). The original mixture con- 
tained 2°2 per cent. H,O and a small amount of excess H, (included in 
the N,). 

t :. Res. Comm. Report, No. 902, April 1924. This pressure ratio 
is 2 per cent. higher than that given by Fenning; in a later paper he 
says that certain earlier results were 2 per cent. too low, but it is not 
clear just which. 

§ J.C.S. p. 15 (1928). Original mixture contained 1°7 per cent. 
H,0 (included in the 


Phil. Mag. 8.7. Vol. 8. No. 53. Suppl. Dec.1925. 31 


814 Mr. B. H. Thorp on the Explosion of 


between 15° C. and the explosion temperature, calculated 
from each explosion. The method of obtaining these and 
the figures used for the specific heats of the diluent gases 
are shown in Appendix A. Allowance has been made in 
all eases for heat-loss on the basis of actual measurements 
by Prof. W. 'T. David (Appendix B). In calculating the 
apparent mean specific heats, it has been assumed that there 
was no dissociation of steam or incomplete combustion at 
maximum pressure. Hence the extent to which they exceed 
the true mean specific heat of steam is a measure of the 
amount of energy existing in the gases at maximum pressure 
in forms other than that of temperature. 

In view of this disagreement, it was decided some three 
years ago to undertake a comprehensive set of explosions 
under the direction of Professor David in this laboratory, 
under the best possible conditions and covering a very wide 
range of H,—O,—N, mixtures. The full results of these 
will follow in the second part of this paper ; for the present 
they will only be considered in relation to the results of 
other workers. With one exception, all the explosions con- 
sidered in this paper have been made at an initial pressure 
of one atmosphere. 


Description of Apparatus employed. 


The explosion vessel is a large cast-steel sphere, 45-4 em. 
diameter, with the inner surface silver-plated in order to 
reduce heat-loss before maximum pressure. It is made in 
two halves flanged and bolted together ; before each explosion 
it is opened up and the plating polished with rouge. Hxtreme 
care is taken in the drying, mixing, and analysing of the 
explosive mixtures; the calculated and observed chemical 
contractions usually agree to 1 part in 1000 (except where 
NO, is formed). 

At first a Hopkinson indicator was used, but this was not 
found suitable for rapid explosions. In these it appears to 
overshoot, owing either to the inertia of its moving parts 
or to the pocket below the piston. A new diaphragm 
indicator (fig. 1) was therefore evolved with a view to 
overcoming these defects. 

A flat diaphragm, 2°8 cm. diameter and 0°75 mm. thick 
with a solid rim, is turned out of Invar. It is clamped down 
very firmly on to a seating turned in the wall of the explosion 
vessel, with only ashallow pocket below it; the joint is made 
with a thin ring of soft copper. The glass mirror is stuck 
on to a small brass pedestal, which is itself soldered on to 
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the diaphragm at about half the radius; the mirror and 
pedestal together weigh only 1 gm. A very long beam of 
light (0°65 m.) is used to obtain a high magnification, and 
some months were spent in overcoming the optical difficulties 
involved. ‘The fine, sharp line now obtained enables the 
deflexions on the films to be measured with an accuracy of 
1 in 500. 

At first steel diaphragms were used, but it was found that 
they were giving misleading results, and a large amount of 
work had to be thrown away. Although protected from the 
flame by mica they were heated and consequently distorted 


Fig. 1. 


Q{ 

: 


New diaphragm indicator. (Half full size.) 


early in the explosion. The distortion due to heating was 
opposed to that due to increase of pressure inside the vessel, 
and so the maximum pressure recorded wis some 1 percent. 
too low. The use of invar has reduced this error to a negli- 
gible quantity. ‘Tests on an invar diaphragm in working 
position showed that a bunsen flame could be held beneath 
it for at least three seconds before any measurable deflexion 
took place ; whereas with a steel diaphragm a_ visible 
deflexion occurred as soon as the flame touched it, 

The use of invar had, however, one quite unforeseen 
effect. While calibrating the indicator it was found that 
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the temperature of the sphere had a marked influence on the 
deflexion caused by a given pressure ; a temperature rise of 
3° reduced the deflexion 1 per cent. This was at first put 
down to increase of the modulus of elasticity of invar with 
temperature, though the data available did not suggest an 
increase of such magnitude. It was of vital importance 
to know the cause, as in an average explosion the mean 
temperature of the diaphragm probably rises some 2° before 
maximum pressure is reached. If the modulus increases 
appreciably at the same time, the scale to which the deflexion 
shall be reckoned is doubtful. The following experiment 
seemed to show, however, that the modulus does not vary 
enough to upset the results. During calibration the deflexion 
produced by a given pressure was determined under the 
following conditions :— 


(a) When the diaphragm was at the same temperature 
as the sphere (15° C.). 


(6) When some 3 c.c. of half melted snow were dropped 
on to the diaphragm a couple of seconds before 
the pressure was applied. 


(c) When the same was done with water at about 
60° C. 


No appreciable difference in the deflexions could be 
measured, and since it can hardly be doubted that the 
temperature of the diaphragm varied several degrees, it 
seems clear that it is the temperature of the sphere as a 
whole whose variations affect the pressure-deflexion scale. 
It is therefore only necessary to calibrate the indicator over 
a sufficient range of initial temperatures. The probable 
explanation is that as the invar diaphragm does not expand 
thermally along with the sphere as their common tempera- 
ture rises, but yet is rigidly clamped, it is mechanically 
stretched and so yields less to the pressure inside the vessel. 

Particular care is taken over the calibration of this vessel. 
The sparking plug is removed from the bottom of the sphere 
and a tube, passing up through the hole, is serewed into the 
pocket beneath the diaphragm, the joint being made with a 
rubber ring. A three-way cock at the bottom of this tube 
enables the pressure below the diaphragm to be raised sud- 
denly from atmospheric to that of a large tank (170 litres) 
close by full of compressed air. A record is taken on a film 
showing the actual jump up of the pressure, the rapidity 
of which is of the same order as that taking place in an 
explosion, and the immediate deflexion measured. The 
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pressure in the tank is measured with a compound mercury 
gauge to 0-1 per cent., and since it is easy to determine the 
corrections for such a gauge, the absolute values of the cali- 
bration pressures are above suspicion. It is of the utmost 
importance that the diaphragm should be calibrated in its 
working position, since the firmness of clamping has a 
marked effect on the pressure-deflexion scale. 


Comparison of Experimental Results with those 
of other Workers. 


Some of the experiments of the following workers have 
been selected as suitable for comparison with those made 
here recently :— 


M. Pier* (diaphragm indicator; 40°6 cm. diam. sphere). 

R. W. Fenning + (diaphragm indicator ; cylinder 17°8 
em. diam., 20°3 cm. long). 

Maxwell and Wheeler} (Petavel indicator; 19:7 cm. 
diam. sphere). 

W. Siegel § (diaphragm indicator ; 40°6 cm. diam. 
sphere). 

W. A, Bone || (Petavel indicator ; 7.6 cm. diam. sphere ; 
3 atm. initial pressure). 


Particulars of these, and the apparent specific heats of 
steam derived from them, will be found in Table III. 

Fig. 2 shows the apparent specific heats plotted against 
the maximum explosion temperatures for four series of 
hydrogen-air mixtures—those of Maxwell and Wheeler, 
Fenning, Bone, and one from this laboratory (Leeds). 
It is evident that there are considerable differences, and 
there are good grounds for suggesting that these are chiefly 
due to the size and shape of vessel used. The Leeds experi- 
ments, made in the largest vessel of the four, a sphere 
45-4 cm. diameter, give the lowest apparent specific heats— 
i. e. the highest pressure ratios. Maxwell and Wheeler with 
a sphere 19-7 cm. diameter get intermediate results, while 
Bone with a very small sphere, only 7°6 cm. diameter, gets 
the highest apparent specific heat. Henning, whose results 
are almost as high as Bone’s, used a vessel of about the same 
size as Maxwell and Wheeler’ s, but cylindrical in shape. 
Referring to Table I., it will be seen that Pier with a 


* + Loe. cit. { Loe. 
§ Zeit. f. Phys. Chem. \xxxvil. p. 641 (1914), 
|| Proc. Roy, Soc. A, cviii, (1925). 
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40°6 cm. sphere got a much lower apparent specific heat 
than Maxwell and Wheeler. 

Experimental confirmation of this effect of size and shape 
of vessel on the maximum pressure developed has been 
obtained in this laboratory. Hach mixture was made up 
in a large tank, drawn off from there to four different 
vessels (all with rough, black, inner surfaces) and exploded 
in them under identical conditions. The indicators were 
calibrated simultaneously, and every precaution taken to 
ensure strict comparability. Two or more explosions were 


Fig. 2. 
iFennin 
QO | | 
Wheeler. 
| 
17-0 
eeds 
Siro, 
& 100 
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Temperature t°C. 


Hydrogen-air mixtures. 


made in each vessel for each mixture, and several mixtures 
were tried. The following results are typical :— 


Taste II. 
Maximum pressure developed in 
Mixture. 46-4 em. 303 em. 15lcem. 15:1 em. 
: sphere. cylinder. sphere.  cyli: der. 
23 % 00O+8 % H,+74 % O,... 773 7°60 7-53 7:43 


20 % 0046 % H,+74% 0,1... 7°69 7-46 7°58 7°36 
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If it be assumed that the absolute rate of heat-loss for 
a spherical vessel varies as the area of the inner surface 
(2. e. as the diam.”), then the rate of heat-loss, expressed as 
a percentage of the heat of combustion (which « diam.*), 
varies inversely as the diameter. On the other hand, the 
time taken to attain maximum pressure varies roughly as 
the diameter, so that the total percentage heat-loss up to 
the moment of maximum pressure is almost independent of 
the diameter. In a cylinder it will be perhaps 50 per cent. 
greater than in a sphere. But it is believed (on the strength 
of Professor David’s measurements, see Appendix B) that 
in these rapid explosions it amounts to less than 1 per cent. 


Heat 
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and that the differences in the pressure ratios must be due to 
some other far more important factor. 

The low pressure ratios recorded by some workers may be 
partly due to their use of “static” calibration. Pier and 
Siegel calibrated their indicators by sudden application of 
compressed air (as is done here), approximating to explosion 
conditions. The others appear to allow the calibrating pres- 
sure to act on the indicator fer some time, so that elastic 
creep may occur before the deflexion is measured. Pier 
measured this creep and found that it might cause an error 
of the order of 1 per cent. 

Fig. 3 shows apparent specific heat curves derived from 


| | 
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the series of excess hydrogen-oxygen explosions by Pier and 
Siegel and a third series from this laboratory, all in large 
spheres of much the same size. It appears that the Leeds 
experiments are in fair agreement with the others, though 
they suggest that the maximum pressures recorded by Pier 
and Siegel are from 1 to 2 per cent. too high. When cali- 
brating they moved their diaphragms to a special vessel (of 
very small volume, so that it might be filled rapidly to the 
required pressure with compressed air), and though these 
may not have been so susceptible to alterations in clamping 
as the Leeds diaphragms, this throws some doubt on the 
absolute values of their recorded pressures. They also used 
steel diaphragms, and our experience has shown (see above) 
that the distortion of these by heating early in the explosion 
may cause an error of 1 per cent. either way (depending on the 
particular diaphragm) in the maximum pressure as calculated 
from the deflexion. Every care has been taken in this 
laboratory to obtain very great accuracy in experimenting, 
and it is believed that the pressures are within -+-4 per cent. 
of the true values. 


Discussion of Results. 


It has been shown how marked an effect the size and 
shape of the explosion vessel have upon the maximum 
pressure developed in it for any given mixture of. gases, 
and that this cannot well be explained by heat-loss before 
maximum pressure. The most likely explanation seems to 
be that it is due to a time lag in the process of combustion. 
If the overall reaction H,+40,—+H,O + 57,290 eal., after 
being started in any small volume of gas by the flame 
travelling through it, takes an appreciable time to complete 
itself, then the larger the vessel and the more closely it 
approximates to a sphere with central ignition the more 
advanced will the reaction be on the average throughout 
the vessel at the moment of maximum pressure (or at any 
moment after the flame has reached the walls), and the 
greater the proportion of the total energy showing itself as 
temperature and pressure. 

There is abundant evidence from other sources to support 
this theory of ‘incomplete combustion at maximum pressure,” 
and it will be discussed more fully in the second part of this 


paper. 
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III. 
Max. temp, 15%, 
13:3 2°0 1417 81 
15:0 20 1592 79 
tas 20°0 18 2017 88 
Siegel. «x°/, (H,.+30,)+y°/, He. 
Max. temp. 15°, t 

Be. Th. O. 

22°2 2205 9°3 

26°8 2535 9°9 

27'8 2593 10°1 


Fenning. #°/, H2+y°/, Oot3°9 y°/, No. 


252 2188 126 
20 2306 13-4 
893 2330 12:8 
346 2314 125 


Bone, 29°6 °/, H,+14'8 °/, O.+55°6°/, Np at 3 atm. 
initial pressure, gives 


Max. temp.=2347°O. and O = 138. 
2 


Maxwell and Wheeler. 2°/, Ho+y°/, Og+3°8y°/, No. 


Max. temp. 15%, ° 
tm® C. AH,O. 
148 1385 
1971 1688 12:9 
24'1 2090 122 
28°4 2312 
33°5 2806 
36°6 2258 11:7 
2190 
43°6 2056 11°5 
45°0 2010 116 
531 
56:2 1666 11:5 


61:2 1477 117 
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APPENDIX A. 


Calculation of Apparent Specisic Heat of Steam. 


Let Q=heat of formation of 1 gm. mol. of steam 
=O 
h=heat-loss before maximum pressure, as fraction 
of Q; 
“= proportion of steam formed ; 


w, y, 2=proportions of steam, H, and (N2+0O,) present 
as diluents ; 


ti, tm=initial and maximum temperatures of explosion ; 


Cyu,=true mean specific heat of hydrogen over the range 
to tin 


Cyy,=true mean specific heat of N, and O, over the range 
t; to tm 


Can,o=apparent mean specific heat of steam over the range 
LO 


Then 
&.Q(1—h)=(tm—t) +y- +2. Cyn, |, 


z FQ(1—-h) — 
tm—ti 


= 


If the original mixture is dry, w=0, and this simplifies to 


The mean specific heats of hydrogen, nitrogen, and 
oxygen are taken from Partington and Shilling’s ‘ Specific 
Heats of Gases’ (1925), and are as follows :— 


0° to 1200° 1600° 2000° 2400° 2800° 
5:27 541 5-69 5:83 
5-22 5°37 B55 5°77 6-02 


It should be pointed out that from 1600° upwards these 
figures have themselves been obtained from explosion 
experiments. It is, however, improbable that they are 
so far incorrect (at any rate relatively to each other) as 
to vitiate the main conclusions reached. 
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APPENDIX B, 
Heat-loss before Maximum Pressure. 


The heat-loss up to the moment of maximum pressure is 
calculated from the formula constructed by Prof, W. T. 
David (Proc. Roy. Soc. A, xeviii. p. 308) from his measure- 
ments of loss by radiation and conduction during explosion 
of coal-gas-air mixtures. It has been assumed that the loss 
from hydrogen-air explosions is appreciably the same, and 
that the loss per cent. per second from two different vessels 
varies directly as the ratio of surface to volume. For the 
18 in. sphere the formula becomes 


Per cent. heat-loss= 


5°57 x (max. temp. ° C.)?° x time of explosion in secs./1000 
L0® x per cent. H,O formed 


For the hydrogen mixtures that have been exploded in the 
45-4 cm. diam. sphere the correction varies from 0:3 per cent. 
for the fastest explosions to 1°4 per cent. for the slowest. 
For Pier and Siegel’s sphere the constant in the formula 
becomes 6°23, for Maxwell and Wheeler’s sphere it is 12°8, 
for Fenning’s cylinder 16°7, for Bone’s sphere 10°8 (at three 
atmospheres, the per cent. loss is assumed inversely pro- 
portional to the initial pressure). 

This correction has been applied, in calculating the 
apparent specific heats, by subtracting it from the heat 
of combustion, which is the only legitimate method. The 
maximum temperatures given are those actually reached, 
as calculated from the maximum pressures. 


APPENDIX C. 


Correction of Pressure Ratios for Initial Temperature. 


This correction is based on the assumption that two 
explosions of the same mixture at slightly different initial 
temperatures, but otherwise under identical conditions, will 
give the same temperature rise, 7. e. for two explosions 
a and b (if T;, T,,=initial and maximum temps. abs.), 


= 
Tia Tia 
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and, if P,,, P; be the maximum and initial pressures and e 
the contraction ratio, 


Pas Pie Tis Ts 


All the pressure ratios are corrected to 15° 1. 


P; J se 288 ef 288 


This correction neglects the effect of the initial temperature 
on the mechanism of combustion. Higher initial tem- 
peratures probably assist the combustion, but within the 
limits of these experiments (10° and 23°C.) the effect must 
be negligible. 


LXXXIX. The Explosion of Hydrogen-Air Mixtures in a 
Closed Vessel—Part II. By B. H. Tuorp, M.Sc.* 


| the preceding paper (p. 813) it was shown that the 

maximum pressure developed in a gaseous explosion. 
is greatly affected by the size and shape of the explosion 
vessel, large spheres with central ignition giving the 
highest pressures. It was suggested that combustion lags 
behind inflammation, and that in the large spheres com- 
bustion has, on an average through the vessel, reached a 
more advanced stage at the moment of maximum pressure, 
and that therefore a larger proportion of the total energy 
is making itself felt as pressure. 

A new diaphragm indicator was also described, and it 
was shown that the maximum pressures recorded by this 
for any particular mixtures in a sphere 45-4 cm. diameter 
were among the highest found by any workers. Pier and 
Siegel alone give maximum pressures from 1 to 2 per cent. 
higher, in a sphere of much the same size, but reasons 
were mentioned for believing that the results obtained in 
this laboratory are closer to the true absolute values than 
theirs. 

Our results will now be considered in detail. They 
consist of six series of explosions at an initial pressure of 


* Communicated by Prof. W. T. David, Se.D. 
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1 atmosphere. In each series the ratio of nitrogen to 
oxygen was kept constant, and the resulting “air ”’ was 
mixed with hydrogen over as wide a range as possible. 
The ratios of nitrogen to oxygen in the six series were 
0, 2, 3-9, 4-9, 6-1 and 6:8. The results are given in 
Table II. 


Taste I. 


The Pressures developed in the Explosion of x per cent. 
H,-+y per cent. O,-++ay per cent. N, Mixtures. 


Time of 
H,O Initial Maximum Explosion in 
/o 2" formed. P; |15°. Temp.°C. Temp. °C. secs./1000 from 


Spark. Rise. 


N,/0,=0. 
150 150 5-46 17:6 1430 12-9 51 36 
189 18-9 6-33 16-0 1742 12-5 fgg 
a 6-76 14-1 1903 12-4 2 19 
22-7 92-7 707 18-5 2027 12:3 15 12 
25:3 25-3 7-49 16-9 2200 12-3 15 ll 
865 26-1 8-13 16-0 2420 10-5 = 6 
87-5 24d 7-84 13-1 2292 10-1 9 7 
88-5 7-52 16-5 2165 9-7 8 
89-4 7-22 14.2 2043 9-4 10 
90-1 19-1 6-91 18-2 1930 9-1 162 
906 18-0 6-71 17:8 1857 8-9 
912 16-9 6-45 14-7 1757 8-7 22 «16 
918 15-6 6-11 17.0 1640 8-5 
92-4 14-4 5-79 18-8 1530 8-3 40-28 

N,/0,=2. 
163 163 5-95 20:5 1600 11-2 41 30 
184 18-4 6-44 16-3 1772 111 34 
203 203 6-82 18-2 1917 11-2 23 16 
22-5 22-5 7-22 14:8 2067 11-4 16 12 
25:2 26-2 7-63 21-4 11-6 12 9 
7-94 12:8 2372 11:8 11 8 
593 273 8-19 17-5 2463 10-9 8 6 
616 258 8-02 16-4 2383 10-6 9 7 
8 
9 
10 
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I. (cont.). 


The Pressures developed in the Explosion of x per cent. 
H,+y per cent. O,+ay per cent. N, Mixtures. 


Time of 

H,O Initial Maximum Explosion in 

2" formed. P;_}150. Temp.°C. Temp.°C. An.c. secs./1000 from 
Spark. Rise. 

N,/O,=2 (cont.). 
69-4 20:5 717 15:8 2030 9-7 15 12 
72-5 18-4 6-72 14:5 1860 9-4 18 14 
74:6 17-0 6-38 13-6 1737 9-2 22 18 
76-5 15°8 6-10 16-7 1637 9-1 32 22 
78-6 14-4 5-77 18-5 1520 8-9 43 30 
N,/0,=3 9. 
15-6 16-6 6-11 14-2 1644 10:7 51 35 
17-4 17-4 6-26 18-6 1710 10:7 44 30 
20-3 20-3 6-88 15-6 1938 10-8 29 20 
22-3 22:3 7:24 18-7 2083 11-0 22 16 
25-0 25-0 7-65 18-9 2250 11-4 19 14 
27-6 27-6 8-01 16-0 2407 11:7 15 12 
29-8 28-9 8-13 20-2 2473 11-9 14 11 
32:5 27:8 8-21 16-6 2476 11-0 12 10 
35-6 26-6 8-13 13-1 2427 10-7 12 10 
38-0 25-6 7-98 17-0 2364 10:5 12 10 
41-0 24:3 7:82 12:3 2290 10:3 12 10 
44:3 23-1 7:59 17-5 2204 10:1 13 10 
49-9 20-7 717 15:3 2030 9:8 16 13 
51:8 19-9 7-00 17:9 1972 9-7 18 14 
53:8 19-0 6-81 16:8 1900 9-6 21 16 
56-7 17-9 6-56 16:8 1808 9-5 23 19 
59-8 16-6 6:26 14:6 1694 9-4 30 24 
62-1 15-7 6-04 18-5 1620 9:3 84 28 
64:8 14:5 5-74 18-2 1516 9-2 42 32 
N,/0,=4:°9. 

12:8 12:8 5:13 15-4 1307 10:9 132 95 
15-1 15-1 5-73 20-4 1517 10:8 70 50 
17-1 17:1 6:20 17-1 1680 10:7 50 35 
19-1 19-1 6-65 17-1 1847 10:7 33 28 
21-9 21-9 7-21 18-6 2060 10:9 27 22 
24:3 24:3 7-64 13:8 2228 11:0 22 20 
25-9 25-1 7°74 16-0 2277 11-1 23 18 
27-5 24-6 7-82 - 158 2295 10-4 19 17 
28-4 24-2 7:78 14:9 2273 10:3 20 17 


31-4 23-2 7-63 16-6 2214 10-1 19 17 
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TABLE I. (cont.). 


Time of 
of H,O Initiai Maximum Explosion in 
Tormed. P; J15°, Temp.°C. Temp.°C. “Ano. secs./1000 from 

Spark. Rise. 
N,O,=4'9 (cont.). 
35:8 21-7 7:33 17:7 2097 10-0 21 18 
39-7 20:3 7:10 12-9 2003 9-8 23 19 
42:5 19-4 6-90 15-4 1926 9-7 24 20 
46.0 18-1 6-64 18-1 1833 9-6 27 22 
50-9 16-5 6-25 13:3 1685 9-4 44 31 
53-7 15-6 6-02 12-3 1605 9-4 pas 40 
N,/0,=6:1. 
13:3 13:3 5-22 19-4 1343 10-8 123 83 
15-7 15-7 5-86 19-9 1564 10:7 72 52 
17:8 17:8 6°35 16-3 1737 10-7 54 40 
19-9 19-9 6-82 18-0 1910 10-7 42 30 
22-4 21-8 7:20 18-9 2057 10:8 32 27 
24-9 21:1 7:16 15-7 2035 10-2 30 25 
20:5 7:04 17:8 1990 10-0 29 22 
29:8 19-7 6-92 19-8 1948 9-9 29 23 
34:9 18-4 6-66 15-2 1842 9-7 32 25 
39-9 17-1 6-36 16-7 1734 9-6 34 28 
44-6 158 6-04 17-1 1618 9-6 43 34 
49-6 14:3 5-65 19-0 1485 9-5 58 45 
§2-7 13:5 5-43 18-1 1407 9-5 75, 56 
54:5 13-0 5-30 17:3 1364 9-5 86 64 
N,/0,=68. 
14:0 14:0 5-46 16-2 1420 10-6 107 80 
14:8 14:8 5-64 16-5 1485 10-6 96 70 
15-7 15-7 5-89 17-4 1570 10-6 72 53 
17:2 17-2 6:23 18-0 1696 10-6 62 48 
17-9 17-8 6-40 18-5 1757 10-6 59 45 
20-2 20-2 6-93 14-4 1943 10-6 43 35 
22:3 20-0 6-95 16-9 1957 10-1 39 32 
24-6 19-5 6-86 17:8 1920 9:9 38 30 
26-2 19-0 6-77 19-4 1890 9-8 36 29 
29-0 18:3 6-63 11-0 1827 9-7 35 28 
33:1 17-2 6-39 14-5 1740 9-6 37 30 
36-4 16:3 6-18 17-6 1667 9-6 42 35 
40-2 154 5-95 15-4 1584 95 55 45 


44-4 14:3 5-66 17:8 1487 9-5 75 60 


Maximum Pressure 


Apparent Mean Specific 


Heat of Steam between I5°C.& 


Inftial Pressure 
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The ratios of maximum to initial pressure (P,,,/P;) are 
all corrected to an initial temperature of 15°C., and the 
“apparent mean specific heat ” of steam (C,,,.,) between 


Fig. 2. 
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Temperature t°C. 
the initial and maximum temperature is calculated for each 


explosion (see Appendices to previous paper). These 
apparent specific heats are shown plotted against the 
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maximum temperatures in fig. 2, curves being drawn 
through the points belonging to each series. There are 
three striking features :— 


1. It will be noticed that each curve consists of two 
almost parallel lines, the points on the lower line being 
obtained from mixtures containing an excess of hydrogen, 
and those on the upper line from mixtures containing an 
excess of “ air.” 

2. In the case of mixtures containing excess hydrogen, 
the replacement of part of that hydrogen by nitrogen 
causes a great increase in the apparent specific heat of 
steam, or, in other words, reduces the maximum pressure. 


3. In the case of mixtures containing excess oxygen, 
the replacement of part of that oxygen by nitrogen 
causes a great decrease in the apparent specific heat of 
steam, 7. e. it increases the maximum pressure. 


Considering 1, three possible explanations might be put 
forward :— 


(a) That there is greater dissociation of steam at 
maximum pressure in the excess oxygen than in 
the excess hydrogen mixtures ; or, 


(6) That the chemical process H,+30.->H,O +57290 
cals. is further from completion at the moment of 
maximum pressure in the excess oxygen than in 
the excess hydrogen mixtures ; or, 


(c) That, contrary to general belief, the specific heat of 
oxygen is much higher an that of hydrogen at 
high temperatures. 


As mixtures containing very little excess oxygen or 
hydrogen differ nearly as much as those containing large 
excesses (c) may be ruled out. For the series in which the 
nitrogen/oxygen ratio is 6:8 the highest temperatures 
reached are only about 1950°C., and since “ dissociation ”’ 
in the strict sense of the word is negligibly small at tem- 
peratures below 2000°C. this cannot account for the 
missing energy. It seems probable therefore that there is 
extensive ‘‘ incomplete combustion (defined in the sense 

Phil. Mag. 8.7. Vol. 8. No. 53. Suppl. Dec.1929. 3K 


830 Mr. B. H. Thorp on the Explosion of 


described in (6) above) at the moment of maximum 
pressure in all the mixtures containing excess oxygen, and 
that at the higher temperatures (over 2000°C.) dissociation 
occurs as well. 

The second feature may be readily explained on the 
principle of mass action, but the third is directly opposed 
to that principle, and apparently some process is taking 
place which completely swamps the mass action effect. 
Comparing for example the mixtures :— 


20 per cent. (H,+40,)+ 70 per cent.”O, giving a maximum pressure of 6-55 
atm., 


20 per cent. (H,+40,)+13-3 per cent. 0,+-56-7 per cent. N, giving @ 
maximum pressure of 6-76 atm., 


20 per cent. (H,+40,)+1-3 per cent. O0,+68-7 per cent. N, giving a 
maximum pressure of 6-82 atm., 


it is evident that here the presence of nitrogen has a 
beneficial effect on the overall process of combustion. 

It is possible that H.O, is a primary product in those 
explosions where there is a large excess of O.; it would 
reduce the maximum pressure both through having a 
higher heat of formation than steam, and also through 
increasing the chemical contraction. 

It has been felt that some attempt should be made to 
find out whether chemical equlibrium is established at the 
moment of maximum pressure. If the exploded gases 
could be very rapidly cooled just after maximum pressure 
has been reached, analysis might show whether chemical 
reaction is still in progress. 

An arrangement has been made whereby a valve in the 
wall of the explosion vessel may be suddenly opened and 
closed again at any desired points during the explosion. 
Through this valve carbon dioxide may be injected in 
large quantities into the exploded gases, cooling them down 
suddenly to 200° or 300°C. Alternatively the valve may 
be opened to allow a large fraction (20 per cent.) of the hot 
gases to be injected into cool carbon dioxide or any 
desired reagent. This apparatus is still in the trial stage 
and no definite results can yet be given. 

It is proposed to continue these experiments at higher 
initial densities, and a start has already been made with 
explosions at three atmospheres. 
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Calculation of the Extent of Incomplete Combustion 
for any mizture. 


If several mixtures which develop the same maximum 
explosion temperature give different values for the 
apparent mean specific heat of steam between 15°C. and 
that temperature, it may be supposed that there is least 
incomplete combustion at the moment of maximum 
pressure in the mixture which gives the lowest value for 
the mean specific heat. It is possible to calculate how 
much more incomplete combustion there is in the other 
mixtures than in that one, but not to obtain absolute 
values from the data here available. 

Fig. 2 shows that the lowest values for the apparent 
specific heat of steam come from mixtures containing 
hydrogen only as diluent. It will be assumed for the 
present that these are the true values, and on that basis 
the amount of incomplete combustion in any other 
mixture may be calculated as follows :— 


Let Q=the heat of formation of 1 gm.mol. of H,O =57290 
Gale. 


a=proportion of incomplete combustion at maximum 
pressure, 


x=proportion of H,O formed by the end of the 
explosion, 


t;, tm =initial and maximum temperatures of explosion, 


Cg,o=true mean specific heat of steam over the range 
to 


=true mean specific heat of the diluent N., H, & O, 
between & ty». 


Then +(1—32/2).Cq] 


2—3 


l—a 


bn—ti 2—3x 


3K 2 
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It must be observed that the calculated maximum 
temperature 


P,, 1 
Ps 


depends on the proportion of incomplete combustion. For 
the case x=0-2, t,, will be 0-6 per cent. lower when the 
incomplete combustion @a=0-05 than when a=0, 

erhaps 12° in 2000°C. This will only be so if ‘‘ incomplete 
combustion ” means that some of the H, and O, molecules 
that should have combined to form H,0O still exist as H, 
and O, molecules at maximum pressure, so that chemical 
contraction is also incomplete. If, however, incomplete 
combustion is not chemical in nature, or if it means that 
some intermediate endothermic product such as H,0, is 
in existence at maximum pressure, the chemical contraction 
may be complete by this time. It has therefore been 
considered advisable to neglect the possible effect on the 
maximum temperature. 


IT. 
Percentage Incomplete Combustion at Maximum 
Pressure. 

Temp.°C. = 2 3-9 4-9 6-1 6-8) 
Excess Hydrogen Mixtures. 
$600 16 2-2 2-7 3-1 3-1 
1800 0 18 2-1 2-2 2-7 3-1 
2000 0 0-9 16 1-7 2-8 
2200 0 0:7 1-2 1-2 
2400 0 0-7 0-7 — = = 


Excess Air’* Mixtures. 


74 6-6 6-6 6-6 58. 
1800S 72 6-0 6-0 6:0 5:3 
2000... 7-1 5-6 5-4 5.0 


These are minimum figures in that they are based on the 
assumption that there is no incomplete combustion in 
mixtures of oxygen with excess hydrogen. 


SC 
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Conclusions. 


It has been shown from the experiments carried out in 
this laboratory on mixtures of hydrogen, oxygen and 
nitrogen, that :— 


(a) There seems to be extensive incomplete combustion 
at the moment of maximum pressure in mixtures con- 
taining excess oxygen. 


(6) Probably the same is true in a lesser degree for 
mixtures containing excess hydrogen ; hence the specific 
heat of steam, as calculated from the maximum pressures 
developed in explosion experiments, must always be rather 
too high. 


(c) By taking the lowest apparent specific heat of steam 
found for a given temperature as the true specific heat, 
the extent of incomplete combustion at maximum pressure 
in other mixtures giving that samemaximum temperature 
may be estimated. 


(d) For mixtures containing excess hydrogen, the 
gradual replacement of that hydrogen by nitrogen causes 
a steady decrease in the maximum pressure, 7. e., an 
increase in the amount of incomplete combustion. 


(¢) For mixtures containing excess oxygen, the gradual 
replacement of that oxygen by nitrogen causes a steady 
increase in the maximum pressure, 7. e. a decrease in the 
amount of incomplete combustion. 


I should like to express my indebtedness and gratitude 
to Professor David for his invaluable direction and advice 
during the progress of this work, and for permission to use 
the results of other workers in this laboratory of the 
Engineering Department, Leeds University (see Part I, 
Table 2; they have also taken a considerable part in the 
development of the new diaphragm indicator). 

I should also like to thank the Department of Scientific 
and Industrial Research for maintenance grants during 
the two years 1925-27, which enabled me to begin this 
work. 
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XO. The Characteristic Numbers of the Mathieu Equation with 
Purely Imaginary Parameter. By H. P. MULHOLLAND 
and 8. GOLDSTEIN *. 


N ETHODS are known for calculating the characteristic 
numbers of integral order of the Mathieu equation 


d’y 
da? 


and tables have been given t for real g. We have calcu- 
lated the first eight characteristic numbers for purely 
imaginary g,and the results are given in Tables I. a and 1.6 
(p. 839). 

The notation is the usual one t. a, and 8, correspond to 
the functions ce, and se, respectively. 

The equatiou with purely imaginary parameter is of 
importance in the physical problem of the alternating flow 
of electricity along conducting elliptic cylinders. In 
addition, it has considerable theoretical interest. 


Numerical Example and Tabular Values. 


As an example we shall give the calculation of @, (corre- 
sponding to the function se,(v),odd in w and $7—a) for 


I= 1°82, 


If == Bey sin 2ra, . . 
r=1 


then the recurrence relations fer the B lead te the formula 
2s 1 5) 
Up—1 = (r> ( ) 

and the equation 


* Communicated by Prof. E. L. Ince. 

t Ince, Proc. Roy. Soc. Edin. xlvi. pp. 20-29 (1925) ; xlvi. pp. 816-822 
(1926); xlvii. pp. 294-801 (1927). Goldstein, Trans. Camb. Phil. Soc. 
xxiii. pp. 803-836 (1927). 

{ See the papers cited above, or Whittaker and Watson’s ‘ Modern 
Analysis,’ Chap. xix. The parameters a=4aand 6=8 are also employed 
sometimes 
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The method of calculation is to take u,=0 fora sufficiently 
large value of 7, and to find the required zero of L by 
Newton’s rule. A first approximation, a, to the zero is found 
by methods to be given later, and a second approximation is 
given by a+6z, where 


To calculate 0L/d we first differentiate (5), obtaining 


2. (2su,_,)=u?,-1 |, (8) 


by repeated applications of which we arrive at the convenient 
formula : 


For q equal to 1°87, we choose 3°6829—3°1779 as first 
approximation to P, (see p. 838). It is sufficient to take ug 
equal to zero. 


Then 
La 3°6 
~ 60°31714+3°1779 


= —0°003:, 


3°6 
“6 45°3171 + 3°1779i + 3°6(0°060 — 00032) 
= 0°0787 —0-0055i, 


3°6 
= 39-3171 + 3°17791 + 3'6(0°0787 —0-00551i) 


=0°10949—0°010602, 


and so on until we get 
ug=0°4702 3365 —0°2213 59312, 
3°6 
+ 3°6(0'4702 3365—0°2213 
= 0°7452 61814—0°8828 483427 


and 
L=1—a+42su,=0°0000 4253 —0:0003 54037. 
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Also 

—OL/da=1— — + ... =0°9380 + 1:13744, 
Hence we take 


0:0000 4253 —0:0003 5403: 
_ 970000 428 7—0-000175i. 
ba 0-9380 + 1°1374i 0000187 


Our next approximation to is therefore 
a= 3'682733 + 3°178075:, 
and with this value of @ we find 
= —0°0000 0002 —0:0000 00012. 


It is to be noticed that the ratios of the coefficients in the 
Fourier series are found in the process of finding the charac- 
teristic numbers, so that the calculation of the functions 
themselves is a simple matter. 

The series of Mathieu * suffice to show that, for sufficiently 
small s, a, and By», arereal. We find that a and «, become 
equal when s=0°1836.... For larger values of s, up to 2:0, 
the tables show a and a as conjugate complex numbers. 
It will be seen later that a certain asymptotic expansion 
gives good numerical approximations to a for the larger 
values of s here considered, and unless this ceases to hold 
with increasing s, a remains complex, and a and a, remain 
conjugate. 

Similar remarks apply to 8, and (4. 

We have found and a, and and to behave in a 
similar way. The numerical values suggest the theorem 
that «, as a function of gq, has an infinite number of branch- 
points on the imaginary axis in the g plane. 

In the neighbourhood of the point where # and e«, become 
equal, Newton’s method fails. Hence, having taken u, as 
zero for some sufficiently large value of 7, we must write 
out the equation L=0 in full, and use another method for 
finding the roots in question—say Horner’s method for real 
roots and the root-squaring method of Dandelin, Lobachevsky, 
and Graeffe for complex ones. [The process was employed 
to find «) and a, for s=0°2. It was sufficient to solve a 
quintic, and it was quicker to find the real roots by Horner’s 
method and divide out than to use the root-squaring 
method. | 


ADA Bom41 are conjugate complex numbers. 


* Journ. de Math, (2) xiii. pp. 187-208 (1868), 
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Asymptotic Expansions. 


Numerical values show that, as s increases, a and ay 
approach «, and and and approach a; and 
« Now a and a, being given on the real axis, their values 
on the imaginary axis beyond the branch-point will depend 
on the path taken from a point on the real axis (or, alter- 
natively, on the cuts made in the g-plane). A full discussion 
of the singularities of a as a function of g is beyond the scope 
of this paper, but for the sake of clarity it is desirable to be 
able to distinguish between a and a. We shall therefore 
suppose the path in the g-plane so chosen that along the 
imaginary axis, for the sufficiently large values of s to be 
considered, # is nearly equal to 6. Similarly we shall 
suppose that a, 1s nearly equal to A. 

Now when gq is real and positive, 


am ~ tad An*, (10) 
where 


m'? +1 m'(m'? +3) 


Am = — 


95 
m! (33m'4 + 410m’ + 405) 
63m’>+ 1260m’* + 294.3m'2 + 486 
m! (527m’> + 15617m'* + 69001m’? + 41607) 
(11) 
Here k=+(82q)?, and m/=2m+1. . . (12) 
Also + 
am ~ Bn . . (1 3) 
where 


* Ince, Joe. cit.; Goldstein, loc. cit. 

+ Goldstein, Proc. Roy. Soc. Edin. xlix. pp. 210-223 (1929), The 
numerical results given here were obtained about the same time as the 
above formula, and provided at first a strong incentive to finding the 
formula, and later a valuable check on its accuracy, since it was obtained 
by a formal process only. 
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and it has been found empirically from numerical results. 
for real g that when m is 0, ¢c; is —0°45, and when m is 1, 
is — 1698. 

We had not proceeded far with the calculation when we 
noticed that the following asymptotic expressions would 
give good numerical approximations when g is a positive 
imaginary : 


where arg qg is given its principal value in interpreting the 
g@ involved. 
We found similarly that the formule 


held for the larger values of s considered. 
(The obvious extensions of these formule are 


an 


Ban +1 44m Bons Ban — Bom—1 (20) 


When q is real and positive, A,, is a detter approximation 
to 4(am+8m41) than to either a, or Numerical 
values show similar results for Aj and A, when g is 
imaginary, so that we have compared Ay with }(a)+,) and 
A, with 4(a3+,) (see Table II. a, p. 840). Also 8,—2% is 
compared with Bo, and 6,—a 3 with B, (Table II. 6, p. 840). 
The values of c, are the same as for real g for m equal to 
0 and 1 respectively. 

Finally, we point out how an approximate value of @ is 
found to start the main calculation. First, when s is suffi- 
ciently small, we use the series of Mathieu previously 
mentioned. When s is moderately large, we take Ap—$Bo 
as an approximation to %, and when a has been calculated, 
take %)+ By as an approximation to Similarly we take 
A,—4B, as an approximation to a3, and a;+B, as an 
approximation to ®,. In this way the approximation to By 
(which is conjugate to 8,) when s is 1°8, namely 3°6829 
—3'1779i, was obtained. As examples we may note that 
when s is as low as 0°4 we get 0°603—0°9907 as an approxi- 
mation to %, and 0°5286—0°93952 as an approximation to 
8,—both very good approximations indeed, with which we 
have no difficulty in completing the calculation. 
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0°8 


0-2 
0-4 
0°6 


10 
1:2 
14 
16 
18 
2°0 


TABLE I. a. 
Se ag. 

0:02 0:003209 
0-04 0:012947 
"0:06 0-029567 
0:08 0:053740 
0:10 0086618 
0-12 0°130223 
014 (188444 
016 0:270636 
018 0°422320 

I. 6, 


0°526248—0°217463 1 
0:599837 —0:992788 i 
0°705302— 1°640667 
0:821211 —2:306374 i 
0-931317 —2°995455 
1:030380—3-699973 
1°119646 —4-413499 i 
1:201620—5-132870 i 
1278205 — 5°856723 7 
1350591 —6°584391 

Age 
2211411 —0-016747 i 
9:130732 —0°141334 
2143946 —0:4068231i 
2+294532 —0°719887 i 


2°541327 — 1065650 
2°852994 —1°475122 7 
3°181293 —1:974374 
3°472882 —2°543309 
3°718564 —3'138044 7 
3°935301—3:737708 


ay. 
0°997324 
0°989185 
0°975226 
0°954772 
0°926663 
0°888869 
0°837489 
0°763155 
0°620331 


Bi: 

0328562 —0°416728 
0°528852 — 0940749 
0°713982—1-602207 ¢ 
0°'838783—2°299441 

0939893 — 2999885 
1:032031—3°704884 ¢ 
1:118567 — 44161827 
1-200182 — 51337154 
1:277220 — 5°856638 7 
1:350120—6°584019 ¢ 


Ba. By 
1:053942 3°978050 
1:223940 3'903906 
1544119 3°742990 
2°162533 3346121 

Bu 


2°895037 —0:906426 
3063353 — 1538181 7 
3'255211 —2:088545 
8464373 — 2°6279197 
3'°682733 —3:178075 2 
3901859 3°7462157 
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Notn.—, and a, are conjugate to a, and B, respectively, and ao and B, 
when complex, are conjugate to #, and 8, respectively. 


0°4 
06 
= 
1:0 
1:2 
1:4 
16 
1°8 
2:0 
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Ila. 

Ao. 

0'2 0°427405—0°317095 0°381560 —0°346229 7 
0'4 0°564544 —0°966768 
06 0°709642 — 1°621437 0°709697 —1°622127 ¢ 
0:8 0°829997 — 2°302907 z 0°829848 — 2°302867 
1:0 0935605 — 2°9976702 0:935618—2°9976347 
1:2 1:031205— 8°702428 
1119106—4-4148402 1119107 —4:414843 2 
16 1200901 —5'138292 7 1-200900 —5°1332927 
18 1277712 —5°8566807 1:277712—5°856680 
2:0 1°350356 — 6584205 7 1:350356 — 6°584205 7 

4). Aj. 

1:0 2°7182 — 0:9860 2 2°6720 —0°9759 

1:2 2°9582 — 1°5067 92-9593 —1:50147 

14 3°2183 —2:03157 8°2234 --2:0809 

16 3°4686 — 2°5856 7 3°4693 —2°0874 7 

1:8 3°7006 — 3'1581 3:'7001—3'15857 

2:0 3°9186 — 3°74207 3:9183—3°7418 7 

TABLE II. d. 

B,—ao. Bo. 
— 0:20 —0'202 —O0'14 —O177 
0°4 —0'0710 +0°5207 —0:0712 +0°05832 
06 +0°0087 +0:0089 +0°03827 
08 +0°01757 +0°00693 +0°01754 +0:00686 7 
1:0 +0:008576 — 0:0044307 +0°008542 — 0:0044392 
1:2 +0:001651 —0°0049112 +0:001640 0:004902 
14 — 0001079 —0:002683 —0-:001081—0°002678 2 
16 —0:001488 —0:000845 — 0:001486 — 0°000842 
1:8 — 0:000985 +0-000085 —0'0009838 +-0:000085 
20 — 0'000471 +0:0003872 —0:000470+0:000372 

Ss. B, 

1-0 +0°35 +0:167 +0187 
1:2 +0°210 —0:063% +0°215 —0:062¢ 
14 +0:0739 —0:11427 +0:0731 —0°11497 
16 —0:0085 —0:08467 —0:0082 —0-0842 7 
18 —0°L8583— 004003 —0:03568 — 0 03988 ¢ 


2°0 —0'03344 —0:00851 ¢ —0-03336—0:00849 
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XCI. Measurement 07 the Flow of Heat. By A. F. Durron, 
M.A., DI.C., and W. G. Marry *. 


[Plates XV. & XVI] 


method recently described + of recording the flow 

of heat across the surface of a wall is only applic- 
able when the flow of heat is toward the wall. It requires, 
moreover, the insertion of thermocouples in the surface of 
the wall, and cannot be used at a glass surface such as a 
window. Simpler apparatus has now been devised which 
does not suffer from these limitations and the transfer of 
heat by radiation and by convection can be recorded as 
separate items. 

2. The instrument shown in PI. XV. is designed to measure 
radiation. It comprises 40 sheets of tinfoil (0°5 in. x 0°3 in.) 
connected in series by wire alternately of 40 s.w.g. constantan 
Zin. long and of 46 s.w.g. copper. Alternate sheets are 
cemented together, with paper insulation in the overlap, to 
form panels (2 in. x 4 in.) and the terminals are connected 
to a recording galvanometer. The obverse of one panel and 
the reverse of the other are sprayed with carbon black. 

3. The radiation instrument was calibrated upon a wall, the 
flow of heat to which was measured by the method referred 
to above. The heating was contrived so that the temperature 
of the air was automatically kept equal to that of the wall 
surface : the heating was purely radiant. 

An E.M.F. of 743 microvolts was recorded by the instru- 
ment when the flow of heat was 10 B.Th.U. per square foot 
per hour. 

This calibration was checked by exposing the instrument 
between blackened copper surfaces, one of which was at 
room temperature. Table I. shows the measured values 
together with the radiation between surfaces of emissivity 
0-92 caleulated from Stefan’s Law (o=5°7 x 107° erg cm.~? 
sec. deg.~*). 

4, The instrument shown in Pl. XVI. is designed to mea- 
sure convection. Itcomprises short lengths of 40 s.w.g. wire, 
alternately of copper and of constantan, soldered together to 
form 16 junctions touching the wall surface and 16 junctions. 
50 mm. from the wall. The terminals are connected to a 
recording galvanometer. 


* Communicated by the Authors. 
t Journ. of Scientific Inst. iv. p. 446 (1927). 
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5. The convection instrument was calibrated upon a wall, 
the convection of heat to which was recorded as the difference 


I, 
‘hemperature Radiation (B.Th.U./sq. ft./hour). 
difference 
Measured. Calculated. 
58 5'6 
4:0 6.6 66 
4:7 81 78 
55 9:2 9:2 
119 
8:0 13°5 13°5 
12°6 21°1 
25°4 
Fig. 1. 
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Calibration of convection instrument. 


between the total flow and the radiant heat measured as 
described above. 


| 
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The measurements show that the convected heat is propor- 
tional to the 1-07th power of the E.M.F. recorded by the 
instrument. Although the relationship is not linear, it is 
clear from fig. 1 that for the range required the error in 
“ making it so”’ is negligible. 

6. The radiation and convection instruments can be com- 
bined to record the total flow of heat. Simultaneous records 
of the flow from a room to a wall by the original method 
and by the new instrument are shown in fig. 2. It will be 


Fig, 2. 
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The flow of heat from a room to a wall simultaneously recorded by 
the original method (above) and by the new instrument, 


seen that the original method ceased to measure when the 
flow became zero and when the uncovered surface of the wall 
became cooler than the surface covered by the square foot of 
cork, On resumption the record was inaccurate for a short 
period until the square foot had regained uniformity with 
the remainder of the wall. 


Note on the Convection of Heat. 


In the calibration of the convection instrument the con- 
vection of heat to a wall was recorded as the difference 
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between the total flow of heat and the radiant heat. A record 
was also made of the temperature difference between the air 
(at 5-0 cm. from the wall) and the wall surface. Corre- 
sponding values are plotted in fig. 3. 

Previous experimenters, using small bodies, have found the 
convection of heat per unit area to be approximately propor- 
tional to the 5/4th power of the temperature difference and 
to be a function of the size of the body. 


Fig. 3. 
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Heat transfer by convection, 


Considerations of Similitude * show that if the convection 
per unit area were independent of the size of the body it 
would be proportional to the 4/3rd power of the temperature 
difference. 

The present measurements show the convection of heat to 
a wall to be proportional to the 4/3rd power of the tempera- 
ture difference 2nd equal to 0°315 B.Th.U. per sq. ft. per 
hour for a temperature difference of 1° Fahrenheit. 


Building Research Station, 
Garston, Herts. 


* Dict. of Applied Physics, i. p. 477. 
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XO. The K-Absorption Disvontinuities of Manganous and 
Chromate Ions. By Don M. Yosr, Fellow of the National 
Research Council and International Education Board *. 


‘T is known that the wave-length of the X-ray absorption 
discontinuity of an element is influenced not only by the 
valence of the element, but also by the nature of the erystal 
in which the element occurs. ‘Thus the wave-length of the 
K-absorption discontinuity of chloride ion in crystalline 
sodium chloride is different from that found in the case of 
potassium chloridet. It was thought that similar effects 
might be found when the discontinuity for ions in aqueous 
solution was compared with that for the same ions present 
in a crystal. 

The solutions of the salts used were placed in fine quartz 
capillary tubes. The outside diameters of the tubes varied 
from 0°04 mm. to 0°l mm.,, the inside diameter being about 
six-tenths of the outside diameter. To fill the tubes it was 
only necessary to dip one end into the solution, the capillary 
action causing the liquid to rise in them. The longer-filled 
tubes were cut into sections of about one centimetre in 
length, and the ends of the latter carefully sealed with a bit 
of hard wax. ‘The tubes, before being used, were placed in 
a spectrograph, which was then evacuated After some 
fifteen minutes the tubes were examined and the defective 
ones rejected. They were also examined at the end of an 
experiment to make sure that no leaks had developed. 

In an absorption experiment a tested tube was placed 
directly in a V-shaped slit in such a way that the walls of 
the tube were tangent to the walls of the slit. Two small 
bits of soft wax served to hold the tube securely in place. 
The vacuum spectrograph used was of the conventional 
Siegbahn type. A piece of 20-w red-coloured goldbeater’s 
skin was employed as a window between the body of the 
spectrograph and the X-ray tube to prevent the light emitted 
by the filament from fogging the photugraphic plate. Both 
filament and anticathode were of tungsten. ‘The tension 
used was about twice that corresponding to the wave-length 
of the absorption discontinuity ; the electron current varied 
from 25 to 100 milliamperes, and the length of the exposures 
from 4 to 8 hours. The analysing crystal was of calcite. 

* Communicated by the Author. 

+ See Siegbahn’s ‘Spectroscopy of X-rays’ Oxford Univ. Press, (1925). 
Also Otto Stelling, ‘ Ueber den Zusammenhang zwischen chemischer 
Konstitution und K-Rontgenahsorptionsspectra.’ Thesis from (Lund 
University, Sweden, 
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Preliminary experiments showed that for the region of 
longer wave-lengths (3°5 A and greater) the absorption of 
the quartz itself was so great that satisfactory photograms 
could not be obtained in a reasonable length of time. The 
final experiments were carried out with solutions of man- 
ganous chloride and potassium chromate, the ions of interest 
being manganous ion «nd chromate ion. The small amount 
of material actually present in the path of the beam of X-rays 
made it necessary to use concentrations of at least 0°7 molal, 
since with lower concentrations the lack of photographic 
contrast at the absorption edge made it impossible to measure 
the plates with any accuracy. 

In the following table are presented some representative 
results. The values there given are the distances in milli- 
metres between the absorption edges and a standard reference 
line photographed on the same plate. The corresponding 
measurements for the solid salts are also given, and were 
obtained with the same spectrograph. 


Results of the Absorption Measurements. 


Radius of spectrograph : 176°27 mm. 
Analysing crystal: Calcite. 


Substance. Concentration, Std. line. white line 
from std, line. 
from std. line. 
MnOl, ... Solid. FeKa, 2°58 mm. — 
MnCl, ... 0699 m, 2°64 ,, 
K,CrO,... Solid, 8:13 mm. 


Within the limits of experimental error the absorption 
discontinuities of manganous and chromate ions in solution 
have the same wave-lengths as the corresponding ions in 
crystalline salts. It is possible that differences might be 
observed in the region of longer wave-lengths, say for the 
case of chlorides ; the experimental difficulties to be over- 
come are rather formidable, however. It is to be especially 
noted that the fine structure of the chromate edge (desig- 
nated as a white line in the table) observed in the case of the 
solid salts is also observed unchanged in the spectra obtained 
with solutions. Moreover, it was also found that no observ- 
able reduction of the’ chromate to lower valence forms had 
taken place. In this case, at least, the fine structure is to be 
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attributed to the chromate chromium itself, and not to reduced 
forms arising as a result of the action of the X-radiation. 


I wish to express here my appreciation of Prof. M. 
Siegbahn’s kindness.in allowing me to make use of the 
facilities of his laboratory at Uppsala. I am also indebted 
to him for suggesting and providing the capillary tube. 


Uppsala, 
June 1929, 


XCIIL. The Scattering of Beta-Particles by Light Gases and 
the Magnetic Moment of the Electron. By Matcoimm 
C. Henperson, Ph. D., Trinity College, Cambridge *. 


1. Introduction. 


N 1925 Uhblenbeck and Goudsmit + showed that the 
multiplet structure of line spectra could be interpreted 

if it were postulated that the electron is not merely a 
point charge of electricity, but possesses also an angular 


momentum of zs o2 due to its own rotation and a 
T 


magnetic moment equal to one Bohr magneton. This 
postulate provides a theoretical basis for the fourth 
quantum number. The problem then arises as to the 
effect the magnetic moment will have when two electrons 
collide with one another. 

If an electron possesses at all times a magnetic field, 
such as is necessarily associated with a magnetic moment, 
the collisions between it and a f-particle, that is, between 
two electrons, should be profoundly modified. For 
example, it may be calculated simply from the law of 
force between magnets in the ‘“end-on”’ position 
(6MM’ /r*=e?/r?, where M=M’=9-2 Gauss-cms.) 
that the distance at which the electrostatic force equals 
the magnetic force is about 510744 cm. As usually 
calculated from the electrostatic force alone, the closest 
possible approach in a head-on collision between a 
344,000 volt B-ray andan electron is about 8-3 x 10-cm., 

x Communicated by Dr. J. Chadwick. The material in this paper 
formed part of a Thesis submitted for the degree of Doctor of Philosophy 
at the University of Cambridge. 

+ Uhlenbeck &\Guoudsmit, Naturwiss., xiii. p. 953 (1925); ‘Nature,’ 
exyii. p. 264 (1926). 
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that is, about 1/60 as much. At less than 5 x10-“em. 
the magnetic force of attraction will exceed the electro- 
static force of repulsion, since it varies as a higher negative 
power of the distance. The deflexion experienced by a 
f-ray on passing an electron at 60 times the minimum 
possible distance of approach is about half a degree. All 
deflexions larger than this amount should be greatly 
influenced by the magnetic force. The angle through 
which a particle is deflected for a given nearness of 
approach should be widely different from the deflexion as 
calculated from the inverse square law only. Conse- 
quently the manner in which the number of particles 
scattered varies with their deflexion should also be 
different from the calculated distribution. It may be 
noted in passing that an inverse fourth power orbit is 
unstable and under some conditions the electrons might 
theoretically coalesce. 

The other kind of magnetic deflexion, namely, that 
experienced by a charged particle moving in a magnetic 
field, is eliminated by the assumption that the particles 
come at once into the end-on position. The motion is 
then entirely along the lines of force and the equation 
Hp=mwv/e no longer applies. If this assumption is not 
made, a rough calculation shows that the maximum 
distance at which the two forces may be equal is about 
8 

From the theory of single scattering, as developed in 
detail by Darwin*, it can easily be shown that for small 
angles the ratio of the nuclear to the electronic scattering 
is N?: N, where N is the atomic number of the scattering 
element. Hence for hydrogen and helium the electronic 
scattering will be respectively 4 and 4 of the total. For 
these elements a discrepancy in electronic scattering should 
therefore be expected to show most markedly. By com- 
paring the scattering of B-rays by these two gases with the 
scattering by nitrogen and argon it should be possible to 
determine whether the electrons do in fact possess some 
field of force in addition to their normal electrostatic field. 

Of the solid elements in which the scattering of B-rays 
has been studied, aluminium is the lowest in atomic 
number. For small angles the electrons of this element 
cause, theoretically, about 8 per cent. of the total scattering. 


* ©. G. Darwin, Phil. Mag. xxvii. p. 499 (1914). 
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In previous investigations * the interest was in the nuclear 
scattering and a correction was made for the scattering 
by the electrons. The agreement between the experi- 
mental results when. corrected and those expected from 
theory suggests that the electronic scattering is in fact 
not very different from 8 per cent. But the correction is 
too small, compared to the total scattering, for this 
evidence to be conclusive. From observations of f-ray 
tracks in a Wilson expansion chamber, Bothe + found the 
electronic scattering in air to be of about the expected 
amount, but the small number of collisions observed 
leaves something to be desired statistically. 

Accordingly, some preliminary measurements were made 
with hydrogen, helium, nitrogen, and argon in a simple 
form of scattering chamber. Use was made of a familiar 
principle, namely, the reduction in the number of particles 
in a beam caused when some of them are deflected through 
more than a fixed angle. The results suggested an increased 

scattering in hydrogen and helium. The following more 
accurate investigation was therefore undertaken. 


2. Apparatus and Method. 


The scattering chamber used was of the “ annular 
ring’ type. It was similar to the one used by Rutherford 
and Chadwick? for measuring the scattering of «- particles 
in helium. 

The dimensions are shown in the figure (see p. 850). 
The diaphragms were originally of graphite, but the natural 
effect when the chamber was evacuated was much too high. 
The diaphragms were therefore backed with 2 mm. of lead 
(not shown in the figure) and the central disks were joined 
and supported by a brass rod only slightly smaller in 
diameter than the disks. The support for the rod was a 
thin sheet of graphite. The scattering from this support 
is unavoidable. It amounted, however, together with 
other residual scattering, to less than 4 of the total current 
and was deducted from the readings. The average angular 
limits of the scattering were 9-2 and 29-7 degrees, and the 
thickness of the scattering layer was 3 cm. 

The gases measured were hydrogen, helium, Reouen. 
argon, and air. The first four were drawn from cylinders of 

* Chadwick and Mercier, Phil. Mag. 1. p. 208 (1925); and others. 


{ Bothe, Zect. f. Phys. xii. p. 117 (1922), 
{ Rutherford and Chadwick, Phil. Mag. iv. p. 605 (1927). 
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the compressed gas. Every gas passed over phosphorus 
pentoxide before entering the scattering chamber, and the 
hydrogen and helium were in addition passed through two 
charcoal tubes cooled in liquid air. A test experiment 
showed no difference between hydrogen that was thus 
purified and that which was taken directly from the 
cylinders. Nevertheless it was always thus purified. 
The nitrogen was taken directly from the cylinders. 
Any impurity was presumably air, and since the ratio 
of the theoretical scatterings of air and nitrogen is 69-6: 
65-8, an impurity of this kind would be quite negligible. 
The argon used was about 99-5 per cent. pure. 

The source of f-rays was a deposit of radium E on a 
nickel button 3 mm. in diameter. The deposit was 


Mica, 


| 


} 


obtained by stirring the button for some time in a hot 
solution of radium D, E, and F. The initial strength of 
the source was equivalent to about # mg. of radium, as 
measured on a f-ray electroscope. The amount of 
radium D on the button was less than 1 per cent. 

The f-rays from radium E are heterogeneous and their 
velocity varies over rather a large range. The average 
energy is in the neighbourhood of 344,000 volts*, B=0-8. 
The deflexion which any particle experiences in a given 
encounter varies inversely as the square of its energy. 
This fact is not important as long as only the relative 
scattering by different elements is measured. A possible 
source of error is that hydrogen may scatter, for example, 
the low velocity rays relatively more strongly than the 


* Ellis and Wooster, Proc. Roy. Soc. A, exvii. p. 109 (1927). 
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high velocity rays in comparison with argon. Such a 
selective effect would mean that the scattering varies with 
the energy of the rays according to a different law for 
different atomic numbers. For such a departure from 
theory there is no experimental evidence. 

The measurements were made with a Dolezalek electro- 
meter with a sensitivity of about 3000 mm. per volt. 
During each experiment the electrometer was held at its 
zero position by a charge induced on a small condenser 
forming part of the insulated system. A measure of the 
ionization current was afforded by the time required for 
the voltage on the outer casing of the condenser to change 
by a fixed amount while holding the electrometer 
stationary. The ionization chamber was cylindrical, 
10 cm. in diameter and 14 cm. deep. The outer casing 
was maintained at 200 volts. Guard rings protected the 
insulated system from leakage across the insulation of the 
ionization chamber and condenser. 


3. The Theoretical Scattering Power. 


In the original annular scattering chamber* the 
scattering substance was in the form of foil. A term 
t sec /2 entered the integrand of the expression for the 
scattering because the foil was not perpendicular to the 
incident beam of particles. In the chamber as modified 
for gases ¢ is measured along the incident beam and is 
approximately constant. The expression for the scattering 
by the nuclei then becomes : 


me? 2 


ann 
=—.,-|—s) | —cosec? 
167? \mv?, 


v, 1s the number of particles received per unit area 
of the final diaphragm perpendicular to the axis of the 
chamber, Q is the number of particles in the incident 
beam, ” isthe number of atoms of the gas per c.c., N is the 
atomic number of the scatterer, m is the mass and v the 
velocity of the f-particle, and r is the average distance 
from the source to the scattering layer. This expression 
differs only slightly, in numerical value, from the more 
complex formula in the original paper. 


* Chadwick, Phil. Mag. xl. p. 734 (1920). 
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If we consider an electron in the atom to be at rest before 
the encounter and take its mass equal to that of the 
incident f-particle, then the amount of the electronic 
scattering will be given by 


2 
Ve=N | — cosee 
The number of electrons projected between the angles 
and will be 


/ e? \? Pa 
Ar® \ mv? ) [sec ¢| 


For the values of ¢, and ¢, used in these experiments 
the projected electrons are relatively so few in number 
that they will be neglected in the calculations which 
follow. 

The sum of these expressions gives the total scattering 
by the gas. 

However, as shown by C. G. Darwin*, each of these 
expressions must be corrected for the change in the orbit 
caused by the variation of the mass of the f-particle with 
its velocity, the so-called relativity correction. For the 
nuclear sip miss this takes the form 


Be osec where Scotw= tan [ 
2 

The mean value of this correcting factor was calculated 
to be 1-2 for the conditions of the present experiments. 

The corresponding correcting factor for the electronic 
scattering has not been worked out. However, elemen- 
tary considerations show that it will have a value less than 
one, that is, it will serve to reduce the calculated amount 
of electronic scattering, for in the encounter with an electron 
the B-particle will be retarded in its motion, not accelerated 
as in an encounter with the positive nucleus. The correc- 
tion will in general be small and it has been omitted in the 
following calculations. 

Adding the above expressions for the phic: and 
electronic scattering and Pas 


* C. G. Darwin, Phil. Mag. xxv. p. 201 (1918). 
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Substituting the values for the experimental conditions, 


VizZ., 
cot? 
yak 
Oni fe? 


or, to a sufficiently exact approximation. 
y= 35°8k( 1:2 N?4N). 
As shown by this formula the number of scattered 
particles should be proportional to 1-2N?-+-N when there 
are the same number of atoms present per c.c. The values 


entered in column 2 of Table II are calculated from this 
expression. 


4. Experimental Results. 


In each experiment readings of the ionization current 
were taken for different pressures of the various gases: 
hydrogen, helium, nitrogen, argon, and air. When the 
pressure is plotted against the ionization current the 
points lie approximately on straight lines, indicating the 
predominance of single scattering, up to a pressure of about 
2 cm. for argon, 6 cm. for nitrogen, 50 cm. for helium and 
70 cm. for hydrogen. There is a tendency toward upward 
curvature as plural scattering becomes more prominent *. 
The observed values of the ionization current in two 
experiments, corrected for natural effect, with the corre- 
sponding pressures are shown in columns 2 and 3 of 
Table 1. The ionization current divided by the pressure 
gives a constant for each gas, and the ratios of these 
constants to each other are the ratios of the scattering 
powers. ‘The quotients of the figures in column 2 divided 
by those in column 3 are entered in column 4. 

Table II. gives the collected results of the last five 
experiments. In column 2 are the theoretical values of 


* It may be noted that the scattering ceases to be single at about the 
point indicated by Wentzel’s criterion. For example, the scattering in 
nitrogen is single up toa pressure of about 6 cm., for which 4@ under 
the conditions of these experiments was about 4'8°, while @ varies from 
9° to 80°. (Cf Wentzel, Ann. d. Phys. Ixix. p. 335, 1922.) 
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the relative scattering power of the gas concerned (see 
section 3), in column 3 are the theoretical percentage 


I. 


Showing the pressures and corresponding ionization 
currents for various gases and the resulting relative 
scattering powers. 


1 2 3 4 5 
P I L/P. Av 
501 23-3 4-65 
N : 4:23 20-5 4-84 4-64 
3-23 14-4 446 J 4-64 765-3 
73-50 17-2 234 >) 
55-00 12-9 235 4-64 65-8 iN, 
| 37:30 10-1 O71. 
1:90 3879 19-9 19-4 
1:21 23-9 19-7 19-4 203 PA 
2:05 16-9 825 
Air 414 35:8 8-65 8: 
612 50-5 8-25 
II. 


Collected Results. 


1 2 3 4 5 
Theoretical Theoretical Experimental Ratio 
scattering percentage scattering experimental 
Gas. power electronic powers. to 


1-2N?--N contribution. A=203. theoretical. 


2-2 45 4-0 18 
3-4 29 4-8 1-4 
65-8 10:6 79 1-2 
69-6 10-4 84 1-2 
203 4-4 203 1 


contributions by the electrons to that scattering power, 
in column 4 are the experimentally determined scatterings 
all reduced to 203 for argon (the ratio of nitrogen to air is 
taken as normal, that is, equal to 65-8 : 69-6), in column 5 
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are the ratios of experimental to theoretical scattering 
power. 
5. Discussion. 

The results in Table Il. show that the lighter gases 
possess a scattering power for B-rays in excess of that 
calculated from the expression 1:2N?+N. This excess 
scattering may be attributed to the orbital electrons, 
since it increases with decreasing atomic number in about 
the same manner as does the theoretical electronic 
scattering. An attempt was made to fit the expression 
1-2N?--kN to the observed scatterings by using various 
values of k. When k=3-3 the relative values of this 
expression lie within eight per cent. of the observed 
scatterings in the gases measured. 

In the collisions between nucleus and f-particle it is 
usually assumed that the magnitude of the velocity of the 
B-particle is unchanged by the encounter, unless there is 
radiation. On the other hand its velocity is reduced in an 
encounter with another electron because of the energy 
imparted to the latter. Calculation of a special case 
shows that when a f-particle of velocity -8c is scattered 
by an electron through an angle of 11-2 degrees, its velocity 
is reduced to -79c and the projected electron moves with a 
velocity of -25c in a direction at 76-7° with the incident 
particle. The loss of velocity suffered by the f-particle 
is thus small. No correction has therefore been made for 
this effect in these experiments. 

As calculated in section 1, if the magnetic moment of 
the electron is one Bohr magneton, the maximum distance 
at which the magnetic and electric forces between two 
electrons become equal is 5x 10-14 cm. This is about 60 
times the minimum apsidal distance for collisions between 
the average radium E f-particle and an electron, calcula- 
ting the apsidal distance in the usual way for a purely 
electrostatic field. The theoretical deflexion for’ an 
impact parameter of 60 times the minimum apsidal 
distance is about half a degree. For larger angles of 
scattering the magnetic force should be the more important 
factor in determining the final path of the particles. It 
would appear that under the conditions adopted in the 
present experiments the electronic scattering should be 
determined almost completely by the magnetic forces. 
A rough calculation on these lines shows that the amount 
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of scattering should be very different from that actually 
found, and it therefore seems unlikely that the electron 
can have a magnetic moment as large as one Bohr 
magneton. The fact that an increase in the number of 
particles was observed, roughly proportional to the relative 
importance of the electron in the atoms, suggests, how- 
ever, that there is an additional field of force present in 
collisions between electrons. Without measurements over 
other ranges of angle the nature of this field of force cannot 
be inferred. 

From the nature of magnetic force the force between 
two magnetons must depend both in sign and magnitude 
upon the orientation of the axes of the particles with respect 
to each other. The magnetons will also in general exert a 
torque upon each other tending to bring them into the 
position of greatest attraction. If, however, the magnetic 
moment is the result of an actual spin, the two particles 
will merely precess about the line of centres. The force 
between them will not then be increased by their rotation 
toward the end-on position. Calculation of the orbits 
described under such conditions is difficult. Hicks * and 
Wessel t have made calculations of this nature for the 
somewhat simpler conditions of a fixed magneton and a 
non-magnetic moving particle. Bothe { also has calculated 
the orbits of a B-particle about a nucleus, using a somewhat 
different assumption. He considers that the magnetic 
B-particle is moving in the inhomogeneous magnetic field 
of the nucleus. In this case there are two possible 
orientations. The axes of the particles may be either 
parallel or anti-parallel. 

An objection to the conception of an electron as 
possessing spin and consequent magnetic moment is the 
relative complexity of this model compared to the point 
charge. Dirac § has shown how to avoid the necessity for a 
magnetic electron by a more complete solution of the 
fundamental equations. The electron according to his 
theory is still a point charge but it displays a duplex 
character, in many respects resembling the behaviour of a 
magneton. Iwanenko and Landau have shown by an 


* Hicks, Proc. Roy. Suc. A, xe. p. 356 (1914). 

J Wessel, Ann. d. Phys. xxviii. p. 757 (1925). 

{ Bothe, Zeit. f. Phys, xliv. p. 548 (1927). 

§ Dirac, Proc. Roy. Soc. A, exvii. p. 610 (1928): exviii, p. 351 (1928). 
4] Iwanenko and Landau, Zeit. f. Phys. xlviii. p. 840 (1928). 
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entirely different method that the assumptions of wave 
mechanics lead to the same result. The problem of the 
collision of a B-particle with an electron has not yet been 
examined on the basis of Dirac’s equations, and although 
it seems possible that the duplex character of an electron 
in an orbit may be interpreted to give a magnetic field 
effective in collisions, the magnitude of the effect eannot 
be estimated without a laboured calculation. When such 
calculations are available it may be possible to give an 
adequate explanation of the results obtained in the present 
experiments, without any further hypotheses. 

In the experiments performed heretofore on f-ray 
scattering, additional fields of force, usually magnetic, have 
often been postulated to explain increased scattering, but 
they have always been discarded later after more exact 
investigation. Whether the increased scattering found in 
the present work must finally be explained by that postulate 
is a matter for further investigation to decide. 


SUMMARY. 


The hypothesis of the magnetic electron put forward by 
Uhlenbeck and Goudsmit in 1925 raises the question as to 
the effect which the magnetic field of the electron will have 
upon the collisions between f-particles and electrons in 
their orbits. Experiments upon the scattering power of 
light gases for B-particles have therefore been carried out to 
determine whether there is any such effect, 

The experiments with the f-particles from radium E 
show that hydrogen and helium possess a scattering power, 
measured between 10° and 30°, about 80 and 40 per cent., 
respectively, in excess of that to be expected from the 
usua! theory. The variation in scattering power with 
atomic number (N) is fitted within 8 per cent. by the 
expression a@N?+k£N, where a is the relativity correction 
and k is 3-3 instead of the theoretical value slightly less 
than unity. Although this result suggests that the 
electron possesses a field of force in addition to the normal 
electrostatic field, it apppears unlikely that it can have 
a magnetic moment as large as one Bohr magneton. 


The author is indebted to Professor Sir Ernest 
Rutherford and to Dr. J. Chadwick for their advice and 
interest during this investigation. 


XCIV. The thermal Expansion of Liquids according to van 
der Waals. By J. Professor in Physics 
at the University of Ghent *. 


Wee Philosophical Magazine has published recently a 
paper by Mr. V. N. Thatte+, in which the author 

intended to show that from van der Waals’s equation of 

state, 

v—b’ 


(1) 

taken at temperatures so low that the vapour pressure p may 

be neglected against the pressure of cohesion =, a formula 


may be derived for the relative coefficient of cubical expan- 


1 dv 


sion C= oT of the liquid similar to that of Davies. 


Putting p=0 in the equation (1), one finds 


oT, ... 


and this relation is indeed similar to the formula of Davies 

when we suppose, as Mr. Thatte does, that the term * is 
v 


(at least nearly) a constant proportional to the critical 
temperature TT, ; only the coefficients are different. Deter- 
mining the coefficient of T, by using experimental data, 
Mr. Thatte finds 
T,—2T. 

I wish to say that [ do not quite agree with this result. 

I do not see why Mr. Thatte deduces the value of the term 


sR rom experiment, instead of deriving it from the equation 
Uv 
of state, by which it iscompletely determined ; that seems to 
me illogical. Now we know that at very low temperatures 
that equation gives v=b, and it is known also that, according 
to the same equation, 
* Communicated by the Author. 


+ V. N. Thatte, ‘Coefficient of Cubical Expansion of Liquids and 
Critical Temperature,” Phil. Mag. vii. p. 887 (May 1929). 
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so it results from the equation itself that the term ss is 
really proportional to T,, but the coefticient of proportionality 
: 
=3°4 instead of being 2°5. 


But this is only a limit value. One must not overlook the 
a 
fact that the term oR MY not be considered as constant, 


because v is a function of temperature, and this circumstance 

modifies also the value of the coefficient of T. In order to 

find this coefficient, we write equation (2) in the form 

a. T 9 8 T 

(4) 


: 
where 6= ~ is the (molecular) density ; for according to 
v 


van der Waals’s equation. 


Now values of : corresponding to different values of ir 
e 
have been calculated by J. P. Dalton *. From his table and 


equation (4) we deduce the following one :— 


3.00 1 _ 3.375 To. 
To 0c 
01 2:91 3:07 
0-2 2-81 2-76 
0-3 2-70 2-44 
0-4 259 211 
O-5t 246 1-76 


and from this it is easily seen that we may put approxi- 
mately : 


1 


1 
But according to van Waals the relation between and T 


is not really linear, as in Davies’s formula, deduced from the 


* See J. P. Kuenen, ‘Die Zustandsgleichung der Gase und Fliissig- 
keiten, p. 94 (1907). 


+ Up to this value of T the vapour pressure remains negligible. See 
Kuenen, loc. cit. 
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law of the rectilinear diameter. ‘The exact relation is to be 


found by substitution in (4) of the fnueben determined by 


86 


the equation : 


RT or a= 
which may be put in the reduced form : 


The greatest root of this quadratic equation must be taken. 
By developing it in a series, one finds 


..... 
with 
T 
(9) 
and consequently 
or 
8 
C= ....), . . 
277. 


which last relation may be deduced directly from (8) by 
observing that 


dlogd 

The relation (11), however, is less useful than (10) on account 
of its less rapid convergency. 

It is obvious that form (6) does not agree well with the 
experimental data; only Davies’s formula does for normal 
substances at low temperatures (at which the density of the 
vapour is nearly zero). Even Thatte’s formula, partly 
deduced from experiment, is in disagreement, save for 
T=+4T., at which temperature both the formule of Davies 


and Thatite give. = 1:5 It may, moreover, be remarked 
that van der Waals’s equation does not give a quite different 
value. 


Natuurkundig Laboratorium, 
Ghent, Aug. 1929. 


XCV. On the Operational Solution of Linear Differential 
Equations and an Investigation of the Properties of these 
Solutions. By BATH. vAN DER Pon, * 


1. Introduction. 


N a former paper? the rules for solving operationally 
linear differential equations with constant coefficients 
and with an arbitrary second member were considered. 
Arbitrary given initial conditions could be taken care of, 
The best-known case is that for which the system is 
originally fully at rest, while at the instant ¢=0 a unit force 
is suddenly impressed upon the system. 

It is the purpose of the present paper (i.) to consider 
some general rules in connexion with the operation:l 
calculus, (ii.) to outline a method similar to the one used 
in the former paper for the solution of linear differential 
equations with variable coefficients, and (iii.) to consider 
simple operational methods for the investigations of the 
properties of the functions defined by those equations. 

Though the subject matter is relatively new in physics ¢ 
the idea of fractional differentiation and integration, which 
forms an important part of the operational calculus, was 
worked out by Riemann (‘ Werke,’ p. 331), a report of 
which is to be found in Pincherle’s article on “ Funktional- 
operationen und -Gleichungen”’ in the Enz. der Math. Wiss. 
ii. A. 11, where again, however, no mention is made of 
Heaviside’s important contributions to the subject §. 

As in our former paper, we will here again take as a 
basis of the development the integral equation of Carson ||, 
of which Bromwich’s well-known integral is the solution. 


* Communicated by the Author. 

+ Balth. van der Pol, “A Simple Proof and an Extension of 
Heaviside’s Operational Calculus for Invariable Systems,” Phil. Mag, 
vil. p- 11538 (1929). 

{ E.g. we do not find any reference to it either in the new edition 
of Riemann—Weber, ‘ Die Differential-gleichungen der Physik’ (1927) 
or in Courant—Hilbert, ‘ Mathematische Physik,’ i. (1924). 

§ The same is true for the article by A. Voss, “ Integralrechnung,” 
Enz. der Math. Wiss. 11. A. 2, p. 116, where, also, fractional differentiation 
is considered. 

\| J. R. Carson, ‘Electric Circuit Theory and the Operational Calculus 
(McGraw-Hill, New York, 1926). 

gq T.J. Va. Bromwich, Proc. London Math. Soc. xv. p. 401 (1916), 
In this paper arbitrary initial conditions are already considered in the 
operational solution of differential equations. See also H. Jeffreys, 
‘ Operational Methods in Mathematical Physics’ (Cambridge University 
Press, London, 1927). 
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In order to expound our method as clearly as possible we 
insert various examples, mostly taken from Bessel fanctions, 
spherical harmonics, and other functions. It is found to be 
a relatively simple matter, once the operational represen- 
tutions of these functions have been obtained, to write down, 
almost at a glance, many relations between those functions, 
of which several may be new. In particular we refer to 
the parallelism between the familiar relations which exist 
between Bessel functions of different order on the one hand 
and the well-known similar relations between spherical 
harmonics on the other hand. It will, for instance, be shown 
that the one set of relations can directly be derived from 
the other. Further, three new expressions will be derived 
for the derivative of Bessel functions with respect to their 
order. In this respect a new function, 


which is analogous to 


x 


will be shown to play an important part. It looks to the 
writer as if the shortest way to expound and explore 
the properties of Bessel functions, for example, is by the 
powerful method of operators. Yinally, at the end of 
the paper a list is given of the symbolic representations 
of many functions, of which several may here be given for 
the first time. Throughout, the notations of the well-known 
tables of Jahnke-Emde®* will be adhered to. 


2. Some General Theorems. 


We assume the symbolic or operational representation 
f(p) of a function A(z) of « to be defined by Carson’s 
integral : 


Hp) =p erm ha) de. 


Thus, on the one hand, when the function h(x) is given 
its operational representation f(p) is found by an integration. 
On the other hand, when the operational expression /(p) 
of the unknown function h(#) is given, (3) represents an 
integral equation for the unknown function h(x). The 


* Jahnke-Emde, ‘ Funktionentafeln’ (‘'eubner, Leipzig und Berlin), 
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solution of the integral equation is given by the following 
complex integral, first obtained by Bromwich :— 


where ¢c >0; this was pointed out by Lévy * and March t. 
Now from the definition (3) of the operational repre- 
sentation f(p) of a function h(x) a set of theorems has been 
derived by Carson. As we shall make frequent use of 
these relations, they will here be summarized without proof. 
The proofs can be found in the above-mentioned treatise of 


Carson. 


Thus, if = h(z) 
then (?)= h(s#)  (s=constant), . . . (5) 


p) = h(@) (provided h(0)=0), (6) 


for a<Xr 
= hw +2) §, (A>0). (10) 
when h(x)=0 for 0<a<n. 
To this list the following relations may be added :— 


= mo), (n>0), . (11) 


p{- ar (n>0), (12) 


oil 


* P. Lévy, ‘Le Calcul symbolique de Heaviside’ (Gauthier-Villars, 
Paris, 1926). 

+ H. W. March, “The Heaviside Operational Calculus,” Bull. Am, 
Math. Soc. xxxili. p. 311 (1927). 

{ When /(p) is the operational representation of h(x), as defined by 
(3), we use the symbol =, and write f(p) = h(x’). 

§ See also Jeffreys, /. c. p. 18. 

3M 2 
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o Pp 


Lim. f(p) = Lim} A), 
poo 
Lim.f(p) = Lim. A(@), 
p->0 


provided the definite integrals and limits have a meaning, 
i.e. they can be extended to the limits 0 and «0, 


If, further, 
= hi(2) 
Sa P) = ha(2), 


and 


then 


== . . 


provided either /,(x) or ha(w) can be expanded in a series of 
positive powers of «. 
Finally, we have the very important theorem 


A(p)-fa(p) \, (w@—&). dE | 


= h(E) .hy(a—&) 


As the theorems (11)-(16) were not given by Car tl 
proof will be indicated here. 5 7 


To prove (11) differentiate both sides of (5) witl pect 
to sand then put s=1. We obtain 
dp) _. ,, | 
A repeated differentiation with respect to s thus yields (11 
When, further, (11a) is multiplied by 


to (7), means integration on the right-hand side between 
zero and x, we obtain 


df 


_ FP) d T(p) 
dp. ipl ) A(x), 
from which (12) follows. 


(17) 


or 
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To demonstrate (13) an integration of both members of 
(5) can be effected between the limits zero and unity. 
If the same integration is taken between tle limits zero and 
infinity (14) is obtained, while the subtraction of (13) and 
(14) 3 yields (13a). Further, (15) follows from (3) when it 
is remembered that 

Lim, pe-?? = 

2 
is a function of the type as used by Dirac* in atomic theory, 
and called by him a 6 function. It has the property 


+a 
{ dx = 1, 


and is zero everywhere except at «=0, where it is infinite. 
These functions have the property 


i} 8(x). h(x) dx = h(0), 


and thus (15) follows. It is interesting to remark that 
these 6-functions were recognized and used extensively by 
Heaviside ft already thirty-five years ago. He called them 
‘impulsive ” functions. 


Thus 
p= 82), 
p= 
(a), 


with the general property 
+o (n) 
fa). dv = (0). 


Further, (16) can be demonstrated with the aid of (12). 
When one of the functions h,(#) or h(a) can be expanded 
as a power series (8) is thus seen to be a special case of (16), 
because we have 


= p+ 
= 


* See, e.g., Proc. Roy. Soe. A, exili. p. 625 (1926). 
+ See, Heaviside, Theory, ii. p. 92. 


== = h(a), 


° . 
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and (15), we find 


ax (Alt ( 


+a) = he(a). 


Finally (17), which is the most important theorem of all, 
is frequently used by Carson. It is of the type of a ‘produit 
de composition,” as defined by Volterra, and it is on the 
basis of this theorem that Lévy * builds up the operational 
calculus. This expression was also at the base of Riemann’s 
investigationt. Further, this theorem was given by Borelf, 
and is closely related to Borel’s “associated function ”’ § and 
his ‘‘ sommation exponentielle ” of divergent series. More- 
over, it expresses the well-known Hopkinson—Boltzmann’s 
superposition principle, and it is used and demonstrated in 
a special case by Lord Rayleigh in his.‘Theory of Sound,’ i 
p- 74. It will be of frequent use further on. 


3. Hxamples illustrating the General Theorems. 


After having obtained a rather complete set of theorems, 
or working rules, we will now consider some examples of 
application of these rules. From the definition of the II(n) 
functions 


II (n) = 


0 
where II(n)=n! for n= positive integer, 
i follows at once from (3) that 
ge 
whicn s vaiid for all values of n for which n> —1, fractional 


values of included. 
Similarly, we have 


(18) 


* Loe. cit. 
+ Enz. der Math. Wiss. ii. A.11, p. 771. 
t E. Borel, ‘Legons sur les Séries divergentes,’ (Gauthier-Villars, 
Paris, 1928). 
§ See, e.g. Whittaker and Watson, ‘Modern Analysis,’ p, 140. 
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These expressions can easily be verified by series expansion. 
Ei. gs 
p 1 1 


4 


+. = COSW. 


Consider (9). Let again f(p)=hA(z). As (3) only 
considers the values of h(z) in the region 0 to ow, it 
cannot give any information as to the behaviour of h(z) ia 
the region x <0. In connexion with (9) it is, therefore, 
appropriate to consider all our functions to be zero for «<0, 

This is demonstrated in fig. 1. A factor e~” therefore 
shifts h(a) over a distance to the right, while a factor e+’? 
shifts it over the same distance to the left, but in the last 
case a part of the function is moved through an imaginary 
vertical cut along the y-axis and slipped under a rigid cover, 
which keeps the function equal to zero whenever «<0. 
That this is the only consistent interpretation follows directly 
from the integral (3). For 


0 


= rf h(x) da = h(s—Ar)ds, 
A 


or 


ao 


fp) = h(w—r) de. 


A 


When 2 is positive this would be the representation of 
h(@—2) if the lower limit of the last integral were zero. 
This can be obtained if 4(w—X) is taken to be zero in the 
region 0<a#<vX, as represented in fig. 16. Again, if A is 
negative and 7) =—2’, the lower limit of the integral can 
again be made zero if h(#+ !) is taken to be zero in the 
region «<0, 1. é., if the function starts as indicated in fig. 1c. 

That the function e~”.f(p) is rightly interpreted above 
can also be proved simply by theorem (17), when e~* is 
interpreted as being unity for «>A and zero for a<n. 

Let us consider an instance. 

1 


h(a)=« is represented symbolically by or - =a (see 
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] 
fig. 2a); fig. 2b represents e~*?. 3 fig. 2a, 
however, e~”.— were developed as 
P 


Ap 


2 


Fig. la. 


ov 


Fig. 18. 


fl tp). 


[et f(p)]. 


and the positive powers of p were rejected, we should 
obtain 


—rA—-2-n, 


| 
x 
Fie le. 
| 
K 
| 
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Wan) 
(x) = p e > 
Fig. 2. Fig. 2 d. 
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A 
Fig. 2e. 
| 
0 x 
il 
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as represented in fig. 2c, which is erroneous. It is, however,. 
interesting to note that 


if n is a positive integer, and when the brackets [ ] are 
understood to mean the omission of all positive powers of p. 
Therefore the difference 


[[e~” f(p) -F(p)], 


(where the double brackets ne, the retention of the 
positive powers of p only, the zero’th power excluded) 
yields (see fig. 2 e) 


In the same way it es that 


Ller?.f(p)]]=0, 


for the factor e+” shifts the function to the left as explained 
above, there being no positive region at present where the 
function is continuously zero. 

These considerations are obviously only valid when h(2) 
can be developed as a series of positive powers of 2. Later 
on, in the treatment of the Legendre functions, we shall 
come back to these properties. The relation of these 
properties with a Laurent development has not yet been 
investigated. 

Returning to the meaning of e~”, it follows at once from 
the development of 


1 


il 
that 5 {1 +coth 


means the “ staircase ”-function depicted in fig. 3 *, because 
each following term adds discontinuously unity to the 
function. 

Again 


represents the meander-function of fig. 4 T. 


* See Heaviside, ‘ Electromagnetic Theory,’ ii. p. 90. 
1 Heaviside, doc. cit. ii. p. 97. 
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As an instance of (13) we take 


Applying (13), we find on the one hand 
dp = = tan7 = cot7!p. 
P p pti P|, 
On the other hand we have 


de = Si(x), 


Fig. 3. 


& 


and thus we find, as the peerstional expression for the 
integral sine, 


Applying theorem (14) to this result, we find at once 


sin 
= Lim. = 9° 


Again, let us consider 


Mp) = 


| 
x 
Fig. 4. 
0 
= cosa = h(x). 
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The application of (13 a) yields, on the one hand, 


=log Vp’ +1. 


On the other hand, by definition we have 


C08 da = —Ci(a), 


and thus we find as the operational expression for the 
integral cosine 


log——, = da (20) 


In the same way we can treat 
‘—z 
= 


For, with (13), (13a), and (14), we find 


and thus 


As another example of the application of (14), we take the 
Frullani integral * 


. 
e7 

0 & 


Writing this integral operationally, we obtain 
and thus we have found operationally 


Before taking up our main subjeet, the solution of linear 
differential equations, attention may first be drawn to the 


* Bromwich, ‘ Intinite Series,’ p. 488. 


| dx + — da+ = 
x Jc & Jo & 
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operational representation of logz. Here we first follow 
Heaviside’s example * 
Above we found from (1) 


a” 
Now differentiate both i with respect ton. We thus 
obtain 
where 


We therefore have the special cases 


loga—wW(0) = + C, 


1 


1 1 18 


(loge + C— 1—}—-4...--) 


With the aid of (189) we can now easily expand, 
e.g.,(20). Thus, 


1 1 
7) 
1 
08 + 358 


* Heaviside, Joc. cit. ii. p. 358. 

+ Here our notation for (mn) differs from that of Jahnke-Emde, 
‘ Funktionentafeln,’ but agrees with that of Schafheitlin, ‘Die Theorie 
der Besselschen Funktionen’ (Teubner, Berlin, 1908). As M(x) coincides 
with nz! for n=integer, we believe that it is advisable to use [-functions 
instead of P-functions. Thus, instead of defining 


d 
b(n) = (log 
we use the definition (18 /). 


lo 
(2). 
Z 2} 4! » ( 
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4, Linear Differential Equations : (a) with Constant 
Coefficients. 

Suppose the following simple differential equation be 

given :— 
Yytay = [cosat |, 

where the sign [ ] bracketing coswt means that the force 
as expressed by the right-hand member is zero for t<0 and 
is coswt for t>0. Let the initial condition be y(0)=0, 
from which it follows that y(0) =9(0)=7(0)=... = 0. 

To solve our equation operationally we must obtain the 
expression for 


f(p)=p\ e~Py(t)dt. . « 
0 

Hence we proceed as in our former paper*, multiply 
throughout with e?, and integrate between zero and 
infinity. Integration by parts of the first member yields 


to 


(p+a)X(p)= \ cos wt (22) 
where X(p)= | y(t) dt. 
« O 
As S(p) =p... 
we can write at once, from (21), 


The same result would have been obtained if we had simply 


replaced . by p and had written the symbolic expression 


for coswt for the second member. It may be remarked 
that the result (23) would also have been obtained if the 
original differential equation had been symbolically 


2 
or, written in ¢, 
yt w'y= 7, Le 
Hence the operational method enables us often easily to 
write down several differential equations which have the 


* Phil. Mag, vii. p. 1153 (1929). 
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same solution. It also answers the following question: let a 
right-hand member [1], representing the sudden application 
of a constant force, be given to the equation 


¥y¥+y=0, 


what must be done in order to have no transient phenomena ? 
It is easily seen that a second-order impulse, i. e., p*, must 
also be impressed on the system ; for if 


(p?+1).f(p)=1 +p’, 
then Ap=1, or y(t)=[1]. 
Let us take another example. The solution of 
y+y=0, with y=y(0) for t=0, 
is easily found to be 


Now let the equation be multiplied with (t—a), and the 
solution again be sought for the same initial condition. 
We have 


or 


ty—ayt+ty—ay=0. . . . . (24a) 
Again multiply by e~” and integrate. The first term thus 
becomes 


e 
dt dp dt 


d 


A similar reduction of the other terms of (24a) yields 


df 
ip y(0))+ + =0, 
or 
d 
(p+ + (ap+ )f=ay(0) p (24a) 
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The interpretation of the first term is simply 
yO) F—=y(0) 


but the second term must be explained with (9) as 


) 0 for t<a, or 
for. 


so that the complete solution is 


y(Oje* for t<a, or 
y(0)(14+ Ce*)e for 


meaning that in (24), through the multiplication by t—a 
the initial condition determines the solution up to the time 
t=a. At that moment, however, the multiplier t—a=0, 
and thus leaves the other factor in (24), y+y, free to take 
any value ; thus our operational solution shows that at that 
moment the amplitude may discontinuously be increased by 
the arbitrary value C, and atter that instant the effect of the 
initial condition is lost. If we had made a=0,.1.¢.,-i£ (24) 
had the multiplier ¢, instead of t—a, the initial conditions 
would have been lost from the beginning. As the example 
shows, we could have taken y(0) =0 from the start and solved 
(24a) simply with the second member equal to zero. Later 
ov, in (25), the integration constant C could have been 
adjusted to the initial conditions, and thus the right solution 
would have been obtained in a simple way. Obviously the 
multiplication of the diet differential equation by ¢ simply 


brings in a differentiation a throughout in the operational 
ap 2 
expression; similarly a factor ¢? causes dp? etc. Very often 
ap 
it is, therefore, much simpler to take 


y(O)=y(0) =9(0) =... =0 
from the beginning. 


5. Linear Differential. Equations :(b) with Variable 
Coefficients. 


Our, last example brought us already from the domain of 
equations with constant coefficients to those with variable 
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coefficients, and to the latter we will now turn our attention. 
Let us consider the equation for the Bessel functions, 


25 2 
Substitute 
and u'=4e; 
thus (26) becomes 
d*y d. 


Now multiply (27) again by e-?* and integrate between 
0 and o in order to find 


x={ e-P* .y(w) . dx, 
0 
S(p)=pX = y(2). 


As integrations by parts are mostly necessary, we can, 
beforehand, make a list as follows. Through the integration 


from which 


y becomes X, 


dX 
” dp 
dp* 
d 
” y(0) +pX 
dy d 
da dp? (pX), 
» —py(0) + p°X, 
d? d 
” y(0) dp (p?X), 
d? 


Phil. Mag. 8.7. Vol. 8. 
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Applying the transformations (28) to (27) we can now at 
once write down (taking 7(0) =0) 


or 
pr=f(p)=Cp™.e 


Adjusting © to the ordinary norinalization of the Bessel 
functions, we thus have found 


n 1 
which is valid for n>—1. Thus the development of 
the second member as a power series in p! gives us the 
familiar series for Bessel functions. 

Now that the operational representation has been obtained, 
we can, with the aid of our theorems (5) to (17), write down, 
almost at a glance, several properties of “pe 

(a) Differentiate the left-hand member of (29), 7. e., 
multiply the right-hand member by b. mh we find 


In( at) + In! (ae) = I(#), 


a familiar relation. 


(b) Again, multiplication of (29) by p” yields 


n n 1 


which is also a well-known relation. However, operationally 
it still has a meaning when n is not an integer. 
(c) Now apply (11) in the form 
py dh(a) 
dp 


which easily gives 
(d) Next apply (12) in the form 


1(p) 


dp 
F n——)pX 
)PX=0, 
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which yields, after reduction, 


(e) If we Poon apply (9), we obtain 
n 
where the development of the first member yields 


Each term represents a function of the form 
rn = 
m! 


and thus, after an obvious transformation, we obtain 


=(i-3)!  (a#>a) (34) 
(*) Next consider 


po 

and interpret the series as a series of Bessel functions. 

This series becomes 


n+2 


which, with 


Vl+a=y and 2Va=s 
and ee s again by a, becomes 


being the theorem for the multiplication of the argument of 
a Bessel function. 


2 


1 

“2 
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(g) Further, let us apply (14). This theorem gives 


1 o 
p at) dp= | dz, 


0 0 


or, after the substitution 2 /w=s, and replacing s again by 2, 


I (n) = (n>—1) (86) 


n+1 


(h) Next, with the definition (3), we can at once write 
down the integral 
1 


0 
which, after an obvious transformation, becomes 
1 
0 


n-1 2 
(j) As 
e? Jy(2 / —ix) =ber (2 2) +1 bei (2 Vz), 
and 


e P= / +ix)=ber (2 Wx) —i bei (2 V2), 
we obtain at once the operational expressions for Kelvin’s 
“ber” and “ bei ”’ functions as 


ber = 
(29 a) 
bei (2/27) = 


(k) Let us finally make use of the very fertile theorem 
(17) of the “ produit de composition.” Taking, e. g., 


gem 


A, (p) (n—m > —1) 


mel 

Ap)=p".e (2/2), (m>0) 

we find 
n 
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The application of (17) thus at once gives 


m—1 


or 


Qm-n 
I(n—m) (a?—£)"-™ e dE = x"J n(x) (38) 
(n—m—>1, m>0) 
Taking in (38), 
m=4, 
and making use of 


(38) becomes 


(cos cos (x sin w) do=V/ 7 . 
0 1 
(n> 


which Fae n=0 gives the familiar integral definition of the 
zeroth Bessel function 


0 
The results (30-39), which could easily be still further 
extended, were all obtained with the aid of the operational 
expression for 


J3,(20/2). 

In order to obtain another operational expression for a 
Bessel function we return to the differential equation (26). 
After the substitution z=u-"y, and writing x for u, we 
obtain 
dy 
= =0. 26 a 
(1—2n) ey (26 a) 
Now we transform (26a) again operationally by multiplying 
it with e~?* and integrating ; this is simplified by the use of 
table (28). Hence we can write at once 


where, as explained above, we took y(0)=0, and where 
again 


x 


X(p) ={" y(a) de. 


\ 
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The solution of (26 }) is easily found to be 
K(p)=pX(p) = 
and again, after adjustment of the integration constant to 
the usual normalization of the Bessel function, we find a new 
operational expression 
II(2n) 2nJI(n—4) p 


a" In(a), 
(40) 


which is valid for »>—4. As particular instances we note, 


for n=0 and n=}, 


P 
aw. . (42) 
wp +l: 


The last relation (42) can be interpreted by (186), giving 
the familiar result, 


Of course (40) can again he developed as a series in PO; 
thus yielding the expansion of J,(z) in a series of positive 
powers of Obviously_one in the former case, 
apply several of the theorems -(17) to (40) again. 
We will, however, only derive a ee ae with the aid 
of (17), in the following way. 

Take 


NP)= Gp 


then 


ACP) FAP) sine. 


Hence, applying (17), we obtain 


= sing, 
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which, after simplification and the use of 


becomes 
cos n (; ) (x—§&) (&).d&=sina, (44) 


which is valid for —}<n<+4. Moreover, the expression 
(40) is well adapted for a derivation (with the aid of (8)) of 
the asymptotic series for Bessel functions, but we shall not 
enter into these details. 


We now turn our attention to still another symbolic 
representation of Bessel functions, which will show us the 
close relation of these functions to spherical harmonics. 
Consider therefore again the fundamental differential 
equation (26). ‘This time, however, we write zx instead 
of z, and n+4 instead of n: 


2 1)2 
Substitute 


which gives 


2 
+20 _ (at 


We now obtain again, with the aid of table (28), the dif- 
ferential equation for ; 


(a) de 
as ide 


ax aX 
ees 1)\X=0 


which is the equation of the spherical harmonics. 


We thus find, choosing the right harmonic and the right 
normalization, 


from which we may conclude at once with (1) 


2 


0 n 


I (n).1(—n)= 
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where the I-function is as usually defined by: 

rao 

II(r) 

(45) shows the close relation between the Bessel function of 
imaginary argument and the Legendre functions, the latter 
being the operational representation of the former. In fact 
the familiar relations between Legendre functions on the 
one hand and those between Bessel functions (such as (30), 


(32), (83)) on the other are easily transformed into each 
other. E.g.: 


P (2) —P'@) = (2n+1) . 


is transformed, with the aid of theorem (12), into 
n—1 n+1 


and similarly for the others, 


Leaving the spherical harmonics for the present, we now 
turn our attention to still another operational representation 
of Bessel functions. We therefore obtain from (26) (taking 
u=a) directly through the usual process the differential 
equation for 


The result is ; 


of which the solution is 
C 


Choosing the proper solution and the right integration 
constant, we thus find 


(46) 


which expression we will study somewhat closer, as it gives 
us several properties of the Bessel functions in a very direct 
way. 


* Watson, ‘ Bessel Functions,’ p. 77. 


d =, dX 
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(a) Writing a= /1+p?—p, 
we have 
= J,(z). 
Further 
n 9 
n+1 1 2 4 
and hence 
r=0 n+1+2r 


(b) Again take the series 
P 
(c) Further, another series, 


(d) Again, as for integer n 


J (2) = J (2), 
—2n +2n 
we easily find the sum 
n=+oa 
> 2014+ 2324+ 204+... 


2a? + 2a* + 


as a matter of algebra. 
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(e) Next apply theorem (13) to (46). It follows that 


g 


The left-hand member equals 


and we thus find the meaning of a” as 


(n>0), Gt) 
0 
and further, 


pa" = p(V1+p?—p)"= (n>0). (51a) 


(f) Further, consider (46) in connexion with (17). 
Let us take as 


( 
Ap) = (x) 
Hence 


and thus we can write down at once the interesting result : 
=sina, (—l<n<-+1). (52). 
More generally we also have 


which, interpreted with the aid of (51a), gives us the relation: 


(n>—1, m>—1, n+m>0). (53) 


(g) Let us further inquire as to the meaning of A(z#),. 
where 


f\p)= log* = p = log (p+ Vp? #1) = Ne). 


] 
= (n>Q), 
nr | 
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To this end we start again with 
Jo(a). 
Next apply (13), giving us 
pt 
or sinh-!p= log p*+1)=—Ji(x), (54) 
where, analogous to 
*COS 
Ci(a) 


which function could be called integral Bessel function. 

This function, which does not seem to have been much 
studied, will appear again later on. 

(h) Finally, the application of (17) directly gives us some 
new and interesting relations, for we can write (46) as 


P n ) 
and thus we find, with the aid of (51a), at once 
m (x) (m>0,,.n>—1). (55a) 


A further interesting Frullani integral involving Jo(z), 
and closely related to our Ji function, can also be easily 
evaluated by the operator method. Let it be required to find 


la=Y. 


Jo(x) — cos 


Here the application of our theorem (14) gives us at once 
the required result, for 


ao 


=log 2, 


and hence 


\0 


2/0 
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6. The Derivation of J,(«) with respect to n. 


In § 5 we obtained, among others, three operational repre- 
sentations of Bessel functions: 


n 1 
(n>—1). . . 


(n> —1). . (46) 


BI, 
On 


entiate all three expressions with respect to n. We thus 


P 


In order to obtain new expressions for , we differ- 


obtain, respectively, 


(56) 
an { log + 
2"II(n—3) P 2 
V0 (p?+ 1)*+i {log 2 +(n—4)—log (p +1)} 

(57) 

OIn(#) P —nsinh—1 
Tors — sinh 


Now above we found already the three functions in (56), 
(57), and (58), namely, 


—logp=loga+C . . . (189) 


— log = Cia), 

— sinh" p= 

and therefore we are prepared to write the right-hand 

members of (56), (57), and (58) with theorem “ay as 

‘produits de composition.” The three equations thus yield, 

after some slight transformations, the following interesting 

relations : 

(n>—-1) (59) 


2° I(n— =) (n> 1), (40) 
Ln nl ZX) = (p 1 
J,( 2) = V1 
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_ 4 ¥( n—3)— logs Jn( x) 


Ji(e—£) . J,'(£) dé, (n>0). . . (61) 


All three relations are, as far as I know, new, and the 
ease with which they can be derived with the aid of the 
operators is remarkable. 


7. Legendre, Laguerre, and Hermite Polynomials. 


The familiar differential equation for the Legendre poly- 
nomials is 
d*y dy 
—x’) +n(n+1)y=0. 
Substitution of 
a=s+1 
vields 


(st +28) + 2(s +1) —n(n-+1)y=0. (62) 


Applying the operational transformation on (62) gives us,. 
with the aid of (28), 
n(n +1) 


dp ap 
where Kp) =y(s) 


or, with the proper normalization, and for integral values of 
n, we obtain 


p dp?) 

Hence the operational representation of a Legendre poly- 
nomial is given by a Bessel function, in the same way as: 
we found before that the operational representation of a 
Bessel function was given by a Legendre polynomial (see 
(45)). Obviously this symmetry is no accident and we will 
return to it later on. It is interesting to consider (63) in 


f=0, 
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relation to what has been found before in § 3 in connexion 
with the operational meaning of e~?. Hence, if we 
multiply (63) by e-?, we obtain 


(63 a) 
where, as before, by the brackets [ ] is meant the omission 
of all positive powers of p (the zero’th power excluded). 
With the aid of (3) we therefore obtain the interesting 
integral 
p) [oH (py), (64) 
—(n+3) 
where in the result the negative powers of p (zero’th power 
included) are to be retained only. 
From (9) and (63) one can also find 


da =p (=): 


We now return to another consideration of the Legendre 
functions. If we change the independent variable of 8, 
where 

x= cos 
we obtain for the differential equation for P, (cos 8), which 
for short we write P,,(@), 


dP 
A further substitution of 
yields 
) 


In this case, where the ccefficients contain exponentials, 
the usual transformation leads to a difference equation 
instead of a differential equation. In fact, we obtain 


(p + mip — n—1)X,(p — 1) 
(66) 


~where, as before, 
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From (66) the following solution follows at once 
—1 
and when we now return to the variable 6 again, and there- 
fore let p relate to 0, instead of to A, rer becomes 


= (68) 


24 


PXn(p)=Yn(p); 
Y,.(p) = P.(@), 


Now, writing 
where, therefore, 


(68) becomes 
pn 
n° 


which expression can be regarded as a “ produit de com- 
position.” Hence the application of (17) gives us immedi- 
ately the interesting and apparently new relation 


sin 


The directness with which this result is obtained with the 
use of the operator is again striking. 

Returning to (69), and taking the integration constant of 
the difference equation (66) or (69), such that 


Po(@) 1, 


we obtain as the operational expressions for the Legendre 
functions 


P,(6) =1, ) 

PO 


+ + p? + #) 
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Letting finally, in (69), n> , this expression becomes 


Lim, Y,*(p) = 


p 24 42? 
or Lim. = 
= 
P 


which second member, as we found before (see (41)), repre- 
sents Jo(n@), and hence we found 


Lim. P, (cos) = 


n> 


Finally, we consider the differential equation for the 
Laguerre polynomials L, : 


It transforms into 
d n 
+ 
where Kp) 


Hence, with the proper normalization; we find 
=a! (1- (n=integer). (74) 


Similarly we obtain from its differential equation, as the 
operational expression for the Hermite’s polynomials He, 


—p? 
(uainteger), (75 
e (5) (n=integer), (75) 


where [ } means the retention of negative powers of p only. 

It is, ee a quite simple matter to obtain from (74) 
and (75) the properties of the Laguerre’s and Hermite’s 
polynomials. 


+nlin=0. . 


8. Final Considerations. 


From the many examples discussed above of the opera- 
tional transformation of a differential equation for y(#x) into 
another differential equation for f(p), where f(p) = y(«), and 
alsofrom a general consideration of the method, it follows 
that there exists a close analogy between the well-known 
Laplacian transformation and the transformation discussed 
above. 
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In a Laplacian transformation we express as an expo- 
nential integral : 


while in our case we calculated 
2/0 


As was pointed out in connexion with (3) and (4), with 
a suitable choice of the limits a and 0 in (76), or with a 
suitable choice of the path of integration in the complex 


plane, 

or, in other words, the solution of a Laplace integral, con- 
sidered as an integral equation, is again a Laplace integral. 
Hence many of the operational expressions f(p) for func- 
tions y(#) obtained above have exactly the same form as the 
integrand of the Laplace integral expression for y(#). In 
fact (76), (77), and (78) can probably be considered to 
constitute a Fourier integral identity. 

A very interesting list of Fourier integrals, published 
recently by G. A. Campbell * can therefore be considered 
also to be a list of operational representations of the functions 
examined. 

The ultimate identity between the operational solution of 
a differential equation and the solution in the form of a 
Laplace integral clearly shows when the operational treat- 
ment of a problem has advantages over the normal analytical 
treatment. When, e. y., functions such as Bessel functions 
oceur, an operational treatment clearly is much simpler and 
very often much more direct than the normal way. 

The above examples, which could easily be considerably 
extended, where many properties, known and new ones, 
e. g. of Bessel functions, could be derived in a most simple 
and direct way, were, in fact, intended to prove this statement. 
If, on the other hand, a physical problem would lead, say, 
to the function y(#), where 


na) 
the operational method would lead to a Bessel function 
because this last expression is itself the operational repre- 
sentation of a Bessel function. In sucha case the opera- 
tional method would have no advantages. Hence the 
* The Bell System Technical Journal, vii. p. 639 (1928). 
Phil. Mag.8.7. Vol. 8. No. 53. Suppl. Dec. 1929. 30 
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operational way of solving physical or technical problems 
shows its full advantages in those problems where complicated 
functions occur which cannot be expressed directly by more 
elementary functions. It often happens that the operational 
representation of those complicated functions is given by 
rather elementary functions, as was demonstrated in many 
of the above examples, and in those cases the directness and 
simplicity with which the properties of those complicated 
functions can be studied is most remarkable. 

Finally we wish to remark that several formule derived 
above are also valid outside the range mentioned. In fact, 
negative powers of a (except negative integral powers) can 
be represented operationally by positive powers of p (except 
positive integral powers). However, under these cireum- 
stances difficulty arises from the integral definition in $3 of 
the II functions. In this ease the II(n) functions have to 
be extended into the domain n< —1 (e. g., by means of the 
Gauss limit of an infinite product). 


In conclusion we wish to express our great indebtness in 
preparing this paper to the writings of Heaviside, especially 
volume ii. of his ‘ Electromagnetic Theory,’ and to the late 
Dr. Bromwich of Cambridge and Dr. Le Corbeiller of Paris, 
and Professor Elias of Delft University, with whom I had 
the opportunity to discuss several of the questions treated 
in this paper. My thanks are also due to Dr. K. F’. Niessen 
(of this Laboratory) for verifying most of the formule. 


List of some Operational Representations derived above 
and of some other Results obtained. 


pl1] = 8 (2) (Heaviside’s “impulsive” 
function p[1]=Dirac’s 


coth = “staircase ” function of fig. 3. 
il: hp ” 
5 tanh = “ meander ” function of fig. 4. 


cot-!p = Si(#) (integral sine). . . (19) 
—log ¥1+p? = Ci(x) (integral cosine). . (20) 
— log (p—1) = Hi@). 


il il a 
(189 


= (log . (18 f) 
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1 n pa 
(n>—1) (29) 


cos (2/2) Kelvin’s “ber” 


and “bei ” ft 


bei (2 /z) tions. . . (29a) 


Qn 
4). Gap Iz). 


e—nsinh—! p = Jn (a) 


/1+p? 


inc- 


(40) 
(45) 


(46) 


Jo(x) 


— sinh 1p = Ji(x)= dx 


(r) | = P,(#) (n =integer). . . (63a) 


n+3) 


nt(1—=) = L (a) (Laguerre’s polynomials). 


= He, (5 ) Hermite’ s polynomials). 


n 


(2>a). 


Il(n—m) ,'0 
302 


(integral Bes- 
z function). 


(54) 


(74) 


(75) 


(37) 


(38) 


r=0 

(34) 
(35) 
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COS nar 1) £).J() dé = sina 
(—i<n<+4). (44) 


dé = sina 
(—l<n< 41), 


m z dE =I (2) (m>0). (65a) 
nm +m 


da = log2.. 


0) + £).(2) 3 @ag. 


(n>—1). (69) 


+2| (dk. (60) 
! 


(64) 


(n=integer.) 


—(n 


sin nO 


(70) 
(n=integer.) 


Lim, 


(cos). (cos (@— p)) ¢ dp= 


0n—1 


Summary. 


§1 considers the earlier work on operators by Riemann, 
Heaviside, Pincherle, Carson, Lévy, and Bromwich. Brow 
wich’s complex integral gives the solution of Carson’s 
integral considered as an integral equation. 

fin § 2 some veneral theoterns 1 in connexion with the opera- 
tional method, and formulated by Carson, are considered, and 
some new are added. “impulsive ”” 
function p|1] is shown to be identical with Dirac’s 8 (#) 


‘ v ) 
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function. The very fertile theorem (17) is considered as a 
“produit de composition’ as defined by Volterra, and used 
by Borel in his ‘“‘sommation exponentielle ” of divergent series. 
It was known in physics as the Hopkinson—Boltzmann 
superposition principle for linear systems. 

In §3 some of the general operational theorems are 
considered in more detail, especially the meaning of e~?. f(p) 
by itself, and further, the meaning of the same expression 
after omission in its series development of all positive powers 
of p. The results are illustrated by figs. 1 and 2. Further, 
the “staircase” function of fig. 3 and the “ meander ” function 
of fig. 4 are found to be expressible in a very simple 
way with the aid of operators. Further, with the aid of the 
above theorems symbolic expressions are derived for Si(«), 
Ci(z), Hi(z), and a” log x. 

§ 4 considers some examples of the solution of differential 
equations with constant coefficients, and a simple operational 
method is considered for finding the conditions for the 
elimination of transients in switching operations. It is 
further shown how the factor (t—«a) in (t~—a)(y¥+y)=0 
fully destroys the effect of the initial condition after the 
instant t=a, thus leading up to 

§5, where linear differential equations with variable 
coefficients are discussed. Starting from the differential 
equation for Bessel functions, the operational representation 


n 


is derived. Several of the familiar relations existing 
between Bessel functions can at once be derived from (29). 
Similarly, the multiplication theorem of the argument of a 
Bessel function (35) can easily be found. Further, the 
operational relations 


ber (2 / x) = and bei(2 = (29 a), 


are obtained, while the application of theorem (17) yields in 
a straightforward way the solution of some fairly complicated 
integrals involving Bessel functions. Next, another opera- 
tional expression for a"J,(#) (40), is derived from the 
differential equation, which enables us to obtain several 
other interesting integrals containing a Bessel function. 
Further, a new operational representation, 
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involving spherical harmonics is obtained, from which the 
analogy of the relations existing between the spherical 
harmonics are shown to be at once derivable from similar 
relations existing between Bessel functions. A further 
operational representation, 


V1+p? —p)"= J,(@x) 
is found to lead easily to the derivation of several series and 
integrals involving Bessel functions. 
Special attention is given to the function 


Ji(x)= dx, 


which is formed in a similar way to Ci(#), and which plays 


an important part in an integral expressing Oot), Next 
some Frullani integrals are shown to be easily solved with 
the aid of operators. 

§ 6 discusses three new integral expressions for the deri- 
vation of a Bessel function J,(«) with respect to its order n. 
These three integrals contain log, Ci, and Ji functions 
respectively. 

In $7 operational expressions for Legendre, Laguerre, 
and Hermite’s polynomials are obtained. Legendre poly- 
nomials, considered as a function of @, lead to a difference 
equation for its operational representative, from which an 
interesting and apparently new integral relation (70) 
containing Legendre polynomials can at once be derived. 
Similarly, the asymptotic behaviour of P, asa Bessel function 
when n> is easily demonstrated. 

§ 8 contains some general considerations on operators. It 
is shown that the operational representation of a function is 
closely related to its representation as a Laplace integral. 
It is further pointed out when in physical problems the 
operational solution is to be preferred to the ordinary solution. 
In those cases, which are many indeed, the operational 
solution leads to the desired result in a most remarkable 
and direct way. 

At the end of the paper a number of operational expressions 
derived are summarized, together with some of the integral 
relations obtained. 


der 
N, V. Philips’s Gloeilampenfabrieken, 
Eindhoven, July 1929. 
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XCVI. On the Damping of a Pendulum by Viscous Media. 
By F. E. Hoars, B.Se., A.R.C.S.* 


Introduction. 


FP XHE following experiments were undertaken in an attempt 

to develop an alternative method for the measurement 
of viscosity. The torsional oscillation developed by O. Meyer, 
in which a horizontal disk, supported by a suspension wire, 
performs oscillations in its own plane, is well known, and 
although the theory is not above suspicion, the usual formula 
gives results for the viscosities of liquids which are in fair 
agreement with those obtained by other methods. The 
theory of this method was first developed by Stokes +, and 
in the same paper he gives the theory of a simple pendulum 
performing oscillations in a viscous medium. 

The equation of motion of a pendulum consisting ot fa 
sphere attached to a fine wire and performing oscillations in 
a vertical plane, the sphere being immersed in a viscous 
medium of density p and coefficient of viscosity 7, can be 
written according to the classical theory f, 


dx g 
de —(M—M');.<. 
where 


x = displacement at time ¢, 


= radius sphere, 
T = periodic time, 
M = mass of sphere, 
M’ = mass of fluid displaced, 
1 = length of pendulum, 
g = acceleration due to gravity. 
* Communicated by Prof. F. H. Newman, D.Sc. 


+ Stokes, Camb. Phil. Soe. ix. p. 8. 
Lamb, ‘ Hydrodynamics,’ p. 635, 


2a 
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This equation is formulated without taking into account 
either the resistance experienced by the wire, which will in 
general be small compared with that experienced by the 
sphere, or the finite are of swing of the pendulum. Thus 
the equation will apply accurately only to infinitesimally 
small oscillations. 

The solution of equation (1) may be written 


M'9rk! 


where X and vy are constants dependent upon the initial 
conditions, and 


2 


The second term under the root sign in the expression 
for @ is small and can be neglected in comparison with the 
first term. Hence the time of vibration is given by 


The periodic time, T,, in vacuo, which is not greatly 
different from that in air, is 


=> 


M—M’° 


By substituting for k and rearranging the terms, it will 
be seen that 


and hence 


If T and T, can be measured with sufficient accuracy this 

affords one method of finding y. The results obtained in 
this way will be given later. 


M-kM' 
M—M’ 
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The logarithmic decrement of the amplitude is, from (2), 
8(M+ kM’) . . . (4) 
The quantities & and k’ depend upon the periodic time T, 
but as this is independent of the amplitude according to the 


theory given above, should be constant and its measurement 
would thus provide another method for ak the viscosity 


If, when substituting for k' in (4), Ha >; is neglected in 
1 
comparison with Br? and 2X in comparison with 7, the 
expression for X reduces to the simple form, 
9aM’ 1 
CD) 
4(2M+\’)" Br 


The relative magnitude of the terms as and * can be 
readily estimated by assuming for water the following values 
as rough approximations : 


=°01 ¢c.g.s., 


Br 


With these values, the approximate values of ze and 


are 


Birt = 400" 
Consequently, no serious error for the present purposes will 


be made by neglecting From the results of the 


l 
-experiments performed, the error produced by neglecting 2A 
in comparison with w is about 2 per cent., which is equally 
negligible. 


EKaperimental. 


Originally it was proposed to determine the coefficient of 
viscosity by finding the period of a pendulum in air and then 


p=, 
Pr 2 

34sec. 
iI il 
Br 20° 
1 1 
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immersed in water, using equation (3). For this purpose 
three brass spheres of diameter 14 inches, 12 inches, and 
2 inches respectively—each accurate to two-thousandths of 
an inch—were obtained, and a fine german-silver wire 
of -027 cm. diameter was fastened to each by drilling a small 
hole, heating the ball, and then filling the hole with soft 
solder and immersing a few millimetres of the carefully 
cleaned wire into the molten solder. It was found that the 
wires so fastened held very well and the method employed 
was not at all likely to disturb the sphericity of the balls. 
A large tank about 3 ft. 6 in, by 1 ft. by 1 ft. was used to 
contain the water in which the bails oscillated. Each ball was 
suspended so that when vibrating it was at a mean depth of 
about 5 inches below the surface of the water. Itis necessary 
to employ a large tank, as in the development of the theory it 
is assumed that the fluid is of infinite extent and at rest at 
infinity, and, if a small tank is used wall effects would vitiate 
the results. By using a tank with the above dimensions it 
was considered that no correction would be necessary from 
this cause. Actually, of course, in calculating 7 a correction 
should be made for the resistance experienced by the wire, 
but as the diameter of the wire was so small in comparison 
with that of the spheres employed, and only a short length 
was immersed, it was considered negligibly small, and 
consequently ignored in the calculations. 

The chief experimental difficulty in measuring the periodic 
time of the pendulum when the ball is immersed in water is 
caused by the rapidity of decay of the oscillations, it being 
impossible to obtain more than about 70 vibrations. Pre- 
liminary experiments using a stop-watch to time the vibrations 
showed that much greater accuracy in measuring the time 

would be necessary before this method could be employed 
for finding 7. Aceccr dingly it was decided to use a chrono- 
graph capable of measuring 1/100 of a second, the instrument 
being controlled by a tuning-fork and the order of accuracy 
1/100 sec. The results of preliminary experiments with the 
1} inch sphere are given below. 


Swings 
OR 60 9:9153 008 
100 26595 
65 . 2°8327 


100 ide 2-5794 
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These results were considered too discordant to make the 
continuance .of these particular experiments of any value, 
either as a test of the theory or as a method for finding the 
coefficient of viscosity, and it was decided to find 7 by 
the measurement of the log. dee. 


Logarithmic Decrement Kxperiments. 


In these experiments the log. dec. was measured for each 
sphere for two different lengths of wire suspending the ball. 
Errors arise in the log. dec. if forced oscillations are set up 
in the support, so a means of suspension which was at once 
rigid and easily adaptable was sought. The method finally 
adopted was to pass the wire to which the brass ball was 
attached through a fine drill chuck which could be screwed 
up tight, so that all four jaws were holding the wire. This 
drill chuck could be fastened, for the longer periods, into a 
bracket made of half-inch square brass screwed to a wooden 
beam in the ceiling, and for the shorter periods it was fitted 
into a bracket mounted on a solid cylinder of iron 2 in, in 
diameter which stood on a massive iron foot on a table. 
In neither case was it at all probable that appreciable forced 
oscillations would be set up. 

The amplitudes were observed on one side only with a 
telescope containing a scale in the eyepiece situated 24 metres 
from the suspending wire. At this distance the image of 
the wire when properly in focus was fairly fine, and little 
difficulty was experienced in estimating to a tenth of a 
division in the eyepiece. The greatest angle through which 
the wire was displaced from the vertical in any experiment 
was less than 1°8°, and in every experiment the pendulum 
made several complete vibrations before any readings were 
taken. 

According to theory the amplitudes should be in decreasing 
geometrical progression, and this was tested before proceeding 
further. The ratios of successive cbserved amplitudes were 
not sufficiently accurate for this law to be tested owing to 
the error of observation in each measurement, and it was 


therefore decided to find log —° for a series of observations. 


In the following table the amplitudes are given in eyepiece 
scale divisions and the last column gives log a , the 
n+1) 


results being for the 2 in. sphere vibrating with a period of 


3°58 sec. 
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50°0 26°5 6349 
25°4 5982 
43°5 24°1 5906 
40°4 22'9 ‘D677 
38:2 5564 
21:0 5418 
384°1 5336 
32°3 1971 “D244 
30°9 18°2 "5293 
29°2 5235 
16°6 


It will be seén that, except for one value of the logarithm, 
there is a progressive decrease in the third column showing 
that the assumption of successive amplitudes being in geo- 
metrical progression, and thus that the resistance experiensed 
by the sphere is directly proportional to the velocity, is 
unjustifiable. 


Modified Theory. 


To account for the variation of the damping as found by 
experiment it was assumed that the decay of the amplitude 
could be resolved into two parts, one proportional to the 
velocity and the other to the square of the velocity. The 
differential equation for the variation of the amplitude can 
then be written 

da 
where @ is the amplitude at time ¢, A and B being constants. 


Integrating this equation, and substituting the end con- 
ditions, 


| 


== 0) 


a=a,att=h, 
we have 
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This solution involving constants in the logarithm is 
inconvenient for calculation, so the following method due to 
Stokes * was adopted. 

Rewriting the equation (6) in the form 


—Adt —Bzdt, 
a 


we can substitute for « in the last term from the equation 


an approximation in the present case, but accurately true 
when the amplitudes decrease in geometrical progression. 
Substituting, and integrating, we have 


a 
ao 


log — 


This equation involving two disposable constants will 
represent the experimental results as was shown by one or 
two trials, and from it an equation for 7 can easily be 
developed by the use of subsidiary hypotheses in the following 
manner. 

The dimensions of the right-hand side must be zero, so 


ing log “° to depend only on tl iti 
assuming log a to depend only on the four quantities, 7, T, 
r, p, Where these svmbols have the same meaning as before, 


it follows from a consideration of dimensions that log oe 
must be of the form 


a 
log = K 
ay 
where s is some as yet undetermined power. If the theory 
already given be assumed to hold as a first approximation, 
we find from equation (5) that 


$= 4, 


Consequently, equation (7) becomes 


* Stokes, loc, cit. 


&, \ 
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where nT = 
n = complete vibrations between a) and a. 
The quantities «) and a, in equation (8) are the angular 


amplitudes. In the term log— the ratio will be the same if 


1 
linear amplitudes or angular amplitudes are used, but for the 
term involving a#—a the angular amplitudes must be 
obtained. This is easily done by setting up a metre scale in 
the same place as that in which the pendulum oscillates and 
reading off the number of centimetres corresponding to 100 
scale divisions in the eyepiece. The mean result of several 
observations showed that 100 scale divisions corresponded 
with 5°2 cm. ; 
div. = em. deflexion of wire. 


Also the angular amplitude corresponding to one scale division 
‘052 


> 


where a is the distance in cm. from the point of support to 
the point on wire where the amplitude is observed. 
“Hence, if “a” is the amplitude of an oscillation measured 
in cm., equation (8) can be written 


log? = Ka / 4 (9) 
log— 
ay 


The method of testing this equation was to substitute in 
it the values obtained for dg, ao, and dg from a single 
experiment, using the 2 in. sphere vibrating with a period of 
3°58 sec. and assuming the value *0110 c.g.s. for the viscosity 
at 17°C., at which temperature the experiments were per- 
formed, and 1 for the density of water. Values were thus 
obtained for K, and B,. The latter was next assumed to be 
some function of the periodic time, and the form of this 
function was obtained by considering the experiments with 
the same sphere, but vibrating with a shorter period. In 
this way it was found that B, could be written 


IBS = B; T?, 

Finally, a further assumption was made, viz., that B, was 
some function of the radius of the sphere. From the experi- 
ments on the 1? in. sphere vibrating with a period of 
3°568 sec. it was found that 


r 
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Rewriting equation (9), we now have 


ao nl a—a 

og 

ay 
To verify all the assumptions the value of the viscosity » 
was calculated from the results obtained with the three spheres, 
each vibrating with a long and a short period. For this 
purpose it does not matter if “a” is given in cm. or in scale 
readings, Ky alone being affected by a constant factor. In 
the following table “a” is therefore given in scale deflexions, 
but the constants utilized in determining the results are in 
those appropriate units, so that ‘‘a” is in em. to facilitate 

comparison with future work. 


Results. 


The mean value of w was taken to the nearest cm. for the 
long and short pendulums. 


= 84 cm. 

216 cm, 

Values used for calculating constants in equation (10) : 
= 0110 c.g.s, 


For short pendulums 


& 


For long pendulums......... x 


p=l1. 

BOC. 

r = 2°54. cm. 

a = 50°0 divisions. 
Ay = 27°7 a 


— 
1 seale division = ‘052 cm. deflexion of wire. 
These values, when substituted, give the following values 
for the constants : 
Ky, = 3786. 
K, = 10°04. 
Hence, if ‘‘a” is measured in cm., 


m2 
a 


ay pl? 


The results are divided into Section A and Seetion B, the 
former relating to the 2 in. and 1? in. spheres which were 
used indirectly in finding the form of equation (10), and the 


ay 
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latter referring to these experiments on the 1? in. and 13 in. 
spheres which have not been used at all in the development 


of the equation for ». 


SECTION A. 


Sphere. T secs. ne.g.s. Means. 
2 in. 3:58 47-0 26:9 ‘0097 
49-0 27-1 ‘0116 $ 0105 
49°5 27°8 0101 
2-48 48°5 27:1 
48°2 26-9 -0099 | 
489 27-2 “0098 
43'8 26°3 ‘0101 | 
47-9 26:8 ‘0098 
49:0 27:1 ‘0103 / 
12 in. 3°568 49:0 0098 ) 
25:0 0098 | 
48-4 25°2 0097 “0101 
46'3 24°] 0106 | 
Section B. 
Sphere. T secs. dp: ao ne.g.s. Means. 
12 in. 2-484 44:0 22°9 0097 \ 
46:0 23:2 ‘0125 
49-0 0127 | 
48:2 24:3 | 
459 23°6 0102 / 
1d in. 3:5538 46'8 21:0 *0125 
49-2 22+] ‘0114 
48°7 21°9 0115 
48°7 219 0115) 
2°52 47-0 20°9 ‘0115 
21-0 0105 | 
46:2 20:5 0110 $ 
460 20°6 0105 | 
47°1 20:9 


| 

| 
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Conclusion. 


The results in Section B show a close approximation to 
the true result, viz., ‘0110 c.g.s. units, and indicate that the 
theory accounts for the observed phenomena. The mean 
value of from the results in this section is ‘0112 ¢c.g.s. units, 
which is in remarkably close agreement with the true result 
when it is remembered that the accuracy of the constants 
K, and K, depends entirely upon the accuraey with which 
merely three readings of the amplitude have been made in 
the case of the 2 in. sphere. 

It is suggested that the physical meaning of K, and K, 
may be found from further experiments with different liquids, 
and these are about to be undertaken. 


XCVITI. On Conditions near the Cathode of a Glow Discharge. 
By Nora M. M.Sc., and K. G. EmMEtius, 
M.A.. Ph.D., Department of Physics, The Queen’s 
University of Belfast *. 


’ hs cathode dark space is not a simple unit of the 

glow-discharge, but includes, when fully developed, 
a number of partially distinct sections, there being encoun- 
tered in turn, on proceeding out from the cathode, the 
primary dark space, the cathode glow, the cathode dark 
space proper, and the transition region between the dark 
space and negative glow. The present communicationis a 
preliminary account of an investigation of the region of 
the cathode glow and primary dark space with an exploring 
electrode, for conditions not far removed from those of the 
normal cathode fall in potential f. 


I. Experimental procedure. 


The form of tube used is shown in fig. 1. The cathode 
C consisted of a hollow, nickel-ended iron cylinder, with an 
iron tail T resting in the end-tube E; with this construc- 


* Communicated by the Authors. 

+ The work described was done early this year. Its completion has 
been unavoidably delayed, and although we hope to proceed again with 
it shortly, it is believed that the results already obtained are of sufficient 
interest to be published at this stage. 


Phil. Mag.8. 7. Vol. 8. No. 53. Suppl. Dec.1929. 3P 
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tion small changes were readily made in the position of the 
cathode with a magnet. The anode A was a nickel disk. 
The exploring electrode was a filament of tungsten, 
0:09 mm. in diameter, or molybdenum, 0-107 mm. in 
diameter, strung across the tube in front of the cathode, 


Fig. 1. 


@cms. 


and kept taut by,a 10 gm. piece of copper tube G. Care was 
taken that the glass sheath of F did not touch F for the last 
6 mm. of its length on each side of the discharge, and that 
the filament was parallel to the surface of the cathode, 
although the latter precaution was found not to be of 
great significance. The metal parts were degassed in a 
vacuum furnace before assembly, and before admission 


| | 
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of any specimen of gas, the tube was baked out with a 
hand-flame at a pressure of less than 10°-*mm. Hg. That 
part of the surface of the cathode which wasimmediately 
opposite to the filament was also repeatedly subjected to 
electron bombardment from the latter. The gases used 
were oxygen, nitrogen, and neon, prepared as has been 
described elsewhere *, hydrogen, admitted through a 
pailadium regulator, and technical argon, containing some 
20 per cent. of nitrogen. Pressures were recorded by a 
McLeod gauge, and condensible vapours were kept from 
the tube by liquid-air traps. 

Simultaneous measurements were made of the current 
to F (z,); the current through the tube (2,) as measured 
by the current to the anode; and the difference in poten- 
tial between collector and cathode (V), and cathode and 
anode. The discharge was always allowed to pass for at 
least an hour before any readings were taken, and was then 
so steady that it was not necessary to correct 7, for 
fluctuations in %. The potential applied between anode 
and cathode ranged from 250 volts to 450 volts, and the 
gas-pressures from 0-2 to 1-5 mm. Hg. The filament was 
always between 0-3 and 0-9 mm. from the cathode, and no 
change in its position was noticeable as its potential was 
varied. The results were not affected in their main 
features by change of the discharge conditions between 
these limits. 

A typical collector characteristic is shown in fig. 2, in 
which a current below the axis of voltage is one in which 
there is a net flow of positive electricity from the gas to 
the filament. The zero of potentials is that of the cathode. 
The appearance of the curve is similar to that of the curve 
obtained when a collector in the negative glow is starting 
to act as a main cathode for the discharge f, but there are 
important differences in detail, the ‘curve of fig. 2 
exhibiting more or less sharp discontinuities of slope at 
the points A, B, and C (B is usually sharper than in this 
particular instance). 

We are not sure that the discharge was not seriously 
disturbed locally when the filament was hot. Very small 
changes in % also probably correspond to the breaks 
A, B,and C. Discussion of these two points is reserved for 


* Emeléus and Beck, Phil. Mag. vili. p. 121 (1929). * 
+ Emeléus and Brown, Phil. Mag. vii. p. 17 (1929). 
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a later paper. All data given in the present communi- 
cation were obtained with F cold. 


II. Interpretation of the Characteristic Curves. 


It is of fundamental importance to know if an ionic 
sheath forms round a collector in the cathode dark space. 
Theoretically, the problem is indeterminate, since in the 
usual integration of Poisson’s equation it is necessary to 
know the field either at the surface of the collector or at 
some point in the sheath, supposed present, or at its 


Fig. 2. 
le 
MICROAMPS 
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1h3 = 


22:7 
28:3 
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V : POTENTIAL OF F RELATIVE TO CATHODE : VOLTS 
Tube voltage=310. F to cathode=0°6 mm. 
Gas-nitrogen. Thickness of cathode dark space=7 mms. 


outer boundary; to assume a value for the field in the 
present case is tantamount to assuming the existence of 
the sheath. The problem can be approached directly 
in two ways from the experimental side; first, by visual 
observation, and, secondly, from the shape of the 
characteristic curves. Visual observation is impossible 
with the tubes that we have used, but has recently been 
made the subject of a special study by W. L. Brown and 
E. E. Thomson in this laboratory *. Their observations, 


* See the following paper, p. 918 
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which are reported upon elsewhere, show that the well- 
known shadows thrown by obstacles in the cathode dark 
space can be used to find space-potentials, but are not 
conclusive as to the existence of sheaths like those formed 
in plasma-type discharges. We have thus to rely upon 
the form of the characteristic curves, and this appears to be 
irreconcilable with the existence of ordinary sheaths. 
There are several reasons for arriving at this conclusion. 
First, simultaneous observation of shadows and of collector 
currents in Brown and Thomson’s apparatus, in the part 
of the cathode dark space in which we are now interested, 
has shown that the space-potential lies at or close to the 
break B. There is in this part of the characteristic no 
sign of the reception of the groups of slow electrons which 
produce a rise in the characteristic curve in other parts of 
the discharge at the space-potential, and, in spite of the 
large electric field in most of the cathode dark space, some 
of these slow electrons might be expected to be present 
near the cathode glow. Secondly, the position of the 
characteristic curve relative to the axis of voltage shows 
that, when insulated, the collector charges up positively 
to the surrounding space, as Aston suggested in his 
explanation of the reputed ‘ Kathodensprung” in 
potential *. This shows that there are present large 
numbers of positive ions, or that the collector is emitting 
electrons, or that both of these effects, the latter of which 
certainly tends to destroy a sheath, are co-existent. 
Positive ions are, of course, known to be present in excess 
of electrons through most of the dark space, although the 
electrons may preponderate near the cathode. From 
such considerations we have come to regard the action of 
a collector in the dark space as more nearly like that of a 
grid wire in an evacuated thermionic triode, that is, its 
effect is not confined to a region bounded by a separate 
ionic sheath, but extends without existence of a definite 
boundary into the pre-existing space-charges of the 
cathode dark space, modified by the ionic shadows thrown 
by the collector. 

Conditions appear to be too complicated to make a 
mathematical formulation of these ideas possible, but we 
can account for the course of the characteristic curves 


* Aston, Proc. Roy. Soe. lxxxiv. p. 526 (1910). 
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qualitatively on lines indicated by Aston*. Fig. 3a, 
which is taken from his paper, shows the probable trajec- 
tories of positive ions (thick lines) and electrons (thin 
lines) near a—positive—insulated collector F. We are 
concerned with what happens when F is made more 
negative than in this case, and it is evident that this will 
be accompanied by a gradual straightening of the lines 
until when F is at the space-potential, they run parallel 
to each other and perpendicular to the cathode. At this 
stage the space between the lines tangential to F will be 


Fig. 3, 


shielded from the streams of ions and electrons and will 
appear as a shadow of these dimensions, this constituting 
Brown and Thomson’s method for finding the space- 
potential. Meanwhile, only small and essentially steady 
changes will have been occurring in the current to the 
collector, so far as these ions and electrons are concerned. 
After passing the space-potential F will affect the positive 
ions which have passed it on their way to the cathode, so 
that their paths are more like the lines in fig. 3c. The 


* Aston, Proc. Roy. Soe. lxxxiv. p. 526 (1910). 


|] 
| it 
a 
| i it 
if | | 
| | i 
{ 1 
F 
| \ / i 
| \ 
i i L 


Conditions near the Cathode of a Glow Discharge. 915 


positive ion current to F will commence to increase rather 
more rapidly than before, together with any secondary 
emission of electrons from F, and we might thus expect 
the net positive current to F from the gas to start to 
increase more rapidly after F has been made negative to 
the space, and so to give the discontinuous increase in 
slope that is observed to the negative side of B (fig. 2). 
As F is made still more negative, the positive current to F 
from the gas must increase more and more rapidly as the 
potential of F approaches nearer to that of the main 
cathode, whatever the precise nature of the processes which 
render an electrode capable of maintaining a discharge 
when cold. It is, of course, unlikely that conditions are 
quite so simple as this, but some such correlation of figs. 
2 and 3 seems reasonable. 

We have still to account for the breaks Aand C. Ofthe 
two, C is the more pronounced, and might conceivably 
arise in at least three ways. 

(i.) The region of the characteristic between B and C 
might correspond to the transition from a discharge of 
normal type on the collector to one of abnormal type, 
these terms being used to denote incomplete and complete 
utilization of the geometrical surface of the collector F ; 
after the discharge is started, sputtering between F and the 
main cathode probably rapidly makes the two effectively 
of the same material, and so the fact that the break C, 
supposed to mark the onset of the normal discharge, is 
negative to the potential of the main cathode, could be 
supposed to be due to the fact that the collector F is a 
wire of large curvature. Against this view is the fact that 
a space-charge has already been established for F by the 
main discharge, and certain regularities in the position of 
C that are referred to again in section III. 

(ii.) The point C might correspond to the closing in of the 
edges of the shadow behind F towards the main cathode 
until at this stage it ceased to reach to the main cathode. 

(ui). The break at C is of a type which is usually associa- 
ted with the reception of a group of fast electrons in a 
retarding field by a collector. Such a group could arise 
either from the gas, or from the surface of the cathode, 
being produced in the latter case photoelectrically, or by 
bombardment with excited atoms, or positive ions, or by 
some other process of the second kind taking place there. 
It is, of course, possible that fast electrons could come both 
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from the electrode and gas, but we should be inclined to 
favour the former source, particularly on account of the 
general similarity between our curves and those published 
by Oliphant for electrons liberated from a metal surface 
by metastable atoms *. On this interpretation of the 
curves, whatever the space-potential near F, the maximum 
energy of emission of the electrons from the cathode will 
be measured by the potential of C relative to the cathode, 
and many of the electrons would appear to have one half 
the maximum energy inthe group. Thisform of apparatus 
is. however, not suited for study of slow electrons from the 
cathode, since these, unlike electrons emitted with a 
considerable initial velocity, will follow fairly closely the 
lines of the positive ions incident upon the cathode, and 
since the latter have escaped collection by F, the slow 
electrons liberated where the positive ions hit the cathode 
—but not necessarily by the positive ions—must likewise 
largely escape detection. There is, in fact, no upward 
break in the characteristic near the cathode potential. 
The break A may be very tentatively attributed to the 
failure of fast electrons liberated from F by a similar 
process to reach the main cathode in the field retarding 
their motion in that direction. 

We are unfortunately unable at present to decide 
definitely between the alternatives II. and III. 

Curves similar to that of fig. 2 were obtained by 
Aston f for a collecting electrode behind a hole in the main 
cathode, but are even less easy to interpret, both because 
of the disturbance of the main discharge by the aperture, 
and because there are rather too few experimental points 
to be sure of the existence of discontinuities in slope at 


Band A. 


Ill. Results of Discussion. 


The interpretation of the characteristic curves remains 
much less satisfactory than could be desired. We may, 
however, indicate some of the consequences of assuming 
that the tentative conclusions of the last section are 
correct. 

(a) The break B occurs at a potential consistently 
closer to that of the cathode than it would be if the fall in 


* Oliphant, Proc. Roy. Soc. exxiv. p. 228 (1929). 
+ Aston, Proc. Roy. Soc. xevi. p. 200 (1920). 
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potential across the whole region between cathode and 
negative glow followed the parabolic law of Aston *, in 
general agreement with the results of Brown and Thomson. 

(6) There are certain regularities in the position of the 
break C. In each gas there appears to be, on the 
interpretation ii. of the preceding section, a fairly definite 
maximum velocity for the electrons emitted from the 
cathode. In the cases of neon and argon, the break was 
double, the position of the more negative discontinuity 
corresponding to the position of the break in nitrogen ; 
in these cases the less negative break has been assumed to 
be that for the inert gas itself. On the interpretation ii. 
this could be connected with the shadow reaching to the 
cathode at one stage, and to the primary dark space at 


TABLE I. 


Gas No. of breaks Max. energy of | Corresponding critical 
: measured. electrons: volts. potential of gas in volts. 


Nitrogen...... 16 22 +5 28: dissociation to two 
singly ionized atoms. 

Hydrogen ... i 25-544 29-7 : process as for 
nitrogen. 

9 15-5: ionization of atom. 

Neon 6 18-0+3-5 21-5: ionization of atom. 


another. The data obtained for the maximum energies are 
collected in Table I., where the errors given against the 
energies are the root mean square departures from the 
mean. These errors are large, which points to the need for 
better control of the purity of the gas and electrodes than 
we have yet been able to attain. It may be significant, 
nevertheless, that the averages lie, as shown, a few volts 
below important critical potentials of the gas in each case J, 
and it is tempting to suppose that fast electrons are libe- 
rated from the cathode by the atomic processes in question 
taking place there, the electrons losing at least energy 
measured by the work-function of the cathode—4:-4 volts 


* Aston, Proc. Roy. Soc, Ixxxiv. p. 526 (1910). 

+ Data from Geiger-Scheele, Handb. d. Physik, xxiii. article by 
Franck and Jordan for Ar, Ne, No, from Vencov, Comptes Rendus, clxxxix. 
p. 27 (July 1929), for H.. 
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for nickel in vacuum *—in leaving it. The processes 
may then be identifiable with neutralization of the normal 
singly charged ions in neon and argon, simultaneous 
neutralization of two protons, with recombination of the 
latter to form H, in hydrogen, and neutralization of two 
atoms followed by recombination in nitrogen. 

No breaks were found corresponding to C on five 
characteristics taken in oxygen. 


IV. Summary. 

The characteristic curves for a collector very close to the 
cathode of a glow-discharge have discontinuities in slope. 
The most negative of these, which is always negative to the 
potential of the cathode, may be taken to indicate that 
there is a group of fast electrons leaving the cathode 
surface, where it is produced by some process of neutrali- 
zation and recombination of the ions of the cathode dark 
space. The second, which is always positive to the cathode 
potential, is attributed to a change occurring in the mode of 
reception of the current at the space-potential, which, on 
this interpretation, is consistently more negative than if 
Aston’s parabolic law held. The origin of the third, more 
positive break, is most uncertain, but it may indicate 
that there is an emission of fast electrons from the 
collector. 


XCVIII. The Potential Distribution across the Cathode Dark 
Space. By W. L. Brown, M.Sc., Demonstrator, and 
E. E. Tuomson, M.Sc., Musgrave Demonstrator, Physies 
Department, Queen’s University, Belfast F. 


[Plate XVII] 


Introduction. 


ce chief contributions to the study of the electrical 
conditions in the cathode dark space of the glow 


discharge have been those of Astont, Brose §, and 
Geddes ||. 


* Marx, Handb, d. Radiologie, iv. 3, p. 299. 

+ Communicated by Dr. K. G. Emeléus, M.A. 

J Aston, Proc. Roy. Soc. lxxxiv. p. 526 (1911). 

§ Brose, Ann. d. Phys. lviii. p. 731 (1919). 

| Geddes, Proc. Roy. Soc. Edinb. xlvi. 2, p. 186 (1925-26). 
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Aston determined the electric fields at different points 
along the dark space in a wide tube, by passing across it a 
beam of cathode rays from a subsidiary tube kept at very 
low pressure. The deflexion of the beam at different 
distances from the cathode was used as a basis for the 
calculation of the electric field. It was not possible to 
take measurements right up to the cathode, but over the 
region investigated it was found that the field at any 
point was proportional to its distance from the negative 
glow. This means that the space charge is constant 
throughout that part of the dark space considered. The 
tube voltages varied between 265 and 610 volts. 

Brose, using high voltages, 1000 to 4000 volts, in 
discharge-tubes of narrow bore, and with gases at low 
pressures, found the distribution of the field by measure- 
ment of the Stark effect. In common with Aston, he 
found a positive space charge throughout the main part 
of the dark space, but his results indicate that this is not 
constant, and also that there is a negative space charge 
in the region near the cathode. 

Geddes applied Aston’s method to discharges at high 
voltages. His results were similar to those obtained by 
Aston. 

Many attempts have been made to measure the potential 
distribution across the dark space by introducing a third 
electrode, and measuring by means of an electrometer the 
potential taken up by it. These methods fail chiefly on 
account of the differences in mobility and concentration 
of electrons and positive ions, which prevent the wire from 
taking up the potential of the surrounding space *. 

It was intended, on beginning the work described in the 
present paper, to take current-voltage characteristics of a 
collector or third electrode inserted in the cathode dark 
space of the glow discharge, and to attempt the analysis 
of these characteristics by some modification of Langmuir’s 
method f. Previous work on the glow discharge has been 
concerned with the negative glow, Faraday dark space, 
and positive column, in which the characteristics obtained 
can be analysed by thismethod. The form of the charac- 
teristics in the cathode dark space is quite different from 
that obtained in other parts of the discharge, and the 


* Vide, e.g., Emeléus, ‘Conduction of Electricity through Gases,’ 
p. 39. 
+ Langmuir and Mott-Smith, Gen. Elect. Rev. xxvii. p. 554 (1924). 
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interpretation of the curves is difficult. The chief 
difference between the dark space and the rest of the 
discharge is the high value of the field in the dark space, 
almost the whole fall of potential applied to the tube 
taking place across it. 

The data obtained were plotted in different ways and the 
resulting curves examined. The relation V against log 
(c+const.) gives a straight line with a break at a voltage 
between that of anode and cathode. This voltage 
approaches that of the anode as the collector is moved 
nearer the negative glow. We compared the values of this 
potential with observations taken by the method to be 
described later, but although agreement was very close 
near the cathode, it became less exact as the collector 
approached the negative glow. 

The characteristics proving unsatisfactory, we turned our 
attention to another phenomenon. That there is a shadow 
cast towards the cathode by an obstacle placed in the 
cathode dark space has long been known. However, so far 
as we are aware, the effect of applying a potential to the 
obstacle has not been investigated. We carried out a 
series of observations, using obstacles (or collectors), of 
different shapes. Metals were used for all electrodes. 

The effect of using a large circular collector, 8 mm, in 
diameter, mounted parallel to a cathode 2:5 cm. in diameter 
was very striking. When the collector was insulated, 
there was a circular patch cut out of the cathode glow 
underneath, and of greater area than the collector, while 
there was a very small circular patch cut out of the 
negative glow. On reducing the potential of the collector 
relative to the cathode, the cathode patch diminished while 
that in the negative glow increased in size, passing through 
the intermediate stage when both were equal in area to the 
collector, to the condition in which there was a large patch 
cut out of the negative glow, while the cathode glow was 
almost undisturbed. 

The results of observations, using as collector a wire, 
0-4 mm. in diameter, mounted parallel to the cathode 
surface, may be summarised as follows :— 


(a) The potential applied to the collector by an external 
battery, necessary to maintain the shadow towards 
the cathode in the same form as when the collector 
is insulated, is always near the anode potential. 
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The form of the shadow is shown in Pl. XVII. 
fig. 1. It will be seen that the shadow is much 
wider than the wire. 


(6) If the applied potential is reduced, the width of the 
shadow decreases. When the shadow to the 
cathode side has the same width as the wire, 
another shadow has appeared on the anode side, 
also of the same width as the wire. This is shown 


in Pl. XVII. fig. 2. 


(c) Further reduction of the applied potential leads to 
the complete disappearance of the shadow to the 
cathode side, while the shadow on the anode side 
increases in size. 


(d) As the potential approaches that of the cathode, a 
bright layer appears on the cathode side of the 
wire, while on the anode side the shadow is almost 
parabolic, with a bright layer round the edges. 
XVIL. fig. 3.) 


The photographs shown were obtained in argon, with a 
tube voltage of 470 volts; the potential of the 
wire relative to the cathode was 228 volts in 
Pl. XVII. fig. 2 and 62 volts in Pl. XVII. fig. 3. 
In Pl. XVII. fig. 1 the wire was insulated. The 
dark space was | cm. in width in all cases. The 
photographs have all been enlarged two diameters, 
and only a small part of the cathode is shown. 


The above facts can be explained qualitatively by 
considering that the potential of the collector when the 
shadows thrown to the two sides of it are equal is the 
potential of the surrounding space. This is reasonable, 
because the disturbance of the discharge, caused by the 
presence of the collector, appears to be least at this 
potential. Apparently positive ions moving towards the 
cathode, and electrons moving towards the anode, in 
paths parallel to the axis of the tube, are not deflected. 
When the potential of the wire is near that of the cathode, 
7.é., negative to the surrounding space, we should expect 
electrons moving towards the anode to be deflected away 
from the collector and positive ions moving towards the 
cathode to be deflected towards it, so that the discharge 
would be more concentrated than normal between collector 
and cathode, and less concentrated on the other side of the 
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wire. The converse would be true when the potential of 
the collector was near that of the anode. 

From such considerations, we have concluded that the 
potential of the collector, when the shadows thrown to the 
two sides of it are equal in width, is the potential of the 
surrounding space, and we have used this as a basis for 
measuring the potential distribution across the cathode 
dark space. 


Preliminary Observations. 


Before considering the experimental arrangements for 
carrying out measurements of potential, it is of interest to 
mention some observations made in connexion with the 
shadow thrown to the cathode side of the collector. It was 
found that if the collector were charged to some potential 
at which a shadow was cast towards the cathode, and the 
discharge left running for a considerable time, perhaps 
an hour, the part of the cathode in shadow became sharply 
defined by a curious bluish marking, bounded by sharp 
lines. A photograph of markings on an iron cathode 
surface, obtained in this way, using a molybdenum 
collector, is shown in Pl. XVII. fig. 4. 

The marking obtained was more intense with longer 
exposure to the discharge, and with higher current 
densities. If the potential of the wire is changed after a 
marking has been formed, so that there is no longer a 
shadow on the cathode side, the marking will disappear 
after a time. 

It is well known that the glass surrounding the cathode 
of a discharge-tube becomes sputtered after a time de- 
pending on the intensity of the discharge, and the cathode 
material. This is due to atoms of the metal being 
liberated from the cathode by the impact of positive ions 
and being deposited on the neighbouring glass surfaces. 
We may attribute the marking described above to the 
same cause. Atoms of metal pass temporarily into the 
discharge, and most of these are returned to the cathode, 
the remainder depositing on the tube walls. We suggest 
that those atoms which deposit on that part of the cathode 
in shadow tarnish the surface, while atoms returning to 
the bright part of the cathode are prevented from tarnish- 
ing it by the continual bombardment of positive ions. 
This explanation is borne out by the fact that the 
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marking disappears when the discharge is allowed to fall 
on it. 

We should also point out here that Oliphant * has 
obtained disintegration of the cathode, due to metastable 
helium atoms. 

The form of the markings produced by using different 
types of collectors was studied next. In every case the 
shape of the marking corresponded roughly with that of the 
collector, small obstacles giving the most sharply defined 
markings. 

The sharp lines on the edges of the markings (Pl. XVII. 
fig. 4) might possibly be due to a kind of focussing effect. 
There will be a tendency for positive ions to move across 
the dark space more or less perpendicularly to the cathode, 
under the influence of the large electric field. The path 
will, of course, be zigzag, but the main direction will be 
normal to the cathode. A positive ion moving in this 
direction, and ‘passing near the wire when it is charged 
positively to the space, will be deflected away from the 
wire. If the number of collisions is fairly large, the path 
of the ion will rapidly return to a direction parallel to the 
original direction, and the effect of the wire will be merely 
to displace the particle sideways. Other positive ions 
passing the wire at greater distances will be displaced 
less, and this will tend to lead to a sharp line of demarca- 
tion between the shadow and the normal discharge. If 
the number of collisions is not large, we should expect a 
more diffuse edge to the marking. 

An interesting case studied was that a) a wire intro- 
duced from a side tube, and bent so that its free end was 
furthest from the cathode. The marking produced on the 
cathode surface is shown in Pl. XVII. fig. 5. It will be 
observed that the marking is not parallel sided, as in 
Pl. XVII. fig. 4. but has a maximum width, diminishing 
to a point towards the free end of the wire. We should 
expect that the point where the marking is the same width 
as the wire corresponds to that point of the wire which has 
the same potential as the space surrounding it. Bearing 
in mind ‘the facts mentioned in the Introduction, it is clear 
that the part of the wire nearest the negative glow may 
not cast any shadow at all. 

On the cathode of Pl. XVII. fig. 5 the marking was 


* Oliphant, Proc. Roy. Soc. exxiv. p. 228 (1929). 
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brown, the collector in this case being of copper-clad nickel 
iron, and the cathode ofiron. The potential of the collector 
relative to the cathode was 240 volts. 

Wehnelt * and Villard} found that the cathode became 
tarnished with a film of oxide on that part of it exposed to 
the discharge. They put this down to the oxidizing action 
of positive rays. The part of the cathode in shadow 
remained bright, while the edges of the deposit were not 
sharp. The differences between their results and those 
considered above must be due to the fact that they worked 


with an induction-coil discharge in air at low pressure, and 
a copper cathode. 


Experimental Arrangements. 


The discharge-tube was run off a battery of small 
accumulators, capable of giving a maximum voltage of 
660 volts. The first circuit used is shown in fig. l.a. 
There was a wire resistance of 1100 ohms in series with the 
tube, to prevent the glow discharge from passing into an 


* Wehnelt, Ann. d, Phys. \xvii. p. 421 (1899). 
+ Villard, Comptes Rendus, cxxvi. p. 1564 (1898). 


Fie. |. 

| 

| 

(6) 
(a) 
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arc. The potential applied to the collector C could be 
adjusted by varying the voltage of the battery A, and fine 
control was obtained by means of the potentiometer P, 
which had a maximum range of 60 volts. Voltages were 
read on the Weston high range voltmeter V. When 
taking current-voltage characteristics, Pye mirror galvano- 
meters were brought into the circuit. 

With this arrangement there was the difficulty that as 
the voltage available for the battery A was often less than 
half that of the main battery, there was a break in the 
series of readings when the voltage which was required 
for the collector was nearly half-way between anode and 
cathode potential. In order to avoid this, in some of the 
later work the battery A was abandoned, and the lead to 
it was taken to a plug, which could be connected to any 
desired point along the main battery (fig. 1,6). The 
potentiometer P was retained, and with this arrangement 
any voltage between anode and cathode potential could be 
applied to the collector. 

The only possible objection to this arrangement was the 
fact that the tube and collector circuits were more closely 
coupled than when the first circuit was used. However, 
it was found that the results obtained with the two circuits 
showed no differences. 

The discharge-tubes were exhausted by a Gaede steel 
mercury vapour pump, backed by a Hyvac pump. The 
tube section of the apparatus was cut off by a tap from the 
McLeod gauge and the pumps, and arrangements were 
made to introduce the different gases as required. Argon 
was purified by running an arc between calcium electrodes 
in a bulb sealed to the apparatus near the discharge-tube. 
For neon, this bulb was replaced by a liquid-air trap 
containing charcoal. The nitrogen used was generated by 
exploding sodium azide in a side reservoir, and was con- 
sidered sufficiently pure without treatment. In all cases 
mercury vapour from the pump and the McLeod gauge, 
and hydrocarbon vapours from the taps, were kept out of 
the discharge-tube by means of liquid-air traps. Before 
making observations the discharge-tube and the glass- 
work near it were well baked by the gas flame for some 
hours. While observations were being made a liquid-air 
trap was always used. 

These precautions allowed us to work with the gases in a 
fairly pure state. In the argon, after the calcium arc had 


Phil. Mag.8S.7. Vol. 8. No. 53. Suppl. Dec. 1929. 3Q 
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been running for several hours, no trace of nitrogen could 
be detected spectroscopically, although the Balmer lines 
of hydrogen were just visible. Both the neon and the 
argon showed a well marked primary dark space. The 
sample of neon used contained 2 per cent. of helium, other 
impurities being removed by the charcoal tube. 

The discharge-tube, of soda glass, was 3 cm. in internal 
diameter. The cathode was a steel tube 4 cm. long, and 
fitted it fairly closely. The end of the steel tube was 
closed by a disk of iron or nickel, which acted as the cathode 
face. A small hole was drilled through this disk, as near 
the edge as possible, and through this was passed a thin 
glass rod, 4 cm. long; this rod carried the anode, which 
was made of light nickel sheet. The distance between the 
electrodes was therefore constant. The glass rod was so 
fine that it did not disturb the discharge appreciably. 

The collector was a fine steel wire, 0:4 mm. in diameter, 
carried by a glass sheath which was introduced through a 
side tube. Steel was used as it was the only material at 
hand stiff enough to remain quite straight when supported 
only at one end. 

The discharge-tube was set up with a slight tilt towards 
the cathode end of the tube. When a set of observations 
was being made, it was begun with the cathode as close to 
the collector as possible, so that after each reading had 
been taken, a slight tap on the tube caused the cathode to 
move a little farther from the collector. When a series of 
observations was completed the cathode could be brought 
back to its original position by means of an electromagnet. 
The distance between the cathode and the collector was 
measured with a travelling microscope, adjusted so thay 
its axis was parallel to the collector. 


Results. 


The observations of the potential at which the shadows 
on both sides of the wire were equal could be reproduced 
to within 5 volts throughout the dark space except near its 
cathode and negative glow boundaries, where the shadow 
on one side became so short that it was more difficult to 
determine the point of balance. This represents quite a 
high degree of accuracy when we consider that a wire 
0-4 mm. in diameter placed in a field of 300 volts per cm., 
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such as exists in the dark space, will be lying in a region 
across which there is a fall of potential of 12 volts. 

When the collector approached close to the negative 
glow edge, it was found that not only did the shadows cast 
by the wire alter in form but the edge of the glow itself 
became somewhat distorted, as the potential applied to the 
wire was altered. However, in spite of this, no great 
difficulty was experienced in obtaining the balancing 
potential, since the glow assumed its normal outline when 
the shadows were balanced. By “normal outline,’ we 
refer to the contour of the glow before the collector was 
moved near. 

Another point of some importance arises here. If the 
negative glow is not parallel to the cathode but presents a 
convex surface towards it, as is often the case at low 
voltages when the pressure is low, then we must consider 
that the equipotential surfaces across the section of the 
tube are curved also. In such cases the shadow cast by 
the wire will vary in width, and it is of importance to know 
the point of the wire for which we have balanced the 
shadows. On one occasion an accidental tap on the tube 
caused the glass sheath carrying the collector to move, so 
that its end passed into the path of the discharge. The 
shadow cast by it towards the cathode (it will be remem- 
bered that an insulated obstacle of any description casts a 
shadow towards the cathode) was wider than that cast by 
the wire when near the balancing potential. Although on 
the side of the tube remote from the eye, this shadow could 
be seen quite well in addition to that thrown by the wire. 
The shadow thrown by the glass may be seen in Pl. XVII. 
fig. 6. 

th such cases the point of balance becomes rather 
indefinite, and the shadow of the same width as the wire 
seems to move along it as the potential is slightly altered ; 
so that the shadows thrown by all parts of the wire 
influence the observer in deciding the balancing potential. 
That there is not any great deviation from plane equi- 
potential surfaces in most cases, however, seems to be 
indicated by the parallelism of the marking on the cathode 
shown in Pl. XVII. fig. 4. 

On completing a series of observations over the dark 
space by the shadow method, we usually estimated the 
space potential in the negative glow by moving the collector 
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into it and increasing the voltage applied until the electron 
current to the collector began to increase rapidly. Points 
on the curves which were thus obtained are marked with 
an L. 

Data obtained with different gases under various 
conditions of pressure and voltage are given in Table I. 


V: Volts Fig. 2 


400 | 
360 


320 


(6) "2 “4 1-0 2 4 
Centimetres from Cathode 


The vertical lines indicate the positions of the negative glow edge. 
Data for these curves will be found in Table J. Nos. 6, 7, and 8. 


Some typical graphs drawn from these data are shown 
in figs. 2 and 3. 

One of the curves in fig. 2 is lettered for the purpose of 
explanation. It risés smoothly to the point C, where the 
collector passes into the negative glow. The vertical line 


- 
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through the curve indicates the position of the cathode 
edge of the negative glow. There is a break in the curve 
near its lower end, and the middle part, AC, when extra- 
polated, does not pass through the origin but cuts the axis 
at a distance of the order of half a millimetre from it. 
This initial part of the curve cannot be found accurately 
by the shadow method, as the cathode shadow here is so 
short; but the fact that AC does not pass through the 


V: Volts Fig. 3 


320 


280 


240 vA 


40 


(0) “4 ‘6 1-0 1-2 
Centimetres from Cathode 


The vertical line indicates the position of the negative glow edge. 
Data for this curve will be found in Table I. No. 26. 


origin confirms the supposition that there is a layer of 
negative space charge near the cathode surface, as shown 
by Brose’s curves for strongly abnormal cathode fall, and 
that the effective cathode is the edge of this layer. 

P. M. Morse* has recently given a comprehensive 
theory of the glow discharge, which is the most satisfactory 
yet proposed, and which it is desirable to test from these 


* Morse, Phys. Rev. xxxi. p. 1003 (1928). 
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data. Starting from Poisson’s equation, and making use 
of Townsend’s relation, 


a=pNe . 
he arrived at the equation 
V=4D(1 —c) —4D(1 + 8D N(1 —c) EH, px/3. (2) 


for the fall of potential across the cathode dark space. 
In the above equations : 


a=the number of ions formed per cm. of path by 
an electron moving under a field E, 


p=pressure of gas, 
1/(1—c)=fraction of total discharge current carried by 

positive ions. 

v=ionization potential of the gas, 

V=space potential at any distance z from cathode, 

D=diffusion coefficient, 

K,=a constant. 

Morse also gives the relation, 

4D(1—c)=V’ approximately, . . . . (3) 


where V’ is the normal cathode fall of potential. 
From (2) and (3), 


V=V'(1 —e~PNz3 INE px/3). 


which may be written 


. . . . . (4) 
where k=2N/3, and K=kEa. 
Hence, 


where S=dV/dz. 
When «=0, 


In the analysis of our curves on this basis, the main part 
of the experimental curve AC (fig. 2) has been extra- 
polated to cut the axis, and the point of section taken to 
be the effective cathode (x=0). The slope of the tangent 
drawn at this point gives a value of §. By substituting 
this value in equation (6), the value of k+K is found, 
p and V’ being known. 
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Putting V=V’ in equation (4), we have 
—K pe. 


The value of « for V=V’ is found from the experimental 
curve and substituted in this equation. The values of k 
and K can then be found, either by trial or graphically, 
by drawing the curve 


where n=opx and €=kpz, 


Fig. 4. 


V: Volts 
600 


ia} 


co) 
S 


Space Potential re/ative to Cathode 
w 


(6) 


Centimetres from Effective Cathode 


The distance between cathode and effective cathode is 0'U2 em. 
The data for this evrve will be found in Tables I. and II. No. 17. 


o@ px is known, and so from the graph we can get kpx and 
hence k. The values of k, K, V’, and 7 are now all known, 
and are substituted in equation (4), and the curve of V 
against computed and drawn. 

Owing to the nature of Morse’s equation, involving as it 
does three terms including an exponential, there is some 
latitude in choosing the value of V’ (see below), and of 
course k and K vary in consequence. This, taken in con- 
junction with the difficulty of drawing a tangent at the 


| 
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effective cathode, rendered it necessary to test the values. 
obtained in each case, by drawing the theoretical and 
experimental curves. If the agreement was not satis- 
factory the values were recalculated. In all cases 
excellent agreement was finally obtained between the 
calculated curve and the experimental points. One of the 
curves is shown in fig. 4, where the line is the computed 
curve, and the dots are the experimental values. 


TaBLeE II. 
No Wi k K. p Ey Gas 
5 270 15-6 45 ‘107 288 ) 
6 290 15-9 4:7 +295 
fi 227 14-0 3:3 
8 340 16-0 4:6 
Argon. 
9 230 26:3 1-7 -055 +293 
10 275 28-7 7-0 243 
11 330 34-0 45 
12 376 33-4 4-46 Ms 133 | 
17 480 48-8 7:3 -026 149 
18 385 46-4 7-6 163 
19 515 41-9 8-1 5 193 
20 206 17-1 1-14 -123 067 ) 
21 248 15:9 2-5 
24 341 22-5 5-4 O74 -239 ‘Neon. 
25 261 4:3 :207 
27 288 21-7 4:5 


V' is in volts and p in centimetres of mercury. 


The cathode fall of potential can therefore be repre- 
sented by Morse’s equation between the effective cathode 
and the negative glow. At the cathode edge of the glow 
the equation breaks down *. 

When several of the potential-distance curves for a given. 
pressure and gas are drawn on the same scale, we find that,. 


Md Morse gives the value of the field at the point of breakdown as: 
V'kp. This is incorrect, as it involyes the assumption that at the end 
= the dark space e~4p is negligible, which we know (Table II.) is not 
he case. 
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as we should expect, the field at the end of the dark space 
decreases as the tube voltage decreases, the point of 
breakdown of the dark space curve on passing into the 
negative glow becoming less sharp as the cathode fall of 
potential approaches the normal fall for the material of the 
cathode (see, for example, fig. 2). 

The values found for the quantities in the equation 
are given in Table II. 

In the table the pressure p is in centimetres of mercury 
and V’ in volts. 

V’ is not the normal fall of potential, as it was found that 
the normal fall was too low to admit of the equation fitting 
the curves. The value used was 11/12 of the voltage 
across the tube, which was less than the actual fall across 
the cathode dark space as shown by Langmuir charac- 
teristics taken in the negative glow, but was greater than 
the normal cathode fall in all cases except those at very 
low voltage. The fraction 11/12 was found by trial to be 
a suitable value. 


Discussion. 


From Table II. it may be seen that k is nearly constant 
for a given gas at a fixed pressure. Beyond this there is 
nothing definite, but there are two possibilities. The first 
is that kp is a constant for a given gas. For runs 5-8, 
kp has an average value of 1-65, while for runs 9-12 the 
average value is 1:68. In nitrogen the average value of 
kp is 1:19. The values of & in neon, however, change 
very little with pressure. 

The second possibility is that k is nearly constant for a 
given gas, and almost independent of pressure. In this 
case we should have to assume that some impurity got 
into the gas when the pressure was changed between runs 
8 and 9. It is highly probable that neither of the above 
suggestions is strictly true, but that, as is indicated below, 
k is influenced considerably by both pressure and the 
degree of purity of the gas. 

Since & is nearly constant, the value of N (=3k/2) is also 
nearly constant for a given gas at a fixed pressure. Accord- 
ing to Townsend, N is the number of ionizing collisions 
made by an electron per centimetre of path, in a field of 
1 volt per cm., the pressure being 1 mm. of mercury. 
Since it is known that pressure influences the rate of 
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removal of excited atoms, we might expect some variation 
from Townsend’s values *. N is also dependent on the 
purity of the gas. 

The observed values of N, expressed in Townsend’s 
units, are of the same order of magnitude as those given by 
him. Values for different gases are given in the following 
table :— 


III. 
Gas. Argon. Nitrogen. Neon. Helium. 
Townsend’s value of N... 13:6 12-4 — 2-8 
Observed value of N ...... 3-4 6:8 2-9 


Observed values of N were obtained from average values of / in Table II. 


Townsend does not give a value for neon, but the value 
he gives for helium is close to that found for neon. Con- 
sidering also that we measure the density of the gas only 
indirectly by means of the pressure, the agreement in this 
table is good evidence in favour of Morse’s theory. 

The value of K shows little regularity and is not of 
obvious theoretical significance. Its value is always 
considerably less than that of k, their ratio being roughly 
as1:5. EH, therefore does not vary a great deal from 0-21. 

The total space charge per sq. cm. between two layers 
parallel to the cathode surface whose distance apart is 
(a—b) can be determined from Poisson’s equation. 

Since 


d? 
we have 


and therefore 


podem (ae), 


If we measure the slope of the potential-distance curve 
at two points distant a and b respectively from the cathode, 


* Townsend, ‘ Electricity in Gases,’ p. 295 1915). 
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we can determine the mean space charge per c.c. between 
these points. The slope can be determined either from 
Morse’s equation by differentiating and substituting the 
appropriate constants obtained from our experimental 
data, or directly by drawing a tangent to the experimental 
curve. 

Referring again to figs. 2 and 3, it is evident that we 
can do no more than suggest the form of the curve between 
the cathode and the experimental points. For this reason 
it is impossible to measure the total space charge between 
the cathode and the negative glow, although the variation 
of space charge from the beginning of the experimental 
curve to the negative glow can be studied. 


TaBLE IV. 
1 2 3 “4, 5 6 
No. 5 |* Beka 510 454 402 362 324 294 270 
(kere 169 -143 -121 -102 -087 -073 -062 
2 3 “4 5 6 7 8 ‘9 
594 508 486 373 323 284 247 221 198 
254 +210 -175 -145 -120 -100 -083 -069 -057 
Noot7..5 1 629 509 413 341 283 240 206 
181 -140 -109 -084 -065 -051 -039 
Pa nssleis 1 2 3 4 5 6 7 8 9 
476 404 347 300 258 226 201 179 161 
+164 -135 -111 -091 -075 -062 -051 -042 -034 
a is in centimetres. E in volts per em. p in €.S.u. per c.c, 


From Morse’s equation we have, by differentiating, 


dV 
E= +Kp) (5) 


= 

When the appropriate values of V’, k, K, and p are sub- 
stituted in the above equations, we can draw up a table 
showing the variation of field and space charge as « is 
increased. Such a table of values obtained from the data 
for several of the given curves is shown in Table IV. 
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Values of field obtained by measuring a series of slopes 
on the experimental curve for run 17 are plotted against 
the distance from the cathode in fig. 5a. In the same 
figure, curve (6) shows the corresponding space charge 
distance curve, the data for which were obtained by 
measuring the slopes at different points on curve (a). 
Leaving out of consideration the region between cathode 
and effective cathode, it will be observed that the field 
has its maximum value a short distance out from the 
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cathode, and falls gradually to less than one-third of this 
value at the cathode edge of the negative glow, when it 
drops rapidly to zero in the negative glow. Curve (6) 
indicates that the positive space charge concentration 
prevailing throughout the greater part of the dark space 
decreases almost uniformly to less than half its initial 
value between the effective cathode and the negative glow 
edge. Between this point and the middle of the negative 
glow the positive space charge varies rapidly, at first 
increasing to from 4 to 6 times its average value over the 
dark space, and then falling to zero. 
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Regarding the space between the cathode and the 
effective cathode, it is evident that there must be at least 
one point of inflexion on the potential distance curve ; 
that is, the space charge curve must pass through zero. 
This is indicated by the dotted line. There may, of course, 
be two points of inflexion, giving a positive space charge 
on the cathode surface, with a negative space charge 
farther out, and then the main dark space. This is 
discussed later. 

The sudden change in space charge concentration in 
curve (b) marks the point where Morse’s equation ceases 
to hold. It is impossible to mark off the boundaries of 
this layer of high positive space charge accurately in many 
cases since the point where the curve turns over to run 
parallel to the distance axis is not very definite. The 
average estimated thickness of the layer is 1 mm., while 
the space charge concentration in it is approximately three 
times the mean space charge over the region to which the 
equation is applicable. 

Miss N. M. Carmichael * in some recent work in this 
laboratory, has shown that on the cathode side of the 
negative glow maximum, and between it and the feeble 
luminosity in the cathode dark space, there exists a dark 
layer. This dark space has not yet been observed in any 
gas other than oxygen. Voltage-current characteristics 
taken near this dark space show a change in form as the 
collector is moved through it. It is probable that if a 
similar dark space should subsequently be found in other 
gases, it will correspond in position to the breakdown of 
the cathode dark space potential-distance curve. 

Emeléus and Carmichael tf have given reasons for 
supposing that there is a positive space charge on the 
cathode surface. Now a positive space charge is indicated 
when the potential distance curve is concave to the 
distance axis. If this is so, consideration of the curves 
for argon in fig. 2 shows that the curve between the cathode 
and the experimental points should have two points of 
inflexion. All the curves drawn from the observations in 
neon, of which fig. 3 is an example, show at least one 
point, that nearest the ¢athode, which does not lie on the 
smooth line through the others. 


* Carmichael, Phil. Mag, viii. p. 362 (1929). 
t Emeléus and Carmichael, Phil. Mag. v. p. 1089 (1928), 
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The dotted lines in figs. 2 and 3 show the most obvious 
ways of drawing the curves between the cathode andthe 
experimental points, taking into consideration only the 
actual distribution of the points themselves. At first 
sight it would appear that the observations in neon would 
justify us in drawing the initial parts of the curves in such 
a way as to be in accord with the idea of a positive space 
charge immediately in front of the cathode surface. The 
absence of a corresponding point in argon and nitrogen 
would, however, be difficult to account for. 

Again, when we consider the difficulty of making an 
observation of balancing potential close up to the cathode 
surface, a task which, while troublesome in all cases, is 
especially so in the case of neon, in which a primary dark 
space is visible (for example, Pl. XVII. fig. 6), we are 
inclined to attribute the position of the first point to an 
experimental error. The presence of the primary dark 
space means that the shadow to the cathode side of the 
collector is shorter by the width of the dark layer, with 
the result that the first reliable reading, other things being 
equal, is somewhat farther from the cathode in cases where 
a primary dark space is visible. On the other hand, the 
point always appears on the same side of the main curve, 
seeming to indicate that any experimental error in deter- 
mining its position is not of great account. 

Taking into account all these facts, we do not think our 
results in this region afford decisive evidence either for or 
against the existence of a positive space charge on the 
cathode surface. 

The presence of a layer or layers of space charge near the 
cathode surface must almost certainly be associated in some 
way with the production of the Aston or primary dark 
space. From the fact that the effective cathode had to be 
taken at some distance out from the cathode in most 
cases, whether a primary dark space was present or not, we 
think we are justified in supposing that some, at least, of 
the conditions necessary for the production of a primary 
dark space must obtain in all glow discharges. Provided 
that a similar argument holds in the case of the dark space 
near the negative glow, it appears that the potential- 
distance curve can be expressed in the form of the equation 
given by Morse, if we consider only that part of the cathode 
dark space proper lying between these two more intensely 
dark regions. The effective cathode is not, however, at 
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the cathode potential, and so there must be some modifi- 
cation of the curve in the transition region between the 
primary dark space and the main dark space. 

The sheath surrounding the collector in neon is very 
interesting, as it is the first indication we have had of 
anything resembling the ordinary positive ion sheath in 
the cathode dark space. The characteristic curve affords 
no evidence either for or against the existence of a positive 
ion sheath round a collector in the dark space, for, 
although the curve obtained at potentials negative to that 
of the surrounding space has a form compatible with the 
presence of such a sheath, we know that the grid charac- 
teristic of a valve, in which the grid is not limited in its 
influence by a positive ion sheath, has a similar form. <A 
photograph taken in neon, when the collector potential was 
negative to the surrounding space, is shown in Pl. XVII. 
fig.6. A further investigation of sheath formation in the 
dark space is contemplated. 

To summarize our results: by means of the shadow 
method, a continuous distribution of potential has been 
found from close upon the cathode to the negative glow. 
This in itself is important, as we do not have to account for 
an anomalous distribution of potential of the type of the 
“ Kathodensprung”’ either near the cathode or at the 
boundary between the cathode dark space and the negative 
glow, where our curve passes continuously into the 
Langmuir curve of potentials for the negative glow itself. 
Although, as indicated previously, there is.some latitude 
in the choice of the values of V’, k, and K in Morse’s 
equation, yet the agreement of the values of N with those 
obtained by Townsend is strong evidence in support both 
of the substantial accuracy of Morse’s theory of the dark 
space and of our methods of observation. Nevertheless, 
owing to the complex nature of the problem of the cathode 
dark space, we have thought it best to give in detail a large 
part of the data obtained, so that, should any further 
theory be worked out, a preliminary test will be available. 


Summary. 


1. Details are given of some markings observed on the 
cathode when an obstacle, or collector, is placed in the 
cathode dark space. 
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2. Anew method of determining space potentials in the 
cathode dark space is put forward, based on the behaviour 
of shadows thrown by a wire inserted init. Values of the 
potential obtained in the dark space by this method and in 
the negative glow by Langmuir’s method give a con- 
tinuous curve from near the cathode potential right into 
the negative glow, for cathode falls in potential between 
220 volts and 560 volts. 

3. The results obtained conform to a theoretical expres- 
sion given by P. M. Morse, connecting space-potential and 
distance from the cathode, provided we measure distance 
from an “‘ effective cathode,” about half a millimetre from 
the cathode surface, and choose suitable constants. 

4. The constants in Morse’s equation involve N, the 
number of ionizing collisions of an electron per centimetre 
of path under given conditions. Values of N obtained 
from our data are shown to be of the same order as 
Townsend’s values. 

5. Space-charge and electric field conditions over the 
dark space are discussed. 


We are indebted to Dr. K. G. Emeléus for suggestions 
made during the course of the investigation, and for help 
in the preparation of this paper. 


Key to Plate XVII. 


Figs. 1-3 show the torm of shadows thrown by a wire inserted in the 
cuthode dark space of the glow discharge and charged to different 
potentials. 

The photographs were obtained with argon, with a tube voltage of 
470; the potential of the wire relative to the cathode was 228 volts in 
fig. 2 and 62 volts in fig. 8. In fig. 1 the wire was imsulated. The dark 
space was lem. in width. The photographs have been enlarged two 
diameters and only a small part of the cathode is shown. 


Figs. 4 and 5 show markings produced on that part of the cathode in 
shadow. In fig. 4 the molybdenum wire causing the shadow was 
parallel to the cathode surface and at a distance from it of 0:48 em. 
The cathode dark space was 0'7 cm. in width. The wire was eba r 2 
to a potential 880 volts positive to the cathode, while the tube voltage 
was 440 volts. The diameter of the wire was 0:01 cm. and the width 
of the marking 0:11 cm. 

In fie. 5 the wire was of copper-clad nickel-iron and was inclined to 
the cathode surface at an angle of 30°. ‘The potential applied to the 
wire was 240 volts, while the tube voltage was 490. 


Fig. 6 shows a bright sheath round the wire formed when its potential 
was nevative to the surrounding space. The discharge in this case was 
in neon and shows a well-marked primary dark space. 


XCIX. Problems of determining Initial and Maximum Stresse 
in Ties and Struts under Elastic or Rigid End Constraints 
By W. H. Brooks, B.Se., Ph.D.( Eng.) Lond.* 


Introduction. 


i practice it may often be very desirable to know the 

minimum and maximum stresses to which any particular 
tie or strut may be subjected when forming part of a given 
structure and consequently stressed by irremovable dead 
loads. 

The stresses taken in designing the tension and compression 
members of many framed structures are commonly based 
upon the assumptions of pin-jointed ends and negligible 
deformation of the configuration of the structure when 
loaded, but faulty workmanship, coupling adjustment, or a 
slight alteration or configuration, due to the loads or to the 
subsidence of a support, may set up stresses only deter- 
minate by direct experiment, and in extreme cases of 
subsidence, even reversals of stresses in certain members 
may result. 

Methods for determining the stresses in thin wires and 
slight swaged rods—as the bracing wires and raf-wires of 
aeroplanes—are known “ * ©, but the application is limited 
to members whose transverse dimensions are small as com- 
pared with their lengths. Instruments are in use ™) to (9 
which are applicable to members of larger cross-sectional @Yea ; 
but these indicate only incremental stresses by the direct 
measurements of strains in the usual cases where the initial 
loads cannot be reduced to zero, which fact often appears to 
be overlooked. 

The objects of the author are (a) to show how, by using 
any suitable one of the three methods introduced later, it is 
possible to find the stress (due to irremovable loading) 
existing in a tie or strut of uniform cross-section having its 
ends constrained by hinges or clamps; (0) to indicate by 
Stress Charts how most of the equations established may be 
readily solved ; and (c) to experimentally verify some of the 
cases discussed. 

The following methods are confined to the determination 
of initial direct stress, set up by a longitudinal end load F, 


* Communicated by the Author. Section of Thesis approved for the 
Degree of Doctor of Philosophy in the University of London. (Copies 
of the complete Thesis are deposited in the University Library, South 
Kensington.) 
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to which—at any particular cross-section—would be added 
any secondary bending stresses due to lateral loading or to 
eccentricity of end load, both of which may co-exist with P. 
Methods for calculating such secondary stresses are fully 
discussed in standard books on materials. 


The direct stress fa at any section = A? where A is the 
cross-sectional area, and fz will have its minimum or 
maximum values according as P has its initial and lowest 
value P, or its highest value when increased by incremental 
temporary loading ; 


P 
fa min. =| 


AY 
Ja max. = (P,+incremental P)/A. 


In this latter expression the magnitude of the second 
term may be determined by the strain gauges in use ; hence, 
when Pa is known, fj max. is also known. Alternatively, 
fa max. may be found direct by the application of one of the 
following methods to the member under maximum load 
conditions. 


Principles and Preliminary Argument. 


The underlying principles herein involved in the deter- 
mination of these direct stresses are briefly as follows :— 

The member, whether tie or strut, is very slightly flexed 
by the temporary application of either lateral known forces 
or known couples of magnitude sufficient to give the smallest 
deflexions or slopes that can be accurately measured. The 
test readings taken are then treated as described under 
‘General Procedure.” 

These deflexion and/slope effects are known to be pro- 
portional to the disturbing causes in the isolated cases of 
members having their ends either simply supported or 
encastré with no end load. They are also proportional—as 
will be seen from the equations established later—when the 
end load remains of constant magnitude during the flexing 
of the member. When, however, the end loads are applied 
through strong end constraints, the process of flexing may 
be expected to result in altered reactions of the constraints 
with consequent alteration of the direct stress induced in the 
member, and, if the initial end load is large (in the case of 
a tie it may be many multiples of the crippling load for the 
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same meinber when used as a strut), this effect may be 
apparent for comparatively very small deflexions, 

Consider a tie A, B under initial tension = P, (fig. 1), 
caused by the reaction of elastic end constraints, for which 
the springs A, C and B, D are simple substitutes. Cand D 
are at a fixed distance apart. When a slight flexing is ia 
to A, B, as shown exaggerated in fig. 2, the effect will be a 
reduction of the original length to, say, (/— so forcing the 
springs to a a aan and thereby i increasing Pa to, | say, 

P,+é6P=P in fig. 2. 

Thus, in the case of a tie under end loads applied through 
elastic constraints, the effect of flexing is increased tension. 
Conversely, in the case of a strut, the effect of flexing is 
reduction in compression, as may readily be seen by replacing 


A, Cand B, D by compression springs. 


Fig. 


I 


Tn the case of a tie or strut having its ends fixed by rigid 
constraints, the flexedjaxis length is “greater than the length 
(1) of the unflexed axis, and, assuming the yield of the end 
clamps to be negligible, it may be shown that the longi- 
tudinal strain of the flexed axis is increased by 


21 (3 de, 


with consequent increase in the tension, or decrease in com- 
pression, according to whether the member is a tie or a 


strut : 
EKA (! 
= dz, 


where the symbols have their usual meaning, and the origin 
is taken at one end of the member. 


4 t 
Tm 9) 
LALAAA™L 
| 
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In all cases it should therefore be assumed at the outset 
that, directly a member is flexed, the initial end load will be 
altered in some unknown varying manner, and the problem 
is to find the magnitude of the initial load that exists before 
any temporary alteration is induced by flexing. 


General Procedure and Further Argument. 


The general experimental procedure is as tollows :—A few 
values of either the disturbing force or the disturbing 
couple are chosen as independent variables “ X,” and the 
corresponding valnes of “ Y,” either of the deflexion or of 
the slope at a particular point, are experimentally determined, 
tablulated, and plotted, as indicated by the curve B, B, in 
fig. 3. 


Fig. 3. 


Now the ratio a established infeach of the various cases 
dealt with later, is a function of the end load P (assumed 
variable with flexure), and holds for all values of P. 

When the initial value of P=P,. required, Y=0 and 


X=0; or the initial magnitude of the ratio X=0 and is 


a definite value =“ a,” say, for each set of conditions. 
The value of “a” is given by the slope of the curve B, B, 
when extended backward to the origin, as shown by the 


broken line in fig. 3; i.e, “a”= tan KOX, and may be 
calculated from the curve B, B, or approximately deter- 
mined graphically. 

Having determined this value of “a ”—-which hereinafter 


Y // 
| 
| d=tanta 
| \ 


Maximum Stresses in Ties and Struts. 947 


will for each case treated be referred to as “the initial deriva- 
tive of the X, Y graph,” or more shortly as “the X, Y 
derivative,’ —it is substituted for x in the equation for the 
particular method used, and the equation is then solved for 
P, by approximation or by graph. (For graphical solution 
methods, see Stress Charts later.) 

For example, in the case of a tie with pinned ends flexed 
centrally through distance 6 by a central load W (as dealt 
with in case A and illustrated in fig. 4, quod vide), the general 
equation is 


Y/X=(L—tanh nL/n) /2P 


1 
and “a” =(L— tanh n,Lg)/2P.. (Tie equation Aj.) 


A X 4 B 
L L 


To facilitate its solution by Stress Chart I., this equation 
is written 


2Pa=L—tanh nL/n, 


the suffix (small a) being suppressed for simplicity. 
In the equations, P, is expressed for simplification in 
terms of 


P,/BL, 


where EI is the flexural rigidity of the member at right 
angles to the plane of flexing, and should be made the 
smaller possible value if the member be not of circular cross- 
section. 

Thus P, required=n?, BI, and during bending P=n’ KI, 


where P i P, according as the member is a tie or a strut. 


To determine Expressions for Y/X for Ties flexed by 
various Methods. 
Method 1 4 :— 
Consider a straight uniform tie A, B of length 2L, pin- 
jointed at the ends and strained by an initial end load P,. 


Kio. 4, 
A 
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When a concentrated lateral load W is applied at mid-span 
to produce a small] deflexion in a transverse direction, assume 
Pa to be increased by elastic constraints at A and B toa 
value P, and the unflexed span to be slightly decreased to 
2L, as shown in fig. 4. 

‘Taking an origin O at mid-span on the line of action of P, 
and considering the deflexion y at + from O as negative, — 
and therefore the bending moment for the curvature shown 
as of positive sign’), 


. . @ 


where the symbols have their usual meaning. 
Let P/EL=n’. 
Then d?y/da?—1?y = WL/2EI— Wa/2EI. . . (%) 
The general solution to equation (2) is 
y=Acosh ne+Bsinh ne + Wa/2P—-WL/2P. . (3) 


From equation (3) it follows that 


émax. at C= W(L— : tanh nL)/2P as measured . (4) 


or 8/W=(L— = tanh 
7 


Here 6=Y and W=X referred to in the ‘ General 


Procedure.”’ 


Thus initially, when 6=0 and W=0, equation (5) 
becomes 


0/0O=a=Lz- tanh (Tie equation Aj.) 
(where L, is the initial unflexed half-length), 


if 
va 
Equation (6) when solved for n, by trial or graph yields 
which is sought. 
Graphical equations “to equation (6) may be readily 
obtained from the wide range of “A values” given in 


C hart I. 


zZPa 


VALUES of Zn2E/a 
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To obtain « solution from this chart, the curves plotted 
from generating equation A are used—which is equation (6), 
wherein the suffix a is suppressed for simplicity. Having 
obtained experimentally the “initial derivative a,” the corre- 
sponding value of Ela is located along the A values up the 
ordinate through n?=0-05 on the abscissee and the polar 
ray drawn—preferably on tracing-paper—through the point 
to the origin. The intersection of this ray with the appro- 
priate L curve, interpolated from the vertical scale on the 
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Stress Chart for Ties A and B by deflexion. 


extreme right of the chart, locates a. point vertically above 
the value of n,? sought along the abscissee scale of A values. 
The initial direct stress (f,) is thus found, since 


PfA=n HK’, 


where K is the radius of gyration of the cross-section of the 
member, 

A simpler expression which may be used for an approxi- 
mate solution is obtained by writing equation (6) as 


8 
1600 400. 14003501200 300, 1000750 2 2 
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since tanhn,l, rapidly approaches unity for values of 
aL, > about 4 ; and this value will be easily exceeded in ties 
of fair length and cross-section under moderate stresses : 
e.g., for « tie-rod of mild steel 1 inch in diameter in which 
JSa=10,000 lb. per square inch, H=30x10 |b. per square 
inch, and L=60 inches, 


nL=4'392. 
For 3 in. tie-rod of the same length under the above stress, 
nL=5'850. 


It should perhaps be here observed that the equations 
obtained in the foregoing case, and in all the following cases 
with hinged ends, are based upon the assumption that no 
friction couples exist about these transverse axes of con- 
straint, and that results arising from this theoretical 
condition may, of course, be easily obtained in experimental 
work where friction slipping about the pins occurs, by taking 
two series of Y values: one series for increasing values of 
X, the other for decreasing values of X, so that the effect 
of friction may be neutralized by plotting the series of mean 
Y values against the one series of X values. 


Fig. 5. 
W 
A B 
P 
W 


Method 1a. Alternative :—- 


If, instead of measuring the deflexion at ©, the slope 
(dy/dz) be measured at either end, as shown at B in fig. 5 
we shall get 


dy/dxz=0(say) = W (tanh nL. sinh nL—coshnL+1)/2P 
= W(1—sech (8) 
or 
d6/dW (initially) = a = (1 — sech 
(Tie equation A ); 
Here “a” is found by treating @ as Y,and Was X, as 


described in the ‘‘ General Procedure.”” A graphical solution 


for n4, and hence for’Pa, may then be found by the aid of 
Chart IT. 
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To use this chart, the experimental value of “a” having 
been found, the polar value of 1000 Ela/L?is next located up 
the ordinate through n?L?= 20, or located along the top of the 
chart where these polar values are extended ; then the polar 
ray is drawn. The point of intersection of this ray with the 
“A curve,” projected on to the abscissa scale, gives the 
solution to x?L2, and thus to the solution of P. 

Now, as nl increases, such nL rapidly decreases from 


CuHart II. 


for dil Values of L, 


4 


Equations— 


-¥= Zn*Ela, 
B.y = ob. 


ES of Y000 Ela 


VALUES of 2n2 Ela.=2Pa. 


T 


POLAR VALU. 


20 30 40 50 60 
VALUES of n2L% 


Stress Chart for Ties A and B by slope. 


unity to a very small fraction, and, when nL>about 6, 
equation (9) may be written 

or Pa==1/2a.. . (10) 
This is clearly shown by the “ A” curve of Chart II. 


Method 1 8B :— 
As in Method 14, but with ends constrained by fixing 
couples M, M in the line of action of P, instead of being 
pin-jointed. (See fig. 6.) 
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Using the same sign conventions as in Method 14, the 
bending moment at +. from the origin O= 


M.= W(L—2)/2+P.y—M=Eld’y/dz?, . (11) 
or d?y/da® —n*y=(WL/2- M)/EI—Wa/2EI. . (12) 
The general solution to equation (12) is 


y=A cosh nz + B sinhna + Wa/2P+ M/P—WL/2P. (13) 


Fig. 6. 


From this it follows that numerically at C, remembering 
that y is negative, 


8=— y= W(L—2n. tanh nL/2)/2n7EI; . (14) 


i.é., initially. 


92 


a 


(Tie equation Bas.) 


Graphical solutions to this equation are viven by the B values 
of Chart 1. In this case the solution to n? is located along 
the abscissa scale of B values by interpolation from the 
appropriate polar ray of Ela(B) values, and the curve of L 
as described tor the corresponding Method 18. 


Fig, 7. 


As in Method 1 4, equation (15) will be found to reduce 
to a simpler approximate expression ; for when 


ngL/2> about 5, tanh nL/2 = unity, 
and the above equation becomes 
Method 18. Alternative : 


By measuring the slope 6 at L/2 from the origin, imstead 
of the deflexion 6 at the middle. (See fig. 7.) 
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From equations (11) and (13) we get 


M,=W(tanh cosh ne—sinh nz)/2n. . (17) 


=0 when e=L/2 ; 


i.é., the point of inflexion and therefore of maximum slope 
=8@ occurs at + span from either end. 
Also from equation (13) 


dy/dxz = (tanh sinh — cosh +1)=0 
== W(l—seehnL/2)/2P.. . . . . (18) 
Thus “a” here=d6/dW (initially) =(1—sech ngLg/2)/2n2EI. 
(Tie éguation By)... (19) 


Graphical solutions to equation (19) may readily be found 
from the Beurve of Chart IJ., which has been drawn to cover 
a very wide range. A solution to n?L? is obtained from this 
chart by the intersection of the ‘“B”> curve for all values 
of L, with the polar ray given by the determined value of 
“q” in a similar way to that described for Method 1a. 
Alternative. 

In passing, it may be noted that for values of nL> 
about 12, sechnL/2 is a very small fraction, and therefore 
equation (19) becomes 


a@=1/2P, or Py=1/2a asin Method la. 


It is important to observe that in all cases, only the 
the smallest possible measurable values of -Y and X are 
necessary for the determination of the XY derivative, so that 
the range of cases as regards rigidity to which these flexing 
methods are applicable is controlled by the problem of 
accurate measurement or magnification. 

The complete Thesis contains detailed descriptions of 
flexing apparatus designed by the author to test the applica- 
bility of the theoretical expressions established here, ani 
others of much more general utility established later. The 
results obtained therefrom show that remarkable agreements 
exist between the theoretical and the practical results. 


Units. 


It will be observed that throughout the series of stress 
sharts (there are twelve in the complete series) units of force 
and of length are suppressed thereon. The reason for this 
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suppression lies in the fact—as a moment’s consideration 
will show—that the charts hold for any system of consistent 
units whatsoever, and, since the quantity E is compounded 
of units of both force and length, the particular units appli- 
cable to any experimental determination will be those taken 
for BE. Thus, if E is taken in pounds per square inch, length 
values (L) must be read off in inches, load values (P) in 
pounds, stress values (fz) in pounds per sq. inch, and second 
moment of area values (I) in inches ‘4, or radius of gyration 
values (K) in inches. If H is taken in kilograms per square 
millimetre, L must be read off in millimetres, P in kilograms, 
Ja in kilograms per square millimetre, I in millimetres ‘4, 
or K in millimetres. 

In this way, and by a judicious choice of units—although 
on some of the charts only a few curves for L have been 
drawn,—these curves may be used for a number of actually 
different lengths by the expedient of altering the units of 
length to suit a chosen curve, thereby considerably extending 
the range of the chart used. 
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C. Oscillations in Low Pressure Discharge-Tubes. By 
E. W. B. Gru, B.Sc., M.A., Fellow of Merton College, 
Oxford *. 


MAY observers have noticed the occurrence of 

oscillations when direct current discharges pass in 
gases at low pressures ; but as no systematic investigation 
appears to have been made of the phenomenon the 
following experiments leading to a simple theory were 
carried out. 

The discharge-tube T in the figure, contained two small 
circular electrodes, about 1 cm. in diameter, whose distance 


apart could be varied from a fraction of a mm. to about 
1 cm. The tube contained air at a low pressure. A 
standard adjustable air-condenser C was connected 
externally across the electrodes, and a discharge was 
passed through the tube from a direct current high voltage 
dynamo. This current was maintained at any desired 
value by the insertion in one of the leads to the dynamo of a 
2-electrode valve R. The discharge current had to pass 
across the valve and could therefore be adjusted to any 
desired value by controlling the temperature of the valve 
filament. 


* Communicated by Prof. J. 8. Townsend, F.R.S, 


| 
(4) R 


946 Mr. EK. W. B. Gill on the 


In this arrangement the potential across the discharge- 
tube is that necessary for the particular current passing, 
and the remaining potential of the dynamo falls across the 
valve. It was found that up to a limiting current of the 
order of 2 or 3 milliamps this discharge, over a very large 
range of pressures, produced oscillations varying in fre- 
quency according to circumstances from a few per second 
to about 150,000 a second. 

The oscillations were detected and their wave-length 
determined by a single valve self heterodyning wireless 
receiver graduated in wave-lengths. This receiver was 
placed a few yards from the discharge-tube, which emitted 
quite enough energy to produce a loud signal. 

The waves emitted were extremely rich in harmonics, 
and care was necessary not to confuse the fundamental 
with some of the louder harmonics. 


Experimental Results. 


In what follows the fundamental wave-length A was 
found to depend on the following factors :— 


1. The pressure p of the air in the discharge-tube 
measured in mm. of mercury. 


2. The distance d between the electrodes. 
3. The current 7 flowing through the tube. 
4, The capacity C of the condenser. 


The results of varying one only of these factors keeping 
the remainder constant were as follows :— 


1, Pressure.—Oscillations were examined from p = 85 
to p=2. It is difficult to get oscillations below p=2, but 
they can be produced for much higher values than p=85. 
The general effect of a reductionin pressure is a reduction 
in the wave length but the variation is not large. 


2. Distance between the Electrodes—The wave-length A 
increases rapidly as the distance is increased: thus for 
p=25, ma. & C=110 cm. ; whend=-7 mm. A=2700 
metres ; d=1-:1 mm. A=4200; and d=1-5 mm. A=6800. 
As the wave-length gets inconveniently large for the 
larger values of d most of the measurements were made 
with small gaps. 
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3. The Current—The current magnitude has a much 
greater effect upon the wave-length than any of the 
other factors. Thus for p=25 mm. C=110 cm., d=-7 mm. 
the corresponding values of ¢ and A are given below. 


i. Xi. 
maa. 10400 5-2 
“75 6850 5:14 
i 5250 5-25 
1-25 4250 5-31 
1-5 3500 5-25 
1-75 3000 5-25 
2 2700 5-4 
2-25 2400 5-4 
2-5 2100 5-25 


The last column gives the product of 7 and X and 
shows that very approximately 7A=constant ; for small 
values of 7 the oscillations are of audible frequency. A 
large number of sets of observations were made for various 
values of C and d, and in all cases 1A was approximately 
constant in each set. 


4. The Capacity—As the capacity is increased the 
wave-length increases and the curve giving the relation 
between C and A is very accurately a straight line. 
The line does not pass through the wee and is given 
by the equation 


A=k (C+e’), 


where & and c’ are constants ; c’ is small and depends on 
the current, being larger for the larger currents. 

For example, with p=6-7 d=1-1 mm.; when the 
current is -21 m.a. the wave-length increases uniformly, 
being 8200 metres for C=15 cm. and 12,600 metres for 
C=40 cm.; here c’ =30 cm. 

For a current of 1:35 m.a. XA increases from 2500 
metres for C=15 to 6000 metres for C=100: in this case 
c’ =45. 

The final experimental result is therefore that 

k(C+e) 
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k depends only on d and p and increases if either is 
increased. c’ is small and depends chiefly on 7 and 
decreases if 7 is increased. 

The following theory appears to be adequate to explain 
the experimental results :—Suppose that the valve fila- 
ment is first of all heated to the temperature necessary to 
allow a current 7 to pass the valve. At the moment when 
the external E.M.F. is applied practically all of it falls 
across the valve owing to its very small capacity. The 
current 7 thus commences to flow instantaneously through 
the valve, and, provided the external E.M.F. is large 
enough, this current 7 remains constant whatever happens 
in the rest of the system. 

At the outset the tube insulates, the potential across it 
being low, and all the current 7 flows into the condenser, 
and the potential across it and the tube increases until 
the sparking potential is reached. At this stage the tube 
begins to conduct and the current through it begins to 
grow from zero value ; but as the constant current 7 is still 
flowing into the combination a portion of it will still flow 
into the condenser until the discharge current has grown 
to the value 7. At the instant, therefore, when current 
ceases to flow into the condenser the potential across it 
and the discharge-tube 2s slightly above the sparking 
potential, say vy. 

For air this state of affairs is essentially unstable ; as 
the sparking potential is higher than that required to 
maintain a current, the potential cannot remain at vy, 
the condenser begins to discharge through the tube, and 
the current through the tube begins to exceed 7. 

The relation between E.M.F. and current in a discharge- 
tube is that for zero current the E.M.F. is the sparking 
potential, and that as the current is increased the E.M.F. 
falls till it reaches a certain value v, corresponding to a 
current 7%, after which, however much the current is 
increased, the E.M.F. remains practically the same. 

The superposition of the condenser discharge current 
upon the current 7 flowing through the tube runs the 
potential across the tube down to this value vy, and the 
condenser then ceases to discharge, as the P.D across it 
can fall no further. If 2 is less than the critical value 7, 
the discharge must go out, as the potential v, at this 
instant across the tube will only maintain currents equal 
to or greater than 7,. If 7 is equal to or larger than 7, the 
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discharge current remains constant, the condenser has no 
further effect, and no oscillations occur, which explains 
why no oscillations could be found unless the current was 
small. 

For the small currents the discharge goes out and the 
current again all flows into the condenser, charging it from 
v, to v1, at which point it again discharges, after which the 
process repeats itself indefinitely. The time for the 
charging is 

a 
and if the time for discharge is t) the total time for one 


cycle is eee +t). But Ais proportional to the total 


time of one cycle, and since Az was found to be constant if C 
is fixed, this shows, as might have been expected, that f, 
is indefinitely small; for vy, and v,. are constant and 
independent of 7. Thus far, then, it has been proved that 


It remains to account for the constant c’. 

At first sight it might appear that stray capacities and 
the capacity of the electrodes might account for this, but 
these corrections cannot be more than a few electrostatic 
units while c’ is about 30. 

It must be remembered, however, that when the 
potential across the condenser reaches its maximum value 
v, the gap is discharging, and its capacity is not solely due 
to the electrodes but is very much increased by the space 
charges. That this is the correct explanation is borne 
out by the fact that c’ increases if 7 is increased. 

Finally, therefore, 


Vv, — V2 is the difference between the sparking potential and 
the minimum maintenance potential. This difference is 
increased if we increase the distance between the electrodes 
or raise the pressure, which is in accordance with the 
observed facts. 

A simple experiment verified the deduction that ¢, is 
very small. A hot wire ammeter was placed to read the 
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current flowing into and out of the condenser. There is no 
possibility of resonance on the wave-lengths emitted, and 
the condenser current cannot therefore exceed the current 
flowing into the system. But when 1 m.a. was flowing 
through the 2-electrode valve the hot wire ammeter 
read as much as 30 m.a. This can only occur if the dis- 
charge current reaches momentarily a very high value 
causing the AC instrument to give a reading greater than 
the DC instrument. 

The impurity of the emitted wave is evident from a 
consideration of the process occurring. The potential 
across the condenser rises uniformly from v, to vy, and falls 
instantaneously to vy, each cycle. 

An application of Fourier’s analysis indicates that in 


such a case the amplitude of the nth harmonic is : of the 
n 


amplitude of the fundamental. 

The presence of the external condenser is not essential 
for the production of oscillations, the capacity of the 
electrodes magnified, as explained by the presence of the 
space charges, being sufficient to produce oscillations. 

It appears, therefore, that any discharge-tube is in a 
state of. oscillation provided that the sparking potential 
is higher than the maintenance potential and that the 
current is less than that at which the maintenance potential 
is independent of the current. At the higher pressures 
there is a fair range of currents over which oscillations 
occur, but at pressures of 2 or 3 multimetres they only 
occur if the current is a very small fraction of a milli- 
ampere. 

It is scarcely necessary to add that the existence of 
oscillations in no way depends on the particular method of 
limiting the current by a valve. If the current is con- 
trolled instead by a large resistance oscillations occur 
equally well, but the calculations are not quite so simple. 


Professor Townsend, in whose laboratory the experi- 
ments were made, has most kindly given me much advice 
and assistance. 


(The Editors do not hola themselves responsible for the 
views expressed by their correspondents. | 
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